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In this dissertation, we use the following notation: 


R, € 

R n , C" 

x m £n x m 

N", P" 

S n , E 
In- ® 

A r , A*, A~ } 
tr A, 71(A) 
He G 


Sh G 


real numbers, complex numbers 

7? x 1 real column vectors, n x 1 complex column vectors 
77 x m real matrices, nxm complex matrices 

77 X 77 nonnegative definite matrices, nxn positive definite matrices 
77 x 77. symmetric matrices, expectation 
77 x 77 identity, Kronecker product 

transpose of A, complex conjugate transpose of A, inverse of A 

trace of A, range space of A 

Hermitian part of arbitrary complex matrix G 

He G= l 2 (G + G*) 

Skew-Hermitian part of arbitrary complex matrix G 

Sh G = \(G - G*) 


For convenience, we define the following terms: 

£ = BRi l B T , £ = C t V 2 -'C. 


We also define the dimensions of the various signals below: 


x(t) € M n is the plant state vector 
x c (t) G R n,: is the controller state vector 
u(t) € R m is the control input signal 

w(t) G R rf is a unit-intensity, zero-mean, Gaussian white noise signal 

d(t) G R s is the uncertainty input signal 

y{t) G R' is the measurement output signal 

z(t) G R p is the performance output signal 

e(t) G R r is the uncertainty output signal 


Furthermore, we define the order of the closed-loop system to be n = n + n c , where 
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77 is the order of the plant and n c is the order of the compensator. We also define the 
following closed-loop matrices: 



A BC c 
B c C A c + B c DC c 



B o 
B C F X 


L 


C 0 = [ Co F 2 C c ] , E = [ E x E 2 C c ] . 


Finally, we define 


R = E r E = 


Ri 0 
0 CJR 2 C c \ ’ 


where R\ = EjE\, R 2 = EjE 2 , and Ri 2 = EjE 2 = 0, and 


r = e t e = 


Vx 0 

o B c v 2 Bj j ’ 


where Ij = D x Dj, V 2 = D 2 Dj , and V n = D x Dj = 0. 


xviii 


A D\ 
B c D 2 


Summary 


The ability to develop an integrated control system design methodology for ro- 
bust high performance controllers satisfying multiple design criteria and real world 
hardware constraints constitutes a challenging task. The increasingly stringent per- 
formance specifications required for controlling such systems necessitates a trade-off 
between controller complexity and robustness. The principle challenge of the minimal 
complexity robust control design is to arrive at a tractable control design formulation 
in spite of the extreme complexity of such systems. Hence, design of minimal com- 
plexity robust controllers for systems in the face of modeling errors has been a major 
preoccupation of system and control theorists and practitioners for the past several 
decades. 

Although the theory for designing linear output feedback controllers is quite ma- 
ture, the actual solution of the design equations can be a daunting task. This becomes 
even more difficult if an optimal reduced-order controller is sought. In this disser- 
tation, we develop a general fixed-structure control design framework that addresses 
the following paradigm: Reduce control law complexity subject to the achievement 
of a specified accuracy in the face of a specified level of uncertainty. 

The control law complexity is reduced by developing a decentralized static out- 
put feedback formulation for fixed-structure controller synthesis. The decentralized 
static output feedback formulation captures a large class of controller architectures 
within a common framework and allows a common numerical algorithm to be used 


for computational purposes. 

We first discuss this decentralized static output feedback framework and demon- 
strate its applicability on the space-based ACTEX control testbed. Since the strict 
architecture constraints of the ACTEX flight experiment preclude standard LQG and 
'Hoc techniques, we show that the decentralized static output feedback formulation 
for fixed-structure controller synthesis can directly account for the controller archi- 
tecture constraints and improve closed-loop performance by designing fixed-structure 
T^-optimal controllers. 

The next problem we consider is fixed-structure stable "^-optimal controller syn- 
thesis using a multiobjective optimization technique. The problem is presented in the 
decentralized static output feedback framework developed for fixed-structure dynamic 
controller synthesis. A quasi-Newton/continuation algorithm is used to compute so- 
lutions to the necessary conditions. To demonstrate the approach, two numerical 
examples are considered. The first example is a second-order spring-mass-damper 
system and the second example is a fourth-order two-mass system, both of which are 
considered in the stable stabilization literature. The results are then compared with 
other methods of stable compensator synthesis. 

Next we use a similar approach to consider fixed-structure % 2 -optimal relative 
degree two controller synthesis. By considering dynamic controllers structured to 
appear as the augmentation of two strictly proper dynamic controllers in series, the 
relative degree condition is guaranteed. Three examples are presented to demonstrate 
the effectiveness of this design technique. 

We then explore the applicability of the implicit small gain guaranteed cost bound 
for controller synthesis. For flexibility in controller synthesis, we adopt the approach 
of fixed-structure controller design which allows consideration of arbitrary controller 
structures, including order, internal structure, and decentralization. Two numerical 


xx 


examples that have been addressed by means of alternative guaranteed cost bounds 
are presented to demonstrate the fixed-structure/implicit. small gain approach to ro- 
bust controller synthesis. 

One of the fundamental problems in feedback control design is the ability of the 
control system to maintain stability and performance in the face of system uncer- 
tainties. To this end, elegant multivariable robust control design frameworks such as 
Hoc control, C] control, and p-synthesis have been developed to address the robust 
stability and performance control problem. An implicit assumption inherent in these 
design frameworks is that the controller will be implemented exactly. In a recent 
paper by Keel and Bhattacharyya, it was shown that even though such frameworks 
are robust with respect to system uncertainty, they are extremely sensitive, or frag- 
ile, with respect to errors in the controller coefficients. Here we extend the robust 
fixed-structure controller synthesis approach to develop controllers which are robust 
to system uncertainties and non-fragile, or resilient , to controller gain variations. 

Next, we develop linear, fixed-order pressure rise feedback dynamic compensators 
for axial flow compressors. Unlike the nonlinear static controllers proposed in the 
literature possessing gain at all frequencies, the proposed dynamic compensators ex- 
plicitly account for compressor performance versus sensor accuracy, compressor per- 
formance versus processor throughput, and compressor performance versus distur- 
bance rejection. Furthermore, the proposed controller is predicated on only pressure 
rise measurements, providing a considerable simplification in the sensing architecture 
over the bifurcation-based and backstepping controllers proposed in the literature. 

Finally, we use a unifying absolute stability result for mixed uncertainty to obtain 
fixed-structure controllers and fixed-order stability multipliers which provide robust 
stability and performance. The robust controller synthesis technique proposed here 
permits the treatment of fully populated real uncertain blocks which may, in addition, 


possess internal structure. The ability to address real uncertain blocks is based on the 
use of an appropriate class of multipliers whose structure is compatible with the real 
block uncertainty. Hence, tailoring the multipliers to the structure of the uncertainty 
not only leads to the ability to address more general uncertainty characterizations but 
can also lead to less conservative controllers than obtained from standard mixed-p 
synthesis. 
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CHAPTER 1 


Introduction 

1.1. Robust Fixed-Structure Control Design: 
Motivation and Overview 

The growing complexity of dynamic systems results in unavoidable discrepan- 
cies between real physical systems and the mathematical models used to describe 
them. These uncertainties , in turn, result in severe degradation of control system 
performance. Thus, one of the main objectives of feedback control system theory 
is to design controllers that are robust with respect to system uncertainties as well 
as guarantee specific performance objectives. The problem of robust control design 
constitutes a significant challenge in mathematical system theory which, at the same 
time, addresses a fundamental issue in the practical implementation of feedback con- 
trol systems, namely, modeling uncertainty. Modeling uncertainty must be character- 
ized and quantified so that it can be accounted for within the control design process. 
For example, the dynamics of large flexible space structures are highly nonlinear (due 
to geometric and material nonlinearities) but are commonly approximated by lin- 
ear models. Furthermore, since flexible structures are inherently infinite dimensional 
systems, any finite dimensional approximation model will exhibit significant error, 
particularly as the modal frequency increases. 



These uncertainties arising due to inexact modeling are referred to as plant uncer- 
tainties and are broadly classified into parametric and nonparametric uncertainties. 
Parametric uncertainty here describes errors that can be translated into uncertainty in 
the physical elements of some time-invariant state-space representation of the design 
model (i.e., perturbations or uncertainties in specific parameters of the physical sys- 
tem). On the other hand, nonparametric uncertainty is best viewed in the frequency 
domain and describes errors that have bounded gain but arbitrary phase (e.g., uncer- 
tainties due to unmodeled system dynamics or system linearization). The distinction 
of parametric versus nonparametric uncertainty in the plant model is of paramount 
importance in robust control design. For example, in the problem of vibration control 
of flexible space structures, if the stiffness matrix uncertainty is modeled as nonpara- 
metric uncertainty, then perturbations to the damping matrix will inadvertently be 
allowed. Consequently, stability and performance predictions for a given compensator 
will be extremely conservative, which limits achievable performance. 

Hence, robust control theory mainly deals with two issues; namely, the qualitative 
issue of robust stability and the quantitative issue of robust performance. Robust 
stability addresses the stability of a given system in the face of uncertainties while 
robust performance addresses performance degradation due to system uncertainties 
over the region of robust stability. Often worst case performance of the system in the 
face of all possible uncertainties is addressed as a measure of robust performance. 

Modern multivariable feedback control theory and application has been one of the 
most rapidly growing areas in the scientific and engineering communities for the past 
several decades. Some of the most fundamental advances in this field can be traced 
back to World War II. During the post- World War II era, the emergent superpowers 
turned their research focus to aerospace technology in order to compete in the race to 
space. One of the central achievements of this research effort was the development of 
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state-space based methods introduced in the West by R. E. Kalman in the 1960’s [71, 
72] which led to the Kalman filter [68,73]. Another pivotal achievement of the early 
1960’s was the development of the Linear Quadratic Regulator (LQR) [4,67,80,101] 
which, when combined with the Kalman filter, led to the formulation of the Linear 
Quadratic Gaussian (LQG) controller design methodology [4,6,67,80,107], These 
revolutionary breakthroughs in optimal navigation and control led to the successful 
launch of the Apollo mission which resulted in the first manned moon landing on July 
20, 1969. 

One of the most attractive features of the LQG controller design methodology 
is the characterization of the optimal compensator gains via a system of two de- 
coupled Riccati equations. A number of commercially available computational tools 
offer efficient solutions to standard Riccati equations. Unfortunately, however, the 
dimension of an LQG controller is always equal to the dimension of the design plant. 
Since on-board processors have limited throughput and uncertainty is always present 
in the system, this necessitates robust minimal complexity controllers which further 
renders LQG controllers impractical for many applications. This has motivated the 
study of optimal fixed-order controllers. Balanced truncation techniques [124] were 
used to reduce the order of the plant or the designed full-order controller, however 
these reduced-order controllers could not even guarantee closed-loop stability, let 
alone performance. Another approach, developed to directly synthesize fixed-order 
controllers, is the fixed-structure control framework developed in [62,63]. In this 
approach, the compensator structure is fixed a priori and the optimization is per- 
formed over the compensator parameters. The application of fixed-structure control 
theory to the fixed-order controller design problem yields a characterization of the 
optimal fixed-order controllers via a set of coupled Riccati and Lyapunov equations, 
each containing a projection matrix which motivated the name “optimal projection 


equations” in [62,63]. 

At the same time as the development of the LQG controller, the work on absolute 
stability theory which forms the basis of modern robust control theory was being pio- 
neered by the Romanian mathematician V. M. Popov [102]. Absolute stability theory 
addresses the stability of feedback systems whose forward path contains a dynamic 
linear time-invariant system and whose feedback path contains a memoryless (pos- 
sibly time- varying) nonlinearity (uncertainty). These stability criteria are generally 
stated in terms of the linear system and apply to every element of a specified class 
of nonlinearities. Research efforts in this direction were vigorously pursued in the 
former Soviet Union by V. A. Yakubovich [116-120] and Y. Z. Tsypkin [109]. In the 
United States, progress along the same direction was made by R. E. Kalman [69], 
E. I. Jury [65,66], R. W. Brockett [19,20], J. L. Willems [20], and K. S. Naren- 
dra [91-94]. The significant progress made towards resolving the absolute stability 
problem is now well documented in research monographs, such as Aizerman and Gant- 
macher [1], Lefschetz [82], and Popov [103]. A more modern treatment of the subject 
is given by Safonov [104] while excellent book treatments are presented in Narendra 
and Taylor [94] and Vidyasagar [111]. 

Since all real-world systems exhibit nonlinear behavior and possess numerous un- 
certainties (due to such phenomena as exogenous disturbances, noise, mathemati- 
cal modeling errors, unmodeled nonlinear dynamics, unmodeled high frequency dy- 
namics, and unmodeled actuator and sensor dynamics), a viable controller synthesis 
methodology must be able to account for system nonlinearities and uncertainties. 
The problem of accurately controlling system performance variables whose dynam- 
ics contain significant nonlinearities and uncertainties poses a challenging problem in 
control system design. The critical issues of system uncertainty can be traced back 
to H. S. Black’s 1927 patent where large loop gains were proposed for addressing 
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the problem of uncertainty [17,28]. In classical single-input, single-output control 
design techniques based on the Nyquist criterion and Bode gain and phase plots, the 
question of robustness to system uncertainty was addressed by requiring that the 
overall closed-loop system possess prespecified relative stability margins, i.e., gain 
and phase margins. The issue of parametric plant uncertainty, however, was largely 
neglected. Even though, as demonstrated by Kalman [70], optimal linear quadratic 
state-feedback regulators possess guaranteed robustness properties, i.e., oo dB up- 
ward and 6 dB downward gain margins and ±60° phase margin, these robustness 
properties are misleading since the guaranteed gain and phase margins are valid only 
for uncertainty at the system input. Furthermore, it was shown in [55] that the ro- 
bustness properties of LQR designs could diminish with increasing parametric system 
uncertainty. Finally, these robustness guarantees are nonexistent in observer-based 
LQG controllers [4,29]. 

In the late 1970’s and early 1980’s, a renewed research interest emerged in ad- 
dressing plant uncertainty which led to an accelerated progress in the field of robust 
control (see [28] for a representative collection of papers on this topic). A main mo- 
tivation for this growing research effort appears to be the ever increasing complexity 
of the modern systems within the engineering field (such as large flexible structures, 
advanced tactical fighter aircraft, and variable-cycle gas turbine engines, to name a 
few), as well as complex economic and biological systems. Thus the predominant con- 
siderations in control law design for modern engineering systems have become control 
law complexity and control law robustness, respectively. Indeed, with increasing sys- 
tem complexity comes increasing (and usually overriding concern with) system cost, 
reliability, and maintainability, whereas with increasing accuracy requirements come 
increasingly complex control systems. Since, generally speaking, the more complex 
the control system, the more it costs, the less reliable it is, and the harder it is 


to maintain, it follows that high accuracy requirements conflict with control system 
complexity requirements in highly complex systems. In fact, they are also in conflict 
with each other through the specification of control law robustness. In an attempt to 
capture robustness guarantees within observer based designs, Doyle and Stein [32,33] 
proposed a two-stage analysis/synthesis recovery framework that led to the LQG/LTR 
methodology [108] yielding dynamic compensators with recovered LQR-type margins. 
However, the LQG/LTR methodology is limited to the recovery of LQR margins and 
employs high-gain feedback. In addition, this design framework is not applicable to 
plants with non-minimum phase zeros and involves inversion of the plant which leads 
to controllers which are extremely sensitive to parametric uncertainty. 

One approach to controller synthesis in the presence of unstructured uncertainty 
is the well known 'H (X > methodology. Several authors have shown that the prob- 
lem can be solved via a pair of modified Riccati equations [11,31,45,99]. However, 
as is also well known, this methodology is highly conservative for parametric uncer- 
tainty. Using a Lyapunov bounding framework, the authors in [12,13,88] addressed 
the problem of robust analysis and synthesis in the face of parametric uncertainty. 
Several bounding functions were considered in [12], while [10] extended the quadratic 
bounding technique of [100] to robust fixed-order controller synthesis. Alternative 
approaches to robust stability and performance in the presence of parametric uncer- 
tainty are discussed in [28,77,100]. However, a major drawback of the conventional 
Lyapunov bounding techniques [48] is their inability to restrict time-variation of the 
parametric perturbations. In a recent series of papers [46,48-50,52], a refined Lya- 
punov function technique was developed to overcome some of the current limitations 
of Lyapunov function theory for the problem of robust stability and performance in 
the presence of constant real parameter uncertainty. The authors in [46,48-50,52] ex- 
tend the rich theory of absolute stability to develop parameter-dependent Lyapunov 
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functions that guarantee robust stability and performance in the presence of constant 
real parametric uncertainty. Finally, using frequency domain stability criteria, the 
authors in [57, 61] demonstrate that the seemingly remote modern robust stability 
and performance tools based on the structured singular value [30, 123] can be directly 
connected to absolute stability theory via selection of stability multipliers for various 
classes of nonlinearities. Hence, [46,48-50,52,57] provide an alternative approach to 
mixed -\i synthesis [122], while avoiding the standard D, N-K iteration and curve 
fitting procedure. 

Prior experience with the state-space techniques and the LQG controllers prompt- 
ed many researchers to seek state-space-based robust controller synthesis methods 
that essentially try to mimic the solution features of the LQG controller, i.e., controller 
gains based on solutions to modified Riccati equations. Recent advances along these 
lines are given in [10,11,13,31,45,46,52,57,76,88,99,100]. Furthermore, the optimal 
projection equation approach of Hyland and Bernstein [62, 63] has been extended to 
the problem of robust controller synthesis in the presence of unstructured system 
uncertainty [11,45,88], arbitrary time-varying parameter uncertainty [10,12,13,88], 
and constant real parameter uncertainty [46,48-50,52,57] via the quadratic Lyapunov 
bounding framework developed in [44]. 

1.2. Brief Outline of the Report 

In this research, we will build on the results of [14,95] in several directions. In 
Chapter 2, we formulate the decentralized static output feedback problem for fixed- 
structure controller synthesis. This formulation captures a large class of controller 
architectures within a common framework and allows a common numerical algorithm 
to be used for computational purposes. This framework will provide the basis for 


much of the subsequent controller synthesis, such as the stable 7i 2 -optimal controller 
synthesis discussed in Chapter 3. In this chapter, a multiobjective optimization tech- 
nique is used to optimize the % 2 cost of the closed-loop system while maintaining 
controller stability. The decentralized static output feedback framework is also ex- 
ploited in Chapter 4 when we consider "H 2 -optimal relative degree two controller 
synthesis. 

Next, in Chapter 5, we consider the applicability of the implicit small gain guaran- 
teed cost bound for controller synthesis. In Chapter 6, the decentralized static output 
feedback framework is used to develop controllers with the ability to maintain control 
system stability and performance in the face of system uncertainties, which could 
even include variations in the controller gains. Then, in Chapter 7, we use the fixed- 
architecture control methodology to develop linear, fixed-order pressure rise feedback 
dynamic compensators for axial flow compressors. Unlike the nonlinear static con- 
trollers proposed in the literature possessing gain at all frequencies, the proposed 
dynamic compensators explicitly account for compressor performance and is pred- 
icated on only pressure rise measurements, providing a considerable simplification 
in the sensing architecture over the bifurcation-based and backstepping controllers 
proposed in the literature. 

Finally, in Chapter 8, we use a unifying absolute stability result for mixed uncer- 
tainty to obtain fixed-structure controllers and fixed-order stability multipliers which 
provide robust stability and performance. The robust controller synthesis technique 
proposed here permits the treatment of fully populated real uncertain blocks which 
may, in addition, possess internal structure. The ability to address real uncertain 
blocks is based on the use of an appropriate class of multipliers whose structure is 
compatible with the real block uncertainty. Hence, tailoring the multipliers to the 
structure of the uncertainty not only leads to the ability to address more general 
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uncertainty characterizations but can also lead to less conservative controllers than 
w 

obtained from the standard mixed-/i sythesis techniques. 

Finally, in Chapter 9, conclusions and recommendations for future research are 
discussed. 
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CHAPTER 2 


Fixed- Structure Control 
Framework 

2.1. Introduction 

This chapter reviews the decentralized static output feedback problem formula- 
tion for fixed-structure controller synthesis. As shown in [14,95], the decentralized 
static output feedback format captures a large class of controller architectures within 
a common framework and allows a common numerical algorithm to be used for com- 
putational purposes. 

Consider the (m + q + l)-vector-input, (m + q + 1)- vector-output decentralized 
system shown in Figure 2.1, where w is the exogenous disturbance input, z is the per- 
formance variable, the signals yi and it*, i = 1, . . . , m, are measurement and control 
signals, respectively, and e* and di , i = 1 ,...,q, are used to account for model un- 
certainty. The decentralized static output feedback multi vector-input, multi vector- 
output system shown in Figure 2.1 is characterized by the dynamics 

^ ; - T " 5 1 ■ » ■ i-. ! ■: 

m Q 

x(t) = Ax(t) «,(«)+£ B dt d k (t)+B w w{t), i€[0,oo), (2.1) 

j= 1 k= 1 

m q 

Ui(i) = ] D yUij Uj(t) ^ T^yd ik dkif) +T3ywi w {t)<> z = 1, 2, . . , , m, (2.2) 



Figure 2.1: Decentralized Static Output Feedback Framework 

m q 

ei{t ) = C ei x{t)+^2T> eUlj iij{t)+^2v edik d k (t)+V eWi w(t), 1 = 1,2,...,?, (2.3) 

j = 1 k = 1 

m 9 

z(t) = (2-4) 

i=i fc=i 

In the above formulation, model uncertainty is represented by the decentralized 
static output feedback map 

di(t) = Aiei(t), i = l,...,g, (2.5) 

where the uncertain matrices A, are not necessarily distinct. To represent decentral- 
ized static output feedback control with possibly repeated gains, we consider 

Ui{t) = JCiyi{t), * = l,...,m, (2.6) 

where the matrices Ki are not necessarily distinct. Reordering the variables in (2.5) 
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and (2.6) if necessary and defining 



Ul(t) 


y i{t) 

u(t) = 

\ 

» y(t) = 

j 


^m(^) 





di{t) 


ei (t) 

II 

; 

II 

•Vi 

* 


dq(i) ^ 


CqW. 


(2.7) 


0„ = [B u , ••• 0 U „], 0< = [0* ■■■ B d ,],V„=[D m ■■■ 


<11 


D = 

j L'yu — 

’ T> 

^ yuw 

V 1 

, 2>y d = 

^X/dll 

' ’ ' 'Dydi,, 


Cym _ 


V 

yuml 

V 

^ y^-mm _ 


^Dydmi 

'Dyd mq _ 



1 

V ... 

1 

, V eu = 

ts eu\\ 

. . . T> 

^euun 

,T> ed = 

1 

... 

CL 

T^ed iq 


C e 

L e< ? J 


V 

_ £y eu q i 

* ’ ’ T^eUqm m 



T^edqq 


V zd = [ v zdl • 

T^zdq ] 7 

V — 

ISyiv ~ 

V 

*Sywi 

V = 

5 ^ew — 

" V 

^ew i 




V 

LSyWm 


V 

t^eWq 


(2.5) and (2.6) can be rewritten 


(2,8) 

(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 


d(t) = A e(t), 
u{t) = ICy(t), 


( 2 . 12 ) 

(2.13) 


where A and K have the form 

A = block-diag [/^, ® Ai, . . . , I^ v <S> A„] , (2-14) 

K = block-diag [/*«£!,..., h g ® /Cj , (2.15) 


where u is the number of distinct uncertainties A* (E O’* x ^‘ or R p< x ^‘, ^ is the number 
of repetitions of uncertainty A*, g is the number of distinct gains G R r ‘ xc ' and 
4>i is the number of repetitions of gain /Cj. Note that /Ci, . . . , JC g are not necessarily 

U 9 

square matrices, and ipi = q and <j>i = m. Furthermore, define the matrices 



Quj and Qn.j t0 be 


0ri</>i xrj 


Oci^iXCj 

Or 2 02 XT, 


0c 2 </>2XCj 

Orj_i0i_ixr, 


XCi 

Or^j-llxri 

Ir t 

> = 

4 

b r,(0, -j)xr t 



9ri +1 </> t+ i XT* 


0« + 1 * +l XCi 

xrj 


m XC( . 


(2.16) 


where i = 1, . . . , v, and j = 1, . . . , <fo. 

With the definitions in (2.7)— (2.11), the transfer function G(s ) from [u T ,d T ,w T ] T 
to [y T , e T , z T ] r of the decentralized system has the realization 


G(s) ~ 


A 

B u \ B d ' B w ] 

_py_ _ 


QJ 1 

T^eu j T^ed , 


T^ZU \ T2 zd \ T^zw . 


(2.17) 


which represents the linear, time-invariant dynamic system 


x{t) = Ax(t) + B u u(t) + B d d(t) + B w w(t), t € [0, oo), (2-18) 

y(t ) = C y x(t) + V yu u(t ) + T> yd d(t) + T> yw w(t) y (2-19) 

e(t) = C e x(t) + V eu u(t) + V ed d{t ) + V ew w(t), (2.20) 

z(t) = C z x(t) + V zu u(t) + V zd d(t ) + V zw w(t), (2.21) 


which is equivalent to (2.1)— (2.4). Furthermore, by rewriting the decentralized control 
signals (2.6) in the compact form given by (2.13), the closed-loop system realization 
from [d T , io T ] T to [e T , z T ] T is given by 



' A 

L4 

' D 

G(s ) ~ 

Co 

jh\ A 


. E 

Ex 

: e 0 , 


( 2 . 22 ) 
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where 


A = A + B u KL- K l C y , B 0 = B d + B u ICL^V yd , D = B W + B u fCL^V yw , 

Co — C e + V eu K.L K l C y , Dq = T> cd 4 - D CU KL^ Dy d , D\ = D ew + D eu JCL ^ D yw , 

E = C z + D zu K.Lf^C y ^ E\ — T> zd + T) zu JCL K l D yd , Eq = D zw + D zu fCL^ D yw , 

and where L>c = I - V y JZ. Note that we assume det(Ljc) 7 A 0 for all /C given by 
(2.15) to ensure the well-posedness of the feedback interconnection. 

2.2. Performance Constraints 

The decentralized static output feedback framework can be used to help synthesize 
controllers that are optimal with respect to user-chosen performance criteria. A 
performance criterion consists of a cost function and one or more constraints. The 
cost function represents some characteristic of the controlled system which a desirable 
controller will minimize, while the constraints represent properties that an acceptable 
controller must have. An example of a cost function is an induced norm of the closed- 
loop transfer function, while examples of constraints include asymptotic stability of 
the nominal closed-loop system or robust stability with respect to uncertainties of a 
certain size and/or structure. 

As an example, consider the T^-optimal centralized strictly proper dynamic con- 
trol (Linear Quadratic Gaussian) problem. Here we will consider the n th -order stabi- 
lizable and detectable plant with noisy measurements, uncertainty signals, and per- 
formance variables given by 


x(t) = Ax(t) -I- Bu{t) + B 0 d(t ) + Diw(t), t 6 [0, 00), 

y(t) = Cx(t) + Du(t) + F\d(t) + .D 2 tu(f), 
e(t) = C 0 x(t) + F 2 u(t), 


(2.23) 

(2.24) 

(2.25) 


z(t ) = E\x{t) + E 2 u{t). 


(2.26) 


To recast this system into the decentralized static output feedback format, we first 
define the decentralized control signals to be 


ui(t) = A c yi(t), u 2 (t ) = B c y 2 {t), u 3 (t) = C c y 3 (t). 


(2.27) 


Note that we can also write the control signals (2.27) in the compact form given by 
(2.13), where, for the centralized strictly proper dynamic control problem, K takes 
the form 

JA 0 0 | 

(2.28) 




A c 0 0 

0 0 
0 0 C c 


Further note that Q Li] and Q r 0 are given by 



1 

o 

1 


0n c xn c 


0n c xm 

II 

►J 

o* 

OflcXTlc 

, Ql 2 i — 

In c 

> Qui - 

0n c xm 


0 mX n c 


0mxn c 


1 m 


Qrh — [ Iji c Ofi c x( O nc xn c ] > Qr 2 i [ ^ixn c Ojxn c ] > Qr. 31 [ ^n c x/ fn c ] • 

With these definitions, the closed-loop system is given by 

3 

i(t) = Ax(t) + ^2 #u , u i(t) + B w w(t), t £ [0, oo), 

1=1 

3 


Vi(t) — B yi x(^) T 'y ^ By Un Uj (t) + 'Dy Wi w[t ), i 1,2,3, 
j = l 

3 

z(t) = C z x(t) + 


i = 1 


where 
A = 


i 1 

O O 

^ o 

II 

' 0 ‘ 
I 

, B U2 = 

r i 

O 

1 1 

II 

i 1 

CQ O 

1 1 

II 

1 1 

Q o 

1 l 


By\ — [ 0 I ] > Cyi ~ [ & 0 ] > B ya — [ 0 I ] , 

By Uu = 0, By Uli = 0, 'Dy Ul3 = 0, By U2l = 0, ^J/tl22 
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Ip* 








Pyu 31 — 0, V yU32 — 0, P yUi3 — 0, Pyw 1 — 0, Pyw 2 — B^i Pyw 3 0, 

C, = [ J5i 0 ] , P JU1 = 0, V zu , = 0, V ZUl = E 2 . 

The LQG problem can now be defined as 

min ||G 2tu (s )||2 , (2.29) 

flea's 

where /C s is the set of all 1C of the form (2.28) such that A is stable. 

If K £ /C s , then 

\\G zw (s)\\l = tvQR, (2.30) 

where Q is the unique, nonnegative definite solution to the Lyapunov equation 

0 = AQ + QA t + V. (2.31) 

Necessary conditions for optimality involve a Lagrangian function that accounts for 
the constraint (2.31). The Lagrangian for the ^-optimal control problem is given by 

C(P, Q, Ac, B c , C c ) = trQR + trP \aQ + QA T + v] , (2.32) 

where PeN" is a Lagrange multiplier, and the partial derivatives of C(P,Q,A c ,B c ,C c ) 
are given by 

= AQ + QA t + V, (2.33) 

dP 

^ = A r P + PA + R , (2.34) 

dQ 

H = 2 Ql tl r [B*PQCj + BlPDVj w +V] u EQCj] L?Ql.., 

i= 1,2,3, j = l, (2.35) 

where r = I + Pj u L^ T /C T . 

Finally, although the LQG problem does not account for model uncertainty, for 
completeness we note that, for centralized strictly proper dynamic control problems 
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where the model uncertainty for the nominal system matrices (A, B, C ,D) is modeled 


as 


AA = B 0 AC 0 , AB = B 0 AF 2 , AC = FiACq, AD = F X AF 2 , 


the matrices B d , C e , V eu , and V ed have the form 


B d = 


B o 
0 


j c e = [ Co 0 ] , V eu — [ 0 0 F 2 ] , T> yd = 


0 

Fi 

0 


2.3. Quasi-Newton Gradient Optimization 


Once the problem has been posed in terms of the gradients of an associated La- 
grangian, a general-purpose BFGS quasi-Newton algorithm is used to solve the non- 
linear optimization problem. The line-search portions of the algorithm were modified 
to include a constraint-checking subroutine which decreases the length of the search 
direction vector until it lies entirely within the allowable set of parameters that yield 
a stable closed-loop system. This modification ensures that the cost function remains 
defined at every point in the line-search process. Numerical experience indicates that 
this subroutine is usually invoked only during the first few iterations of a synthesis 
problem. For details of the algorithm, see [35]. 

One requirement of gradient-based optimization algorithms is an initial stabilizing 
design. Here, this was usually accomplished by using either the LQG controller for 
full-order designs or a balanced truncation of the LQG controller to obtain a reduced- 
order design. 
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2.4. Structural Dynamics Modeling of the ACTEX 
Flight Experiment 

To demonstrate the applicability of the decentralized static output feedback for- 
mat for controller design, we present the ACTEX flight experiment. This testbed 
provides a unique opportunity for users to implement and test controllers on a space- 
based platform. However, the hardware environment has several features that must 
be accounted for when specifying control algorithms. First, the feedback control al- 
gorithms that can be implemented on ACTEX are fixed gain, and thus adaptive con- 
trollers cannot be used. Furthermore, these fixed-gain controllers are analog, which 
avoids sampling effects. Finally, the implementable analog controllers have a prespec- 
ified structure in which only the filter gains and natural frequencies can be modified. 
Since this constraint does not permit implementation of dynamic compensators of 
arbitrary structure, standard LQG and methods cannot be applied. 

We thus show that the decentralized static output feedback formulation of fixed- 
structure controller synthesis can directly account for the control-structure constraints 
of the ACTEX flight experiment. Specifically, we show that the ACTEX controller 
structure can be written as a decentralized static output feedback problem. Having 
done this, we then proceed to apply our techniques to obtain "^-optimal feedback 
controllers for suppressing broadband disturbances. 

The ACTEX flight experiment consists of a plate connected to a satellite by three 
struts, as shown in Figure 2.2. Each strut is equipped with its own control piezo- 
actuator as well as a colocated and nearly colocated sensor. A disturbance can be 
introduced to the experimental package through each of the three control actuators, 
or through a disturbance actuator on the plate. In addition, each of the three control 
actuators has an independent decentralized analog controller. 



Figure 2.2: ACTEX flight experiment 


Several experiments have been run on the ACTEX package, with telemetry re- 
turned for identification purposes. Measurements were taken from onboard accelerom- 
eters, sensor inputs, and actuator outputs, sampled at 4 kHz. Using this data, iden- 
tification was performed on the plant from the strut 1 actuator to the strut 1 nearly- 
colocated sensor. For frequencies below 200 Hz, these dynamics can be represented 
by (2.23)-(2.24) where 


-1.56 

77.98 

0 

0 


' -1.10 ' 


' -1.10 

0 ' 

-77.98 

-1.56 

0 

0 

, B = 

-0.21 

, D,= 

-0.21 

0 

0 

0 

-0.25 

24.99 

-0.11 

-0.11 

0 

0 

0 

-24.99 

-0.25 


0.03 


0.03 

0 



C — [ —0.15 0.70 1.83 6.88 ], D = 1, 


D 2 =[ 0 0.1 ] . 


Finally, we define the performance signals (2.26) based on the output measurements 
so that the closed-loop system reduces the vibrations seen by the sensors. Specifically, 
we have 


E ! = 


-0.15 0.70 1.83 6.88 
0 0 0 0 ’ 


E 2 = a 


0 

1 * 
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2.5. Fixed- Structure Synthesis for the ACTEX 
Flight Experiment 


Since continuous-time controllers are implemented on ACTEX, hardware con- 
straints place a limit on the form that the controllers may take. Each of the three 
struts on ACTEX has a decentralized controller, available in three fourth-order con- 
figurations and a sixth-order configuration, as shown below 


Gci(s) 

G c 2 (s) 

Gcz(s) 


k 3 u\ 

k 4 u>% 

s 2 + 0.3u;is + Ui 

+ s 2 + 0.3u 2 s + cu 2 


WlS 

s 2 T 0.3w 2 s + w 2 

S 2 + 0.3CU! s + cu 2 ’ 

k 3 wl 


S 2 + 0.3cd2^ + 

+ 0.3u x s + uJi 


(2.36) 

(2.37) 

(2.38) 


/ V ^3^3 . ^2 . 

c4 S s 2 + 0.3o;3S + s 2 + 0.3W2S + w 2 s 2 + 0.3u;is + ai 2 ’ 
where wi, w 2 , u 3 , fc 3 , and k 4 are subject to the constraints 


(2.39) 


0 < u}\,u)2,Uz < 1024, 

-16 < h < 16, (2.40) 

-8<k 4 < 8. 

Now, we consider the controller associated with strut 1, and express the controller 
in the decentralized static output feedback framework. Controller configuration 1, 
given by (2.36), can be expressed in dynamic compensator form as 



0 

1 

0 

0 


' 0 ' 

A c = 

-«? 

0 

— 0.3cui 
0 

0 

0 

0 

1 

, B c = 

1 

0 


0 

0 


—0.3w 2 


1 


C c = [ k z u\ 0 k 4 u\ 0 ] . 


Controller configurations 2 through 4, given by (2.37)-(2.39), can be expressed simi- 
larly and are given in the Appendix. 



Having specified the form of the controller, we consider the K 2 synthesis problem. 
The W ,2 norm of the closed-loop transfer function G 2tl ,(s) is given by (2.30). In order 
to design T^-optimal controllers for the ACTEX experiment, we pose the following 
optimization problem: Determine K that minimizes 

J(IC)=tvQR, (2.41) 


where QgN" satisfies (2.31). The necessary conditions for optimality can be derived 
by forming the Lagrangian (2.32). The partial derivatives with respect to Q and 
P in (2.32) are given by (2.33) and (2.34). To obtain the partial derivative of the 
Lagrangian with respect to the free parameters in the controller gains, we first specify 
the controller configuration. As an example, we consider controller configuration 1, 
given by (2.36). The settings for each of the other controller configurations are given 
in the Appendix. 

For controller configuration 1, the block-diagonal matrix JC has the form (2.28). 
Note that this controller has four free optimization parameters; namely, Wj, o> 2 , k 3 , 
and £4. Thus, we construct the matrix 


u\ 0 0 0 

0 0J2 0 0 

0 0 k 3 0 ’ 

0 0 0 k 4 


and note that 


K = K 0 + L X KR , + L 2 KKR 2 + L 3 KKM 3 KR 3 , 
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where 


K 0 = 
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0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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0 

0 

0 
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" 0 

0 
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'0 

0 

0 

o' 
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0 
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-1 

0 

0 
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0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 
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0 

0 

0 

0 

0 

- 0.3 

0 

0 


0 

-1 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

, L 2 = 
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0 

0 

0 
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0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 
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0 

0 

0 

0 

0 

0 

0 
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0 

0 

0 
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0 

0 
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0 

1 

0 

0 

0 

0 
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0 

0 

0 


0 

0 
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M 3 = 


0 0 
0 0 


0 0 0 0 
0 0 0 0 


Now with K. in terms of K , we can take the derivative of the Lagrangian with respect 
to K. Define the following notation 

r = I + VlL?K\ A k rBlPQCjL^, 

D k tBZPDV^L?, E k rVlEQC^L 


The derivative of the Lagrangian with respect to K is given by 

iff = diag {lJArJ + l]Ar*k + klZarZ + lJArJkmJk + klJArJkmJ 
+MJKKLJArJ + LjDRj + LlDR^K + KL^DRZ + L jDRjKMjK 



+KLjDRjKMj + MjKKLjDRj + LjERj + LjER^K + KLjER% 

+lJerJkmJk + klJerJkmJ + mJkklJerJ}. 


2.6. Simulation Results 

As mentioned in Section 2.3, one requirement of gradient-based optimization al- 
gorithms is an initial stabilizing design. Initial designs showed that for large values of 
u>i, i — 1,2,3, the cost function depended very weakly upon these values. Therefore, 
Ui, i = 1,2,3, were initially chosen to be 48, 72, and 96, respectively. Furthermore, 
since the open-loop system was stable, an initial stabilizing design could be obtained 
by setting k 3 = k 4 = 0. With these values set, the BFGS quasi-Newton algorithm 
was applied to find the T^-optimal solution for a given fixed-structure controller 
configuration. 

The ACTEX experiment is a lightly damped flexible structure, as can be seen by 
the impulse response shown in Figure 2.3, where x(0) = [ 0 0 0 1 ] T . For each of 
the four different controller configurations, three different controllers were designed 
by setting the control weighting matrix R 3 to 1, 0.01, and 0.0001. The cost of 
these controllers can be seen in Table 2.1. It is seen that the R 2 performance is best 
for the first controller configuration. In fact, the fourth controller configuration ran 
into parameter limits (2.40), as explained in Section 2.3. Therefore, the optimization 
routine terminated due to a boundary constraint rather than a small gradient condi- 
tion. It can also be seen from Table 2.1 that the cost of the closed-loop system 
decreases with increasing controller authority. 

The optimal controller parameters for the first controller configuration are given 
in Table 2.2. We can see from the output signals of the impulse response in Figure 2.4 
that with i ?2 = 1 the controller does not attenuate the vibrations significantly. How- 
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r 2 = 1 

R 2 = 0.01 
R 2 = 0.0001 


Config. 1 
0.5825 
0.2444 
0.1910 


Config. 2 
0.6774 
0.5601 
0.5495 


Config. 3 
0.6774 
0.5601 
0.5495 


Config. 4 
0.6782 
0.5613 
0.5507 


Table 2.1: R 2 costs for various controller configurations and weightings 



Wi 




fl 2 = i 

26.6460 

80.0602 

0.1352 


R 2 = 0.01 

36.4434 

99.1179 

0.3195 

law 

R 2 = 0.0001 

37.3664 

123.8404 

0.2081 

0.6314 


Table 2.2: Optimal controller parameters for controller configuration 1 
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ever, as the controller authority is increased, the attenuation becomes greater, as seen 
in Figure 2.5 and Figure 2.6, though this does increase the control effort expended, 
as shown in Figures 2. 7-2. 9, which could lead to actuator saturation. Finally, Fig- 
ure 2.10 shows the Bode plots of the open-loop and closed-loop systems. It is seen 
that the size of the first peak in the closed-loop response is decreased as the con- 
troller authority is increased. Another important feature of the closed-loop frequency 
response is the high-frequency roll-off. Since the ACTEX system model has a relative 
degree of zero, designing a controller without high-frequency roll-off would result in 
a closed-loop system which possesses gain at all frequencies, thus unmodeled high- 
frequency dynamics could destabilize the closed-loop system, whereas these dynamics 
would be attenuated with a strictly proper dynamic controller. 
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2.7. Conclusion 


In this chapter we introduced the decentralized static output feedback format. 
Specifically, we showed how dynamic output feedback control problems can be trans- 
formed into a decentralized static output feedback form. By using this format, a 
numerical optimization scheme can be used to optimize the controller gains with re- 
spect to a given cost function and constraint equation. Furthermore, we demonstrated 
the effectiveness of this framework on the ACTEX flight experiment. 


CHAPTER 3 

Stable ^"Optimal Controller 
Synthesis 

3.1. Introduction 

It is well known that even though LQG synthesis will stabilize the nominal closed- 
loop system, it can produce controllers with unstable dynamics, especially at high 
authority levels. Of course, for certain plants, specifically those that do not satisfy 
the parity interlacing property [121], only unstable controllers are stabilizing. How- 
ever, even for stable plants, LQG synthesis often produces unstable controllers, thus 
requiring Nyquist encirclements of the critical point. These encirclements and the 
resulting multiple gain margins, which must be maintained by the input actuators, 
can be jeopardized by actuator saturation and startup dynamics [87]. Therefore, 
whenever possible, it is desirable to implement only stable controllers. 

Several modifications of LQG theory have been proposed to obtain stable com- 
pensators. Several of these techniques involve either modified Riccati equations 
[64,112,113] or constrained weights [59,60]. Thus the resulting controllers may sacri- 
fice performance for controller stability. In [43,86,87], an augmented cost technique 
was proposed to obtain stable controllers without unnecessarily sacrificing perfor- 
mance. However, even though the authors in [87] give an excellent discussion on the 


implementation issues surrounding stable versus unstable controllers, they focus on 
multiple-model control, and therefore the cost function to be optimized is a weighted 
average of a number of system costs and does not give any insight into the trade-off 
between system performance and controller stability margin. 

The purpose of this paper is to provide a control-system design framework for % 2 - 
optimal strong stabilization. To achieve this goal we formulate the ^-optimal stable 
control problem within the context of decentralized static output feedback control 
which provides a general framework for fixed-structure dynamic controller synthesis 
[14,34], In particular, in order to guarantee stable stabilization, a multiobjective 
problem, reminiscent of scalarization techniques for Pareto optimization, is treated 
by forming a convex combination of the ^ norm of the closed-loop system and a 
weighted % 2 norm of the controller. It is shown that as the trade-off parameter is 
varied to obtain better H 2 system performance, the controller eigenvalues approach 
the imaginary axis. Thus the control engineer can decide if additional performance 
improvements warrant the resulting reduction in the stability margin of the controller. 

Two examples from the stable stabilization literature are considered here. The 
first example is a second-order spring-mass-damper system and the second example 
is a fourth-order two-mass system involving two flexible modes. The H 2 cost of the 
stable controllers developed for the first example, though larger than that of the LQG 
controller, was comparable to the lowest cost possible by a stable controller. For the 
second example, the difference between the cost of the stable controller and the 
unstable LQG controller is negligible. 
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3.2. Stable H 2 - Optimal Control 


In this section we state the "^-optimal stable stabilization problem. Specifically, 
given the n th -order plant 



x(t) = Ax(t ) + Bu(t) + Diw(t), t e [ 0 , 00), 

(3.1) 


with noisy measurements 



y(t) = Cx(t) + D 2 w(t), 

(3.2) 


and performance variables 



z(t) — Eix(t) + E 2 u(t), 

(3.3) 


determine an n* h -order strictly proper dynamic compensator 



— A c x c (t) T 

(3.4) 


w(t) — 

(3.5) 

— 

such that the % 2 performance criterion 



1 

J(A C , B c , C c ) = lim -E / 2 T (s) 2 (s) ds, 
t Jq 

(3.6) 


is minimized and the compensator dynamics matrix A c is asymptotically stable. 


The closed-loop system (3.1)-(3.5) is given by 



x(t) = Ax(t) + Dw(t), t G [ 0 , 00), 

(3.7) 


z(t) = Ex(t). 

(3.8) 


The closed-loop transfer function from disturbances w to performance variables z is 
given by 

G„{s) = E(sI a - A)~ x b. 


v 
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Next, we define a weighted controller transfer function from plant output y to plant 
input u by 

G c (s) = E 2 C c (sI nc - A c )~ 1 B c D 2 . 


Hence, the fixed-structure "^-optimal stable control problem is defined as: 

min 

(4c,Bc,C c ) 


subject to 


WGcWWl < °°- 


3.3. Design Equations 


The n 2 norm of G zw (s) is given by 

\\G zw (s)\\ 2 2 = tvQR, (3.9) 

where Q G N" is the unique nonnegative definite solution to the algebraic Lyapunov 
equation 

0 = AQ + QA t + V. (3.10) 

Furthermore, if A c is stable, then the % 2 norm of the weighted transfer function of 
the controller G c (s) = E 2 C c (sI nc - A c )~ l B c D 2 is given by 

\\G c (s)\\l = tvQ c CjR 2 C c , (3.11) 


where Q c € N" c is the unique nonnegative definite solution to the algebraic Lyapunov 
equation 

0 = A C Q C + Q c Aj + B c V 2 BJ. (3.12) 


To obtain (3.11) and (3.12) in terms of /C, given by 



A c 0 0 

0 B c 0 , 

0 0 C c 


(3.13) 


34 



we define the matrices Q Uj and Q KiJ , i = 1,2,3, j = 1, as given by (2.16), so that 

A c = Ql.ICQl,, 

Be = Ql n ICQl n , 

Cc = QI.ICQ1 

Thus, (3.11) and (3.12) become 

IIGcWll! = tr QcQrv&Ql^QLICQI,,, (3.14) 

and 

0 = Ql u ^-Qr u Qc + QcQku^Qu. + Qh-n ^-Qr 2 i ^2<5r 2 i ^ Qh 2X , (3.15) 

respectively. 

In order to design T^-optimal stable controllers we pose the following multiobjec- 
tive optimization problem: For p G [0, 1], determine /C that minimizes 

J{K) = (1 -p) tr QR + ptv Q c Qr 31 JC t Ql 31 R2Q1 31 K:q1 31 , (3.16) 

where Q e N”, Q c 6 N" c satisfy (3.10) and (3.15), respectively. Note that (3.16) 
involves a convex combination of the H 2 norm of the closed-loop system and the 
weighted W 2 norm of the controller. By including the %2 cost of the controller in 
the objective function, we can guarantee that the controller is stable as long as the 
objective function is finite. By varying p € [0, 1], (3.16) can be viewed as the scalar 
representation of a multiobjective cost. To achieve the best closed-loop performance 
with a stable controller, we only want to use the 'H 2 cost of the controller as a con- 
straint, and thus we set p > 0 to be small so that the contribution to the multiobjective 
cost due to the H 2 cost of the controller is negligible compared to the contribution 
due to the 'H .2 cost of the closed-loop system. Thus, the optimization routine will 


minimize the cost of the closed-loop system and not attempt to minimize the H 2 
cost of the compensator. However, increasing the trade-off parameter p will increase 
the controller stability margin. Finally, note that by letting p — > 0, we recover the 
^-optimal control problem. 

The necessary conditions for optimality can be derived by forming the Lagrangian 


£(P, Q, P c , Q c , K) = (1 - p)trQR + ptvQ c Q R31 IC T Q L31 R 2 Ql 3l JCQl 3l 


+tr P [AQ + QA r + V] + tr P c [Ql u fCQl u Qc 


+QcQRn£ T QLn+QL£Ql 21 V2QR 2l JC T Q Ui ), (3.17) 


where Peff and P c € N" c are Lagrange multipliers. The partial derivatives with 
respect to the free parameters in (3.17) are given by 


dC 

dA c 


f£ = A t P+PA + (1-p)R, 
oQ 


d£ 

dP 

dC 

dQ c 

dC 

dP c 

- qT 

^Ln qjq Mu ? 


= AQ + QA T + V, 

= AjP c + P c A c + pCjR 2 C c , 


= A C Q C + Q c Aj + , 


dC 

dB ' 


= Q 


dC 


Ql ... 


C\S~1 ^vLsi R31 ’ 


Ls, a/C Rsi ’ dCr h31 dlC 


where 


dC 


^ = B^PQCJ + B^PVVy W + (1 - p)Vj u EQCj 


+pQ h3l R 2 Ql 31 ICQl 3l Q c QR il +QLnPcQcQRn+QL 2l Ql 21 ICQl 2i V2QR n . 


3.4. Optimization Algorithm 


As noted in Section 2.3, one requirement of gradient-based optimization algorithms 
is an initial stabilizing design. For plants satisfying the parity interlacing property, 
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initialization can be accomplished by using sufficiently low authority compensators 
[24]. This was accomplished here by multiplying the control weight E 2 by a scalar 
77 > 1 . At sufficiently low authority, the LQG controllers tended to be stable. These 
low authority, stable, full-order controllers generally can then be truncated using 
an appropriate model reduction technique without destroying closed-loop stability. 
For decentralized control, this technique can be implemented in a sequential manner 
for each channel to obtain initializing gains with the given structure. These low 
authority LQG designs are used to initialize a low authority optimization algorithm. 
The optimized controller gains are then used to sequentially initialize higher authority 
problems until eventually the desired controller authority level is regained. At this 
point, the trade-off parameter p is varied until the best H 2 performance is attained 
in the face of a desired controller stability margin. 


3.5. Spring-Mass-Damper Example 


Consider the spring-mass-damper system given by the state space realization [42, 

60] 


x{t) = 

V(t) = 


0 1 
-3 -4 

„ *J 

[ 21 ] x(t). 


x(t ) + 


0 

1 


u(t), 


The matrices D x , D 2 , Ei, and E 2 are chosen to be 


D x 


35 0 
-61 0 ’ 


D 2 = [0 1 ], Ey = 


52.9150 8.9443 
0 0 


, £2 = 


0 

1 


For the given data, the LQG controller is unstable. To initialize the stable 

control problem, the control weighting was increased by multiplying E 2 by 77 = 16, 
as described in Section 3.4. This stable LQG design was used as a starting point 
for the quasi-Newton algorithm which found optimal stable compensators as 77 was 


decremented back to unity, returning the control authority to its original value. 



Figure 3.1: Bode plots of LQG and stable 7i 2 controllers ( p = 0.0288) 


As can be seen in Figure 3.1, increasing the authority toward a critical level 
(77 = 8.58) causes the gain of the LQG controller to approach infinity, at which point 
the low frequency phase jumps —180° and the gain begins to decrease, though the 
controller must now be unstable to maintain closed-loop stability. At this point, it 
can be seen that the gain of the ^-optimal stable controllers increase as well, though 
not as drastically as the LQG design, and the T^-optimal stable controllers always 
have a phase of 0 ° at low frequencies. The Nyquist plots of the LQG controller and 
the stable controller at full control authority, as seen in Figure 3.2, show the poor 
gain margins of the unstable LQG controller. 

Since the loop gain with the "Hs-optimal stable controller in feedback is much 
larger than that with the LQG controller at full authority, the impulse responses of 
the LQG controller and the stable controller were simulated to compare the actual 
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Figure 3.2: Nyquist plots of the loop gain for LQG (left) and stable “Hz (right) 
controllers (77 = 1 , p = 0.0288) 

control effort needed to bring the closed-loop system back to the equilibrium. These 
comparisons are shown in Figure 3.3. As expected, the performance of the system 
with the LQG controller is better than the performance of the system with the Tiz- 
optimal stable controller. However, note that even though the loop gain of the stable 
controller is much larger than that of the LQG controller, the control effort used by 
the stable controller is significantly less than that used by the LQG controller and 
hence is less likely to saturate the system actuators, which could cause closed-loop 
instabilities to occur when an unstable controller is used. 

Once the control authority was increased to the desired level, the value of the 
parameter p was varied to study the trade-off between the Hz cost of the system and 
the stability level of the compensator. Figure 3.4 shows the position of the controller 
(n c = 2 ) eigenvalues as a function of p as well as the Hz cost of the closed-loop system 





Time (sec) 



Figure 3.3: Impulse response of closed-loop system with LQG and stable I-L 2 con- 
trollers ( 7 / = 1, p = 2.88 x 10 -8 ) 


40 





! V> 


V 




v-' 



1/p 



Figure 3.4: Location of controller eigenvalues and 7<2 cost versus p 


as a function of p. By studying the resulting trade-off, the control engineer can decide 
if subsequent cost reductions justify bringing the controller eigenvalues closer to the 
stability boundary. 

For p = 0.5, the controller transfer function is given by 


-274.21s - 4183.8 
Gc ^ “ s 2 + 38.324s + 43.539’ 


(3.18) 


which has eigenvalues at Ai = —37.152 and A 2 = —1.1719, while, for p = 5 x 10 7 , 
the controller transfer function is given by 

-504.93s - 9744.5 


G c (s) = 


s 2 + 60.543s + 0.086523’ 


(3.19) 


which has eigenvalues at Aj = —60.541 and A 2 = —0.0014292. 

Since, as stated in [42], the lowest possible cost via stable stabilization for this 
example is given by a 4 th -order controller, we used the present framework to obtain 
stable 4 th -, 6 th -, and 8 th -order controllers to quantify the benefits of expanded-order 
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control. Specifically, these expanded-order K 2 -optimal stable controllers were initial- 
ized by adding one, two, and three stable modes to a full-order stable LQG controller, 
at which point the optimization algorithm was applied. The corresponding closed- 
loop costs, computed at various levels of control authority, are shown in Figure 3.5. 
At lower levels of controller authority, the stable controllers have nearly identical costs 
to the LQG controllers. As the authority is increased, the closed-loop cost associ- 
ated with the stable controllers becomes noticeably worse than the LQG controllers, 
however it is noticeably better than the best stable LQG design (i.e., the LQG de- 
sign with rj = 8.58). At the specified control authority (77 = 1), the dependence of 
the augmented cost function on the controller cost was decreased by decreasing the 
variable p. 



Figure 3.5: "H 2 cost versus control weighting for various-order stable % 2 controllers 

(p = 0.0288) 

As seen in Figure 3.6, the full-order controller has the highest H .2 cost, followed 
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by the controller with one extra mode. The other two controllers, however, have 
nearly identical closed-loop costs, suggesting that arbitrarily high-order controllers 
may not achieve significant performance improvements. Since decreasing p decreases 
the dependence of the cost function J on the weighted H 2 cost of the controller, the 
controller loop gain increases greatly as p becomes smaller, as shown in Figure 3.7. 
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Figure 3.6: Stable H 2 controller cost versus p 



2 nd-Order 

4th-Order 

6 th-Order 

8 th-Order 

Fixed - Structure 

627.11 

622.20 

621.73 

621.67 

Ganesh [42] (Optimal) 

N/A 

622.73 

N/A 

N/A 

Ganesh (Sub-Optimal) 

N/A 

628.40 

N/A 

N/A 

Halevi [60] (First Result) 

678.97 

N/A 

N/A 

N/A 

Halevi (After Tuning) 

637.18 

N/A 

N/A 

N/A 

Stable LQG (77 = 8.58) 

713.02 

N/A 

N/A 

N/A 


Table 3.1: %2 costs for various stable stabilization techniques 


As can be seen in Table 3.1, the T^-optimal stable controllers obtained here com- 





pare favorably to earlier results. In [60], the design weights were constrained in such a 
way as to yield stable controllers. Even with a tuning procedure, however, those costs 
are larger than what was obtained using the present framework. In [42], a nonlinear 
programming approach was used to obtain 7{2*°pt* ma l controllers, but only for SISO 
systems. However, as can be seen, the % 2 cost is slightly larger than the fourth-order 
controllers designed here. In fact, the optimal controllers obtained in [42] had two 
poles located at the origin, and thus the controller was merely conditionally stable. 
The stability boundary was then pushed back to s = —0.5 to yield a stable controller. 
It should be noted that the cost obtained by the second-order controller synthesized 
using our method was 627.1, whereas in [60], it is stated that the minimum cost 
possible by a second-order stable compensator is 628, though this may be simply a 
numerical artifact. Also listed is the U 2 cost of an LQG design with t] chosen as low 
as possible while still yielding a stable controller, which is significantly larger than all 
other methods considered. 


3.6. Two- Mass Example 

Consider the dynamic system [113] shown in Figure 3.8. The equations of motion 

. *1 =y . *2 


U 



k 


rti\ 


m 2 


w 


Figure 3.8: Two-mass system 



for this system are given by 


rriixi(t) + k (xi(t) - x 2 {t)) = u(t), 
m 2 x 2 (t) + k {x 2 (t) - X\ (t)) = 0. 


Here we consider the case of a colocated sensor and actuator pair, where the output 
is given by y = Xy. Letting my = m 2 = k = 1 yields the plant state space realization 


x(t) 


— 1 

o 

o 

o 

1 


' 0 ' 

0 0 0 1 

x(t) + 

0 

-1 10 0 

1 

1-10 0 


1 

o 
1 

1 0 0 0 ] x(t). 



u{t), 


As in [113], the matrices Dy, D 2 , Ey, and E 2 are chosen to be 
0 0 

D = 00 

Dl 0 0 

_ 68 0 

For this example, a full-order and a reduced-order "H 2 -optimal stable controller 
were designed. The control authority was chosen to be sufficiently low so that the 
LQG controller was stable. This controller was then used to initialize the optimization 
algorithm. 

The same general trends can be observed here as in the first example, though in 
this case, even when full authority is achieved, the % 2 cost of the full-order stable 
controller rivals the performance of the unstable "^-optimal LQG controller, which 
would be apparent in Figure 3.9 if the curves were not directly on each other. In fact, 
the relative difference in the Ti 2 costs is merely 9.7114 x 10 -4 %. For this reason, 
expanded-order controllers could not give further cost improvements to justify their 
increased complexity, and thus were not designed for this example. Also, the initial 
value of p was sufficiently small that further reductions did not improve the H 2 cost of 


D 2 = [ 0 1 ] , E y = 


10 10 
0 0 0 0 


E 2 


0 

0.01 
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Figure 3.9: 7/ 2 cost versus control weighting ( p = 0.04752) 


the closed-loop system. Table 3.2 compares the performance of the 7/ 2 -optimal stable 
controller designs with the LQG controller, the best stable LQG controller (i.e., the 
LQG controller with 7 ? chosen as low as possible while still yielding a stable controller), 
and the full-order controller designed in [113], which used an over-bounding approach 
along with parameter tunings to obtain stable controllers. 





LQG 

N/A 

16.175703 

Fixed - Structure 

27.9340 

16.175782 

Wang [113] 

N/A 

16.261138 

Stable LQG (77 = 12.21) 

N/A 

34.334324 


Table 3.2: H 2 costs for various stable stabilization techniques 


The ability of this stable controller to achieve an 7^2 performance nearly identical 
to that of the LQG controller was then explored. After running numerous examples, 





it appears that a minimum phase open-loop stable plant (such as the first example) 
will yield a significant performance degradation when the controller is constrained to 
be stable. However, a minimum phase, open-loop unstable plant does not seem to 
exhibit this lack of performance, as demonstrated by this example. Further investi- 
gation seems to show that T^-optimal stable controllers designed for non-minimum 
phase, open-loop unstable plants will also show performance degradation over an LQG 
controller, whereas no appreciable loss of performance occurred when a %2-optimal 
stable controller was designed for a non-minimum phase, open-loop stable plant. 

3.7. Conclusion 

In this chapter we investigated a scheme to synthesize "^-optimal controllers by 
including the %2 cost of the controller in the Lagrangian and using a multiobjective 
optimization technique. The problem was formulated in a decentralized static out- 
put feedback framework, which facilitated the use of a quasi-Newton optimization 
algorithm. This technique was applied to two numerical examples. It was numeri- 
cally shown that for some systems, namely minimum phase, open-loop unstable or 
non-minimum phase, open-loop stable plants, a stable controller can rival the perfor- 
mance of an unstable T^-optimal LQG controller and yet not be constrained by the 
loop margins of unstable controllers. For other systems, however, there could be a 
significant degradation in performance by requiring the controller to remain stable, 
although this technique provided controllers yielding the minimal W2 closed-loop cost 
for all stable linear controllers. 
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CHAPTER 4 


"^"Optimal Synthesis of 
Controllers with Relative 
Degree Two 

4.1. Introduction 

It is well known that modern multivariable control design frameworks such as H? 
and H 0 0 control yield dynamic compensators with relative degree zero or one. Hence, 
the structure of the dynamic feedback controller is such that the measured system 
output appears explicitly in the control signal or in the control rate signal [85]. In 
the single-input/single-output system case, the resulting controller transfer function 
is non-strictly proper with relative degree zero or strictly proper with relative degree 
one. In this case, the Bode plot of the controller transfer function at best rolls off at 
20 dB per decade. Alternatively, for relative degree r controllers, the Bode plot of 
the compensator has a high frequency roll-off of 20r dB per decade. 

High frequency roll-off is particularly useful when the system under consideration 
is a lightly damped flexible structure. Since flexible structure models are by necessity 
truncated to a finite number of modes, it is desirable for the frequency response to 
roll off as quickly as possible after the gain crossover frequency so that unmodeled 
high frequency system dynamics are not excited by the controller dynamics. 


For single-input/single-output systems, where the H 2 norm corresponds to the 
area under the Bode plot and the Hoo norm corresponds to the maximum magnitude 
of the Bode plot, roll-off rates cannot be specified solely by minimization techniques on 
these norms. Loop shaping weighting functions can be used in the controller design 
process, but these specify the frequency where the roll-off starts, not the roll-off 
rate. Furthermore, these techniques also tend to result in high-order controllers when 
frequency weighting is included in the design process. In this section we extend the 
fixed-structure controller design framework of [14,95] and [34] to design "^-optimal 
relative degree two controllers for multi-input/multi-output systems. Since we cast 
the relative degree two design problem within the fixed-structure control framework, 
fixed-order (i.e. , full- and reduced-order) controllers can be designed with increased 
roll-off rates at the gain crossover frequency. Even though the proposed framework 
can be easily extended to include desired weighting functions for loop shaping, we do 
not do so here to facilitate the presentation. 

The proposed 'H 2 -optimal relative degree two controller design technique is applied 
to several structural control problems, showing that the resulting relative degree two 
controller incurs minimal increase in % 2 performance over the optimal LQG controller 
while enforcing a 20 dB per decade increase in the roll-off rate at the gain crossover 
frequency. 

4.2. % 2 -Optimal Relative Degree Two Control 

In this section we state the "^-optimal relative degree two control problem. Specif- 
ically, given the n th -order stabilizable and detectable system 

x(t) = Ax(t) + Bu(t) + Diw(t), t € [0, oo), 
y{t) - Cx(t) + D 2 w(t), 
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(4.1) 

(4.2) 



z(t ) = Eix(t) + E 2 u(t), 


(4.3) 


determine an n‘ h -order relative degree two dynamic compensator 


id (t) = A Cl x c i(t) + B Cl y{t), (4.4) 

±c 2 (t) = A C2 x c2 {t) + B C2 v(t), (4.5) 

v(t) = C Cl x c i(t), (4.6) 

u{t) = C C2 x c2 (t), (4.7) 


where x cl (t) G M ncl , x c2 (t ) G R Uc2 , and n c = n c i + n c2 , such that the U 2 performance 
criterion 


J(A Cl ,A C2 ,B Cl ,B C2 ,C Cl ,C C2 ) = lim jE [ [x T (s)fi 1 z(s) + u T {s)R 2 u{s)]ds, (4.8) 

t—>OC l J Q 


is minimized. 

Note that the dynamic controller (4.4)-(4.7) corresponds to a cascade interconnec- 
tion of the two controllers in the feedback loop (see Figure 4.1) so that the controller 
transfer function realization is given by 


G c (s) = G C2 G Cl (s)~ 


A c Be 
C c 0 


(4.9) 


where 


A r = 


A r . 0 

A 


*Cl 

Br.Q 


J C2^Cl 


C2 


B c = 


-'Cl 


, c c = [ 0 c H ] 


(4.10) 


Note that since B C2 is always multiplied with C Cl , B C2 C Cl can be considered a single 
free parameter, thus leaving only five controller gains over which to optimize, instead 
of six. Finally, we note that this framework can be easily extended to address the 
design of (vector) relative degree r controllers by considering a cascade interconnection 
of r dynamic controllers in the feedback loop. 



Figure 4.1: Relative degree two controller set-up 

4.3. Relative Degree Two Controller Synthesis 

In this section we use the decentralized static output feedback framework as de- 
scribed in Section 2.1 to design "^-optimal relative degree two dynamic compensators. 
Specifically, we define 

/C — block-diag [.4 Cl > A C3 , B Cl , B C2 C Cl , C C2 \ . 

Now, if A is asymptotically stable for a given feedback gain 1C € ]R( 2n c+m)x( 2 n c +i) ^ 
follows that the H 2 performance criterion (4.8) is given by 

J{K) =||6 w (s)||l = trPV, (4.11) 

where P € N” is the unique nonnegative definite solution to the algebraic Lyapunov 
equation 

0 = A t P + PA + R. (4.12) 

Now, the necessary conditions for optimality can be derived by forming the La- 
grangian 

C(P, Q,JC) = tr {PV + Q [i T P + PA + p] } , (4.13) 
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where Q € N™ is a Lagrange multiplier. The gradient expressions with respect to the 
free parameters in (4.13) are given by 


^ = AQ + QA T 4- V, 
dP 

= A T P + PA + R, 
dQ 

m = 


= 2QljT [BlPQCl + + Vj„EQC „ T ] L^Ql,, 

j=l, 2 , 3, 4, 5, j = l, 


(4.14) 

(4.15) 

(4.16) 


where r = I + and 


/Ci = A Cl , K ,2 — 44 C2 , /C 3 — Ben ^4 — B C2 C Cl , K- 5 C C2 . 

4.4. Control Design Process 

The initial stabilizing design for the relative degree two controller was obtained 
by computing the balanced truncation of the full-order LQG controller to obtain a 
reduced-order LQG controller corresponding to the controller in the feedback loop 
with the lower order. Note that since two controllers are being synthesized, it is not 
necessary for this first truncated controller to stabilize the system. This first designed 
controller was augmented to the plant. If G Cl is designed first, the augmented plant 
realization is ’ 


(4.17) 


' A 

0 

B ' 

| 

A ' 

B Cl C 

4* 

B C1 D\ 

B Cl D2 

0 

c 

[_o_; 

0 

. E x 

0 

e 2 ; 

E 0 


whereas if G C2 is designed first, the augmented plant realization is 



A 

bc C2 

0 



0 

^C2 

B c 2 

0 

_l 

rsj 

C 

dc C2 

0 



. E x ' 

E 2 C 'c 2 

' 0 
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Note that in the first case, the D 2 term is identically zero, whereas in the second 
case, the E 2 term is zero. Thus these augmented matrices result in a singular control 
problem. This was overcome by replacing these terms with nonzero matrices struc- 
tured such that DiDj = 0 or EjE 2 = 0, as appropriate. Once an LQG controller 
was designed on this “artificial” system, the dynamics of the original system were 
tested, and if the closed-loop was asymptotically stable, these designed controllers 
were used as the initial controllers for the gradient search algorithm. For details of 
the algorithm, see [35]. 


4.5. Two-Mass Example 


Consider a two-mass-spring-damper system with a colocated sensor/ actuator pair 
and state space realization in real normal coordinates given by 


" -0.0002 

0.2208 

0 

0 


‘ -0.1439 ' 

-0.2208 

-0.0002 

0 

0 

x(t ) + 

0.2168 

0 

0 

-0.0103 

1.4320 

-0.0426 

0 

0 

-1.4320 

-0.0103 


1.1890 


y(t) = [ -0.0545 0.0819 -0.0352 0.8181 ] x(t). 


The weighting matrices D\, D 2 , Ei, and E 2 were chosen so that LQG synthesis would 
place a notch at the second mode. This is accomplished when [40] 

' 0 0 

D - 1 0 

Dl ~ 0 0 

. 0 0 

For this system, the two controllers in the feedback loop (G Cl and G Ci ) were chosen 
to be of order two. Initializing reduced-order LQG controllers were designed, and the 
gradient search algorithm was initiated. The Bode plots of the loop gain of the full- 
order LQG controller and of the relative degree two controller, along with dotted lines 
representing the respective high frequency asymptotes, are shown in Figure 4.2. The 


D 2 = [ 0 1 ] , Ei = 


10 0 0 
0 0 0 0 


E 2 = 


0 

1 


54 




V 


v 


V 






v 






Figure 4.2: Bode plots of LQG and relative degree two controllers 


"^-optimal LQG cost is 3.9734 and the % 2 -optimal relative degree two controller is 

4.0743 which corresponds to only a 2.5% increase over the % 2 -optimal LQG controller 

cost. This marginal increase in the % 2 cost is not surprising since I-L 2 - optimal relative 

degree two controllers are sought. 

The transfer function for the LQG controller is 

-0.9090s 3 - 0.1073s 2 - 1.8660s - 0.1817 
G c( lqg)(s) “ s 4 + 0.3629s 3 + 1.2425s 2 + 0.6170s + 0.1443’ 

which has natural frequencies at 0.345 rad/sec and 1.10 rad/sec. The relative degree 

two controller transfer function is given by 

-99.9827s 2 - 2.0960s - 205.2481 

Gc(Rei.Deg.2)(s) - fi4 + 64.6193s 3 + 46.1170s 2 + 34.9138s + 77.6907' 

This transfer function has a natural frequency at 1.02 rad/sec and break frequencies 
at 0.862 rad/sec and 0.015 rad/sec. Thus it is seen in this case that constraining the 
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controller to be of relative degree two does not push the controller poles out to such 
high frequencies that the extra 20 dB/decade roll-off is not useful. 


4.6. Three Mass Example 


Consider the three-mass system given by the state-space realization ( [52]) 


x(t) = 

y(t) = 


— , 
0 

0 

0 

1 

0 

0 
1 


1 

0 
1 

0 

0 

0 

0 

1 

0 


0 

0 

0 

0 

0 

0 

1 

x(t) + 

0 

-1 

1 

0 

0 

0 

0 

0 

1 

1 

to 

1 

0 

0 

0 


0 

0 

1 

-1 

0 

0 

0 _ 


1 

1 1 

0 

0 0 

0] 

X 

(t). 




The weighting matrices D x , Z? 2) E x , and E 2 are chosen to be 


D x = 


0 0 
0 0 
0 0 
1 0 ’ 
0 0 
0 0 


D 2 — 0 1 ] , E x — 


1 1 0 0 0 0 
0 0 0 0 0 0 



For this system, two different relative degree two controllers were designed. The first 


was designed with the first controller in the feedback loop of order two and the second 
of order four. The Bode plots of the loop gain of the full-order LQG controller and 
of this relative degree two controller are shown in Figure 4.3. Again, the increase in 
the H 2 cost is only 9% above the optimal value. 

The transfer function for the LQG controller is 

—2.4070s 5 — 0.1479s 4 - 10.1266s 3 — 1.3616s 2 — 8.33745 - 2.000 
G C (lqg)(s) ~ 5 6 + 3.1221s 5 + 8.8738s 4 + 15.8141s 3 + 21.5012s 2 + 17.8508s + 13.3444’ 

which has natural frequencies at 1.30 rad/sec, 1.55 rad/sec, and 1.82 rad/sec. The 

relative degree two controller transfer function is given by 

-332.2810s 4 - 11569s 3 - 679770s 2 - 1186900s - 293090 
G c(RD2)(«) - s 6 + 182.4514s 5 + 12699s 4 + 392630s 3 -I- 4634100s 2 + 3606700s + 2212100' 


56 


w 


V* 








Figure 4.3: Bode plots of LQG and relative degree two controllers 

This transfer function has natural frequencies at 0.715 rad/sec, 34.6 rad/sec, and 
60.1 rad/sec. In this case, the relative degree two controller does have much higher 
frequency poles, so that the extra 20 dB/decade roll-off does not occur until higher 
frequencies are encountered, as seen in the figure. 

The second relative degree two controller was designed with the first controller in 
the feedback loop of order four and the second of order two. The Bode plots of the 
loop gain of the full-order LQG controller and of this relative degree two controller are 
shown in Figure 4.4. For this configuration, however, the H 2 cost of the relative degree 
two controller is less than that of the first relative degree two controller designed. In 
particular, the increase in the H 2 cost for this controller is only 2.9% above the 
optimal value. 

The transfer function of this relative degree two controller is 

— 130.3438s 4 - 582.6117s 3 - 258.0451s 2 - 679.4487s - 176.8302 
G C (rd2)(s) “ s e + 56.0271s 5 + 811.7752s 4 + 1118.0s 3 + 2196.9s 2 + 1630.2s + 1186.4’ 


v 


57 





Figure 4.4: Bode plots of LQG and relative degree two controllers 

which has natural frequencies at 0.975 rad/sec and 1.30 rad/sec, and break frequencies 
at 0.0417 rad/sec and 0.0326 rad/sec. In this case, the relative degree two controller 
does not have higher frequency poles, so that the roll-off is more pronounced than in 
the first configuration. 


4.7. Coupled Rotating Disk Example 


Consider the coupled rotating disk problem given in [22] with state space realiza- 


tion in real normal coordinates given by 


±(t) 


-0.1610 1 0 0 0 0 0 

-6.0040 0 1 0 0 0 0 

-0.5822 0 0 1 0 0 0 

-9.9835 0 0 0 1 0 0 

-0.4073 0 0 0 0 1 0 

-3.9820 0 0 0 0 0 1 

0 0 0 0 0 0 0 

0 0 0 0 0 0 0 


0 ' 


0 ' 

0 


0 

0 


0.00640 

0 

x(t ) + 

0.00235 

0 

0.07130 

0 


1.00020 

1 


0.10450 

0 


0.99550 
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V 


y(t) = [1 0000000] x{t). 

The weighting matrices D 2 , E\, and E 2 are chosen to be 



O 

O 


0 0 ' 


0 0 


0 0 

V 

0.00640 0 


0 0 


0.00235 0 

II 

q 

r*H 

0 

II 

IN 

Q 

0 0 

Di = 

0.07130 0 

0.00055 0 


1.00020 0 


0.01100 0 


0.10450 0 


0.00132 0 


0.99550 0 


0.01800 0 


For this system, two different relative degree two controllers were designed. The first 
was designed with the first controller in the feedback loop of order two and the second 
of order six. The Bode plots of the loop gain of the full-order LQG controller and of 
this relative degree two controller are shown in Figure 4.5. Here, the increase in the 




Figure 4.5: Bode plots of LQG and relative degree two controllers 


cost is a scant 0.0238% above the optimal value. 
The transfer function for the LQG controller is 
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„ , x -0.1102s 7 - 0.0210s 6 - 0.6637s 5 - 0.0843s 4 - 1.1105s 3 - 0.0793s 2 - 0.4500s - 0.0152 

G C (lqG)( s ) = s 8 + i. 3183 s 7 + 6.8596s® + 7.3937s 5 + 14.1279s 4 + 11.4920s 3 + 9.2199s 2 + 4.5193s + 1.3859’ 


which has natural frequencies at 0.548 rad/sec, 0.820 rad/sec, 1.41 rad/sec, and 1.86 
rad/sec. The relative degree two controller transfer function is given by 

-28.0400s 6 - 8051.5s 5 - 1122.6s 4 - 20797s 3 - 1486.7s 2 - 9550.1s - 321.8893 
Gc(rd 2 )(«) - s 8 + 569.3318s 7 + 75307s 6 + 93453s 5 + 244920s 4 + 215410s 3 + 186610s 2 + 95718s + 29437' 

This transfer function has natural frequencies at 0.549 rad/sec, 0.819 rad/sec, and 
1.40 rad/sec, and break frequencies at 0.00485 rad/sec and 0.00276 rad/sec. 

The second relative degree two controller was designed with the first controller in 
the feedback loop of order six and the second of order two. The Bode plots of the 
loop gain of the full-order LQG controller and of this relative degree two controller are 
shown in Figure 4.6. For this configuration, however, the I-L 2 cost of the relative degree 




Figure 4.6: Bode plots of LQG and relative degree two controllers 
two controller is slightly more than that of the first relative degree two controller 
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designed, but the increase in the 'H .2 cost for this controller is only 0.864% above the 
optimal value. 

The transfer function of this relative degree two controller is 

—0.2916a 6 - 0.1042a 6 - 0.7773s 4 - 0.1187a 3 - 0.3535s 2 - 0.0118s - 2.9883 x 10~ u 

Gc(RD2)(s)= s 8 + 3.1635a 7 + 7.3610a 6 + 10.8948s 5 + 11.4236s 4 + 8.2630a 3 + 3.9038a 2 + 1.0973a + 2.0906 x 10“ 9 ’ 

which has natural frequencies at 0.548 rad/sec, 0.820 rad/sec, 1.41 rad/sec, and 1.86 
rad/sec. 

4.8. Conclusion 

In this chapter we proposed a scheme to synthesize T^-optimal relative degree two 
controllers by cascading two controllers in the feedback loop and optimizing over the 
five free controller parameters. The problem was formulated in a decentralized static 
output feedback framework, which facilitated the use of a quasi-Newton optimization 
algorithm. This technique was applied to three numerical examples. It was shown that 
constraining the controller to have a relative degree of at least two only marginally 
increased the cost of the closed-loop system, though it was noted that the order 
of the separate controllers in the feedback loop does significantly affect the H .2 cost 
of the closed-loop system and the natural frequencies and break frequencies of the 
controller dynamics. 




CHAPTER 5 


Robust Fixed-Structure Controller 
Synthesis using the Implicit Small 
Gain Bound 

5.1. Introduction 

One of the principal objectives of robust control theory is to synthesize feed- 
back controllers with a priori guarantees of robust stability and performance. In 
/i-synthesis [37, 97] these guarantees are achieved by means of bounds involving 
frequency-dependent scales and multipliers which account for the structure of the 
uncertainty as well as its real or complex nature. An alternative robustness approach 
involves bounding the effect of real or complex uncertain parameters on the K 2 per- 
formance of the closed-loop system. These guaranteed cost bounds take the form of 
modifications to the usual Lyapunov equation to provide bounds for robust stability 
and U 2 performance. 

A diverse collection of guaranteed cost bounds have been developed. An overview 
of many of the early guaranteed cost bounds can be found in [12], while positive- 
real-type guaranteed cost bounds are discussed in [47]. More recently, Popov-type 
guaranteed cost bounds have provided links with frequency-dependent scales and 
multipliers while providing reliable bounds for the peak real structured singular value 


[15,52,106]. Finally, the introduction of shift terms has been shown to reduce the 
conservatism of guaranteed cost bounds [38, 54, 74] for structured real uncertainty 
without requiring frequency-dependent scales and multipliers. 

The goal of this research is to explore the applicability of the implicit small gain 
guaranteed cost bound of [54] to controller synthesis. As shown in [54], unlike the 
quadratic stability bounded-real type bound of [76,100], the implicit small gain bound 
can distinguish between real and complex uncertainty and is particularly effective in 
capturing internal uncertainty structure. For flexibility in controller synthesis, we 
adopt the approach of fixed-structure controller synthesis [62] which allows consid- 
eration of arbitrary controller structures, including order, internal structure, and 
decentralization [34], 

To demonstrate the fixed-structure/implicit small gain approach to robust con- 
troller synthesis, we consider two examples that have been addressed by means of 
alternative guaranteed cost bounds. The first example, which involves two flexible 
modes, was used in [40] to illustrate the Maximum Entropy technique, while the 
second example, which involves three flexible modes, was considered in [52, 106] to 
illustrate fixed-structure Popov synthesis. 

5.2. Robust Stability and Performance Problem 

In this section we state the robust stability and performance problem. This prob- 
lem involves a set U C R nxn of uncertain perturbations da of the nominal system 
matrix A. The objective of this problem is to determine a fixed-order, strictly proper 
dynamic compensator (A c , B c ,C c ) that stabilizes the plant for all variations in U and 
minimizes the worst-case W 2 norm of the closed-loop system. In this section and the 
following section no explicit structure is assumed for the elements of U. In Section 5.4 
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the structure of U will be specified. 


Robust Stability and Performance Problem. Given the n th -order stabiliz- 


able and detectable system 

x(t) = (A + AA) x(t) + Bu(t) + Diiv(t), t € [0, oo), (5.1) 

y(t) = Cx(t) + D 2 w(t), (5.2) 

determine an n* h -order dynamic compensator 

x c {t) = A c x c (t) + B c y(t), (5.3) 

u(t ) = C c x c (t), (5.4) 


such that the closed-loop system (5. 1)— (5.4) is asymptotically stable for all AA € U 
and the performance criterion 

1 ft 

J(A c ,B c ,C c )= sup limsup-E / [x T (s)i?ia:(s) + u r (s)R 2 u(s)] ds, (5.5) 
A AeU t-> oo t J 0 

is minimized. 

For each uncertain variation AA € U, the closed-loop system (5. 1 )— (5.4) can be 
written as 

x{t) = (A + AA)x(t) + Dw(t), £ € [0, oo), (5.6) 

A,4 O' 

0 0 ‘ 

5.3. Sufficient Conditions for Robust Stability and 
Performance 

In practice, steady-state performance is only of interest when the undisturbed 
closed-loop system is robustly stable over U. The following result is immediate. 

Lemma 5.1. Let ( A c , B c , C c ) be given and assume that ^4 + is asymptotically 


where AA = 


stable for all AA E U. Then 


J{A c ,B c ,C c ) = sup tr P aA V, (5.7) 

AAeU 

where P A A E N” is the unique nonnegative definite solution to 

0=(i + Ai) T P a/i + P a/i (i + A4) + «. (5.8) 

The key step in guaranteeing robust stability and performance is to replace the 
uncertain terms in the Lyapunov equation (5.8) by a bounding function 51. 

Theorem 5.1 [12]. Let ( A c , B c ,C c ) be given, let 51 : N" — ► S n be such that 

A,4 T P + PAA < S1(P), AAeU, P E N*, (5.9) 

and suppose there exists Pelf satisfying 

0 - A T P + PA + 51(P) + R. (5.10) 

Then (.4 + A A, E) is detectable for all A.4 E U if and only if A + AA is asymptotically 
stable for all A.4 E U. In this case, 

P AA < P, AA E U , (5.11) 

where ■^A/4 is given by (5.8), and 

J(A C , B c , C c ) < tr PV. (5.12) 

5.4. Uncertainty Structure and the Implicit Small 
Gain Guaranteed Cost Bound 

We now assign explicit structure to the set U and bounding function 51 (P). Specif- 
ically, the uncertainty set U is defined by 

U - |a, 4 E R nxn : AA = Y^SiAi, |5»| < 7 -1 , * = !,•• • ,r j, (5.13) 
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where 7 is a positive number and, for i = 1, . . . , r, A{ G R nxn is a fixed matrix denoting 

the structure of the parametric uncertainty and Si is an uncertain real parameter. Note 

that U given by (5.13) includes repeated parameters without loss of generality. For 

example, if <5i = S 2 then replace A x with A x + A 2 and discard S 2 and A 2 . Furthermore, 

" S 5 ~ 

U includes real full block uncertainty. For example, if AA = ^ ^ , then Av4 = 

where A x = * and likewise for A 2 , ^3, and A 4 . Finally, for i = 

1 , . . . , r, letting A { = BiCi, where B { G R nx<, ‘, Q G R fl<xn , and q x < n, and defining 
B 0 = [ B x ••• B r ] and Co =[ Cj Cj ] T , U can be written as the real 
parameter uncertainty set considered in [37] 


U = { A.4 G R nxn : A ,4 = B 0 AC 0 , A G A 7 } , 


(5.14) 


where 

A 7 = (A G R ,x<? : A = block-diag[5i/ 9 , . . . ,S r I Qr ] , |<5*| < = 1, . - - , r*} , (5.15) 

and where q = X)i=i9«- Since an uncertainty set of the form (5.14) can always 
be written in the form of (5.13) by partitioning B 0 and C 0 as above and defining 
Ai = BiCi, i = 1 , . . . , r, robust stability of A 4- A A for all Aj4 G U is equivalent to 
the robust stability of the feedback interconnection of 6r(s) = Co(sIa — A)~ l Bo and 
A, where B 0 = [ Bj 0 ] T , C 0 = [ C 0 0 ] , and A G A 7 . 

With the uncertainty set U given by (5.13), the closed-loop system (5.6) has 
structured uncertainty of the form 


A A = Y i 6 i A i , 


(5.16) 


1=1 


where 




Ai 0 
0 0 


i = l,...,r. 
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We now introduce a specific choice of the bound f2(P) satisfying (5.9) for the 
structure of U as specified by (5.13). For i = 1, . . . , r, let Si G R nxn and define 

^ = [(S, + S t T ) 2 ]*. 

Note that -Zi < a(S t + Sj) < Zi for all a G [-1,1]- If Si is skew symmetric 
then Zi = 0. Furthermore, for i — define = [5j A?][.Fi where 

(•) t denotes the Moore-Penrose generalized inverse. Note that Ii is symmetric and 
idempotent, that is, /* = Ij = If. Furthermore, since IpSi Aj] = [Si 
follows that Ii§i = Si and AJi = A*. If Si = A* and Ai is an EP matrix [21], that is, 
7Z(Ai) = IZ(Aj), then = A\Ai. Recall that normal matrices (and thus symmetric 
and skew-symmetric matrices) are EP. 


Proposition 5.1 [54]. For i = 1, . . . , r, let a* G K, $ > 0, and let S) G R nxn . 
Then (5.9) is satisfied with Q(P) given by 

H(P) = it [7" 2 (« l 5 i +A^) T («.5i+Ai^)+7" 1 /3rV»l^+/5r 2 4] ■ (5-17) 

t=i 

Remark 5.1. As discussed in [54], if Si is assumed to be complex for some z, then 
it can be shown that f2(P) given by (5.17) does not satisfy (5.9). Hence, unlike the 
quadratic stability bound of [76,100], the bound (5.17) can distinguish between real 
and complex uncertainty. 

Combining Theorem 5.1 with Proposition 5.1 yields the following result. For 
convenience define the shifted dynamics matrix A S7 = A + 7 _2 X)[ = i otifiiAiSi. 


Theorem 5.2. For i = 1 , . . . ,r, let a, gR, A > 0, and let £ K nxn - Further- 
more, suppose there exists a matrix PGff satisfying 

0= A^ AbPA 87 +^ [7 -2 ( a <‘5< r Si+PfPAiAj P)+7 -1 ^ t -1 |ai|^i+/3“ 2 /ij+R. (5.18) 

i=l 
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Then (A+AA, E) is detectable for all A A G U if and only \iA + AA is asymptotically 
stable for all AA G U. In this case, 

P aA <P, A AeU, (5.19) 

where Paa satisfies (5.8), and 

J(A c ,B c ,C c )<trPV. (5.20) 

Remark 5.2. If T s7 is asymptotically stable then the existence of a matrix P G 
N* satisfying (5.18) is equivalent to the existence of a frequency-domain condition 
guaranteeing robust stability of A + AT, AA G W, in terms of an implicit small 
gain condition involving the shifted dynamics matrix T S7 which is a function of the 
uncertainty set bound 7 . For details see [54]. 

To apply Theorem 5.2 to controller synthesis, we use the Riccati equation (5.18) to 
guarantee that the closed-loop system is robustly stable. This leads to the following 
optimization problem. 

Optimization Problem. Determine ( A c , B c , C c ) that minimizes J (A c , B c , C c ) = 
tr PV, where P G N" satisfies (5.18) and such that (A c ,B c ,C c ) is controllable and 
observable. 

The relationship between the Optimization Problem and the Robust Stability and 
Performance Problem is straightforward, as shown by the following proposition. 

Proposition 5.2. Let (A c , B c ,C c ) be given. If P G N* satisfies (5.18) and (^4 + 
Aj4, E) is detectable for all A,4 G U, then A + A A is asymptotically stable for all 
A.4 G U, and J(A C , B c> C c ) < J{A C , B c , C c ). 

Proof. Since (5.18) has a solution P G If and (^4 -h AA,E) is detectable for 
all AT G U, the hypotheses of Theorem 5.2 are satisfied so that robust stability and 


robust performance are guaranteed. Now, J(A C , B c ,C c ) < J{A C , B c , C c ) is merely a 
restatement of (5.20). □ 

It follows from Proposition 5.2 that the satisfaction of (5.18) along with the generic 
detectability condition leads to robust stability along with an upper bound for the Hi 
performance. Hence, by deriving necessary conditions for the Optimization Problem 
we obtain sufficient conditions for characterizing dynamic output feedback controllers 
guaranteeing robust stability and performance. 

5.5. Robust Controller Synthesis via the Implicit 
Small Gain Guaranteed Cost Bound 

In this section we state constructive sufficient conditions For characterizing fixed- 
order (i.e., full- and reduced-order) robust controllers. These results are obtained 
by minimizing the worst-case 'Hi cost bound (5.20) subject to (5.18). The following 
result is required for the statement of the main theorem. 

Lemma 5.2 [62]. Let Q, P G N” and suppose that rank QP = n c . Then there 
exist n c x n matrices G, T and an invertible matrix M e R" cX ” c , unique except for a 
change of basis in M" c , such that 


QP = G t MT, rG T = J„ c . (5.21) 

Furthermore, the n x n matrices r = G T T and r± = I n — r are idempotent and have 
rank n c and n — n c , respectively. 


To apply Theorem 5.2 to robust controller synthesis, let Si, i = 1, . . . , r, have the 
form 

*5i 0 n xn c 


Si = 


0n c xn 0 nc xn c 


(5.22) 
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where S* G R nxn . With Si, i = 1, . . . ,r, given by (5.22) it can be shown that 




w 


v 


/t On x n c 

£ — 

^z 0 nX n c 

0 nc xn 0 nc xn c . 

> Z/j 

0 nc xn 0 Uc xn c 


where R = [S* AL][S,- A?]* and Z i = [($ + S?’) 2 ]*. Furthermore, for convenience 

in stating the main theorem, define the notation 

r 

A S7 = A + 7 ^ ^ aAASi. 

i - 1 

Theorem 5.3. Let n c < n and suppose there exist matrices P, Q, P, Q E N 1 


satisfying 

0 = Aj y P+PA sy + R i + Y^ [l~ 2 {o? i SjSi+p?PA i AjP)+7- 1 p- l \a i \Z i +P^i t ] 

i= 1 L 

-PEP+tJPEPtj., (5.23) 

0 = [A sl + ^r 2 fiAAj(P+P)]Q+Q[A sy +j2'y- 2 fiAAj(P+P)) T +V l 

i=l i=l 

-QEQ+r x QEQrJ, (5.24) 

0 = (A s ,-QY:+j2r 2 PlA t AjP) T P+P(A s ,-Q^+J2 1 - 2 P!A i AjP) 


i= 1 


Z= 1 


+^7" 2 ^ l 2 P/l i >l7F , +PSP-rJPEPr x , (5.25) 

1=1 

T T 

0 = (As,-VP+Y / r 2 PlWlP)Q+Q(A,-,-'ZP+Y,r‘ 1 ffA i AjP) T 


1=1 


i=l 


-{QT,Q - txQEQtJ , 


(5.26) 


rank Q = rank P = rank QP = ra c , (5.27) 

and let A c , B c , and C c be given by 

r 

Ac = r(A S7 - QE - EP + 7 -2 Y, P 2 AiAjP)G T , (5.28) 

1=1 

B c = TQC t V 2 -\ (5.29) 

C c = -P 2 - 1 B t PG t . (5.30) 
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Then (^4 + AA } E ) is detectable for all A A G U if and only if A+ A A is asymptotically 
stable for all AA € U. In this case, the worst-case % 2 performance criterion (5.7) 
satisfies the bound 


J(A C , B c ,C c ) < tx[PV l +Q{PEP-TlPi:PT L -r 2 Y, pA ^ P ^ ( 5 - 31 ) 

i= 1 

= tr \QR X + P(Q£Q - t l Q£QtJ ) 

~Q(f^l- 2 ^SjSi + r l 0r>i\Zi + ^ 2 ii)]. (5-32) 

*=i 

Proof. The proof is constructive in nature. We first obtain necessary conditions 
for the Optimization Problem and show by construction that these conditions serve as 
sufficient conditions for closed-loop stability. Specifically, it can be shown (see [13] for 
a similar construction) that the existence of P, Q, P, Q G N” satisfying (5.23)-(5.26) 
implies the existence of P G N" satisfying (5.18) where P is given by 


P + P -PG T ' 
-GP GPG T ' 


Now, the proof of robust stability and the upper bound on / Ri performance (5.5) for 
all uncertain perturbations AA G U follows from Theorem 5.2. 

Next, to optimize (5.20) subject to constraint (5.18) over the open set 

r 

5 = {( P,A c ,B c ,C c ) :PG r, A sl + r 2 Y, A ^P^ P 

i= I 

is asymptotically stable and (A c , B c ,C c ) is minimal} 


form the Lagrangian 


~ r IT ~ 


C(A c , Be, C c , Q, P, A) = tr [A PV + Q{<P + PA S , + R 

+ E + ffPAiAJP) + r'K>,\Zi + ‘/i] }] . (5.33) 


1=1 
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w 


v 


where the Lagrange multipliers A > 0 and Q € R n * n are not both zero. We thus 
obtain 


dC , 7 _ 2 

5? = (A7+t 


J^pfA,AjP)Q + < 3 ( 4 ,, 

i=l 


+ 7 - 2 £Mp) t + au. 

i=l 


(5.34) 


Setting fp = 0 yields 

0 = ( 4 , + r 2 E 0 2 A,AjP)Q + < 3 ( 4 , + r 2 E 0 2 A<AJ Pf + XV. (5.35) 

1=1 t = l 

Since A S7 + 7 -2 £)[ =1 ffiAiAjP is assumed to be asymptotically stable, setting A = 0 
implies Q = 0. Hence, it can be assumed without loss of generality that A = 1. 
Furthermore, Q is nonnegative definite. The remainder of the proof follows as in [13]. 
Briefly, the principal steps are as follows. 


Step 1. Compute and 

Step 2. Partition (5.18) and (5.35) into six equations (a)-(f) corresponding to the 
n x n, n x n c , and n c x n c sub-blocks of P and Q , respectively. Next, since 
the compensator triple ( A c ,B c ,C c ) is controllable and observable, using a 
minor extension of the result from [2] and Lemma 12.2 of [115], (c) implies 
that the lower-right n c x n c block of P is positive definite. Using similar 
arguments we can show that the lower-right n c x n c block of Q is positive 
definite. See [13] for details. 

Step 3. Form (b) times the n x n c sub-block of Q plus the n c x n c sub-block of Q 
times (c) to define the projection matrix r and the new variables P , Q , P, 
Q, G, and T. 
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Step 4. Use the result of Step 1 and Step 3 to solve for the compensator gains 
(5.28)-(5.30). 

Step 5. Manipulate (a), (b), (d), and (e) to yield (5.23)-(5.26). 

Step 6. Use the results of Step 3 to show that (5.20) is equivalent to (5.31). 

For a detailed exposition of a similar proof, see [13]. □ 

Remark 5.3. In the full-order case, set n c = n so that G = T = t = I n and 
Tx = 0. In this case the last term in each of (5.23)-(5.26) is zero and (5.26) is 
superfluous. 

Theorem 5.3 provides constructive sufficient conditions that yield dynamic feed- 
back gains A c , B c , and C c for robust stability and performance. When solving (5.23)- 
(5.26) numerically, the values of 7, Qj, Pi, and Si, i = l,...,r, can be adjusted to 
examine trade-offs between P 2 performance and robustness. As discussed in [54], to 
further reduce conservatism, one can view the scalars a* and /?, as free parameters 
and optimize the worst-case H 2 performance bound J = tr PV with respect to a t 
and The simplest case to consider is the case where Si is skew-symmetric or, 
equivalently, Z* = 0. In this case 

= 2&7“ 2 tr PA t S t Q + 2ai7" 2 tr SjS t Q = 0 (5.36) 

oai 

and 

% = 2a i7 - 2 tr PAiS t Q + 2^~hv PA.AjPQ - I t Q = 0, (5.37) 

dpi 

where Q satisfies 

0 = (i, 7 + 'Er*P?A i A7p)Q + «(^ 7 + 'pr'ftMAjpf + V. (5.38) 

i=l t=l 
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Furthermore, in this case, 



d " J 2 =2r 2 tvS t QSj>0, 

oai 

(5.39) 

d 2 J 

2 7 ~Hr AjPQPA, + 6/3“ 4 tr l',Q > 0, 

(5.40) 

d 2 J d 2 J l 

don 2 dpi 2 \ 

' fL ) = lfry- 2 /3f 4 tr UQtvSiQS! > 0, 
S, daidPiJ 

(5.41) 


which imply that (5.36) and (5.37) provide necessary conditions for a local minimum. 

In the case Z* ^ 0 we need to consider the cases c** = 0 and c*i ^ 0 since J is not 
differentiable at a* = 0. First, let a { = 0. In this case 

% = 2& 7 - 2 tr PA t AJPQ - 2/3-hr hQ = 0, (5.42) 

dpi 

where Q satisfies (5.38) with a* = 0. Next, consider the case where a* # 0. In this 
case 

^ - = 2A7 _2 tr PAM + 2ai7 _2 tr SjSiQ + 7 _1 yS i ‘ 1 sgn a^tr Z { Q = 0 (5.43) 

oai 

and 

= 2a i7 - 2 tr PAM + 2&7 _2 tr PAAJPQ - 7 - 1 £f 2 |a i |tr Z { Q - 2/T 3 tr hQ = 0 

OPi 

(5.44) 

where Q satisfies (5.38) and sgn <*» = By using (5.36) and (5.37) (or (5.43) and 
(5.44) for non-skew-symmetric Si) within a numerical optimization algorithm, the 
optimal robust reduced-order controllers and scaling parameters a*, i = 
can be determined simultaneously. 

5.6. Two- Mass Example 


Consider the dynamic system shown in Figure 5.1, which represents a flexible 
structure with uncertain high-frequency dynamics [40]. The equations of motion for 



this system are 


m 1 X 1 + C\±i - c 2 (x 2 - ii) + kiXi - k 2 (x 2 - xj) = u, 
m 2 x 2 + c 2 (x 2 - iq) 4- k 2 (x 2 - Xj) = 0. 


krl 


h 1 ' *1 K ; 


-A/\n- 

m=l 


3 

* 

II 

— 3 — 


— 30 — 

£ 


c,=0.0/ I ^ M 

I ^ J'eol ■*/ 


c 2=0.01 


4noncol *2 


Figure 5.1: Two-mass system 


We first consider the case of a colocated sensor and actuator pair, where the output 
is given by y co \ = iq. Letting mi — 1, m 2 — 10, k\ = k 2 = 1, and C\ — c 2 = 0.01 and 
transforming to real normal coordinates yields the plant state space realization 


0.0002 

0.2208 

0 

0 


' -0.1439 ' 

0.2208 

-0.0002 

0 

0 

X + 

0.2168 

0 

0 

-0.0103 

1.4320 

-0.0426 

0 

0 

-1.4320 

-0.0103 


1.1890 


y co] = [ -0.0545 0.0819 -0.0352 0.8181 ] x. 


As in [40], the matrices D j, D 2 , E\, and E 2 are chosen to be 


A = 


0 0 
1 0 
0 0 ’ 
0 0 


D 2 = [ 0 1 ] , 


Ei = 


10 0 0 
0 0 0 0 ’ 


Di = 


0 

1 ’ 


so that the LQG compensator places a notch at the second modal frequency. Uncer- 
tainty in the damped natural frequency of the second mode u>d 2 = 1.432 is modeled 
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by choosing 


0 0 0 0 

0 0 0 0 

Al ~ 0 0 0 1 ■ 

_ 0 0 -1 0 

The quasi-Newton optimization algorithm discussed in Section 2.3 was used to 
compute full-order controllers ( n c = n) that minimize the cost bound J for several 
values of 7 . The actual ?{ 2 cost was computed for a range of values of the damped 
natural frequency of the second mode for the LQG controller and for the implicit 
small gain (with S\ = A\ and aq and /?i obtained by (5.36) and (5.37)) and scaled 
Popov controllers [106] corresponding to 7 = 15, 7, and 2. The cost dependence is 
shown in Figure 5.2. As 7 decreases, the U 2 cost of the nominal closed-loop system 
increases while the Ti 2 cost of the perturbed closed-loop system remains near the 
nominal value for a larger range of perturbations. The LQG controller stabilizes the 
closed-loop system for only small perturbations in the damped natural frequency of 
the second mode, while the implicit small gain controllers stabilize the closed-loop sys- 
tem and provide performance close to the optimal level even for large perturbations. 
Hence, robust performance over a large range of the uncertain parameter is achieved 
for only a small increase in the H. 2 cost above the optimal. Also note that the robust- 
ness/performance tradeoffs of the implicit small gain controllers are comparable to 
those of the scaled Popov controllers which are obtained using frequency-dependent 
multipliers [106]. 

The frequency responses of the LQG controller and the implicit small gain con- 
trollers with 7 = 15, 7, and 2 are shown in Figure 5.3. The LQG controller is unstable 
and achieves closed-loop stability and nominal performance by placing a notch at the 
nominal damped natural frequency of the uncertain second mode. Hence, closed- 
loop performance degrades considerably with the damped natural frequency of the 
second mode is perturbed. The implicit small gain controller with 7 = 7 has only 



Figure 5.2: Dependence of the % 2 cost on the damped natural frequency of the 
second mode: Colocated case 

a shallow notch near the damped natural frequency of the second mode, while the 
controller with 7 = 2 has no notch near that frequency. Hence, these controllers sacri- 
fice nominal performance for improved robust performance over a larger range of the 
uncertain damped natural frequency. As 7 decreases, the controllers guarantee robust 
performance over a larger range of S. Note that the controller obtained with 7 = 2 
is positive real. Since the plant is a model of a flexible structure with a colocated 
sensor and actuator pair, it is also positive real, and thus the closed-loop system is 
asymptotically stable for all values of the uncertain damped natural frequency. 

Next we consider the two-mass system of Figure 5.1 with a noncolocated sensor 
and actuator pair by choosing 

2/noncoi = [ -0.1063 0.1597 0.0018 -0.0419 ] x. 

As in [40], we increase the matrix E\ by a factor of ten to enhance the notching 
characteristics of the LQG controller. 
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Figure 5.3: Frequency responses of implicit small gain controllers: Colocated case 

The implicit small gain synthesis technique was used as before to compute full- 
order controllers ( n c = n) that minimize J for a range of 7. The actual cost was 
computed for a range of the damped natural frequency for the LQG controller and 
for the implicit small gain and scaled Popov controllers [106] corresponding to 7 = 15 
and 4. The cost dependence is shown in Figure 5.4. In the noncolocated case, the 
LQG controller stabilizes the closed-loop system for a smaller range of the uncertain 
parameter than the unstable LQG controller for the colocated plant. 

The frequency responses of the LQG controller and the implicit small gain con- 
trollers with 7 = 15 and 7 = 4 were computed and are shown in Figure 5.5. Since the 
plant is not positive real, robust performance cannot be achieved simply by obtaining 
positive real controllers, as in the colocated case. Instead, as seen in Figure 5.5, the 
controllers widen the notch at the nominal frequency of the uncertain mode, and the 
controller with 7 = 15 deepens the notch as well. 
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Figure 5.4: Dependence of the %2 cost on the damped natural frequency of the 
second mode: Noncolocated case 



Figure 5.5: Frequency responses of implicit small gain controllers: Noncolocated 
case 


80 




5.7. Three-Mass Example 






V 


V 




V 


V 




Consider the three-mass, two-spring system shown in Figure 5.6 with mi = m? = 
m 3 = 1 and an uncertain spring stiffness &2 [52, 106]. A control force acts on mass 
3 while the position and velocity of mass 1 are measured resulting in a noncolocated 
control problem. The nominal dynamics, with state variables defined in Figure 5.6, 
are given by 


W 




*2 



M= 1 


-AAAi ^ =1 KWH ^ =1 


Figure 5.6: Three-mass system 
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c = 

[1 

1 0 0 

0 
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]. 

D 2 = [ 0 

1] 

| 5 




and = A; 2 nom — 1. The actual spring stiffness of the second spring can be written 


as k 2 = & 2 nom + S so that the actual dynamics are given by A A = A + 5A\, where 


M = 


0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

0-1 10 0 0 
0 1-1000 


Furthermore, let 


E x 


1 1 0 0 0 0 
0 0 0 0 0 0 ’ 


E 2 = 


0 

1 




81 





82 






A full-order (n c = n ) implicit small gain (with Si = A\ and a\ and /?i obtained by 
(5.43) and (5.44)) compensator was designed with a value of 7 = 20. The actual %2 
cost was computed for the LQG controller and for the implicit small gain and scaled 
Popov controllers corresponding to 7 = 20 for a range of the uncertain parameter 5. 
The cost dependence is shown in Figure 5.7. Though a value of 7 = 20 corresponds to 
a parameter uncertainty of <5 = 0.05, it is seen that the implicit small gain controller 
is robust over a much larger range. Finally, the frequency responses of the LQG 
controller and the implicit small gain controller are shown in Figure 5.8. 

5.8. Conclusion 

This research extended the implicit small gain guaranteed cost bound [54] to con- 
troller synthesis. Specifically, the implicit small gain guaranteed cost bound was used 
to address the problem of robust stability and W 2 performance via fixed-order dy- 
namic compensation. A quasi-Newton optimization algorithm was used to obtain 
robust controllers for several illustrative examples. The design examples considered 
demonstrated the effectiveness of the implicit small gain guaranteed cost bound. Fi- 
nally, we note that the conservatism of the proposed implicit small gain guaranteed 
cost bound is difficult to predict and will depend upon the actual value of P deter- 
mined by solving (5.18). 




CHAPTER 6 


Robust Resilient Dynamic 
Controllers for Systems with 
Parametric Uncertainty and 
Controller Gain Variations 

6.1. Introduction 

It is well known that unavoidable discrepancies between mathematical models 
and real-world systems can result in the degradation of control-system performance 
including instability. Thus it is not surprising that a considerable amount of re- 
search over the past two decades has concentrated on analysis and synthesis of feed- 
back controllers that guarantee robustness with respect to system uncertainties in 
the design model (see [27] and the numerous references therein). These robust con- 
troller synthesis frameworks include the Youla parameterization of all stabilizing con- 
trollers [110], "H 2 and Tioo (including desired weighting functions for loop shaping) 
synthesis [5,39,126], C\ control design [25], //-synthesis for structured real and com- 
plex uncertainty [126], and robust fixed-structure controller synthesis [16]. Almost all 
of these techniques yield very high order controllers in relation to the original system 
order. A notable exception is the fixed-structure controller design methodology [16] 


which directly accounts for controller complexity constraints, including controller or- 
der, within the control-system design process. However, an implicit assumption in- 
herent in all of the above mentioned design frameworks is that the resulting robust 
controller will be implemented exactly. But in most applications (and, in particular, 
aerospace applications), reduction in size and cost of digital control hardware results 
in limitations in available computer memory and wordlength capabilities of the digital 
processor and the A/D and D/A converters. This further results in roundoff errors 
in numerical computations leading to controller implementation imprecision. Hence, 
any controller that is part of a feedback system must be insensitive to some amount 
of error with respect to its gains. 

Within the context of robust controller synthesis, the above issues were first 
pointed out in the enlightening and very interesting paper entitled, “Robust, Fragile, 
or Optimal?” [75]. Specifically, the authors in [75] very elegantly point out that the 
powerful (weighted) Hi, (weighted) "Hoc, £ 1 , and fi controller design approaches, even 
though quite robust with respect to system uncertainty, are surprisingly very sensi- 
tive with respect to errors in the controller coefficients resulting in vanishingly small 
stability margins. Of course, since the control system is part of the overall closed-loop 
system, the authors of [75] show through a series of examples that most of the elegant 
multivariable robust control frameworks discussed above destabilize the closed-loop 
system for extremely small perturbations in the controller coefficients. Hence, even 
though these controllers are robust (with respect to plant uncertainty) and in some 
cases optimal, they are extremely fragilel This further implies that the resulting con- 
trollers preclude the control-system designer from tuning the controller gains around 
a designed nominal controller which has been the creed of practicing control engi- 
neers to capture performance requirements which are not directly addressed within 
the original design problem. Finally, it is interesting to note that numerical experi- 


86 



ments seem to indicate that the fragility or brittleness of the controller is exacerbated 
with increasing controller order. 

In this chapter, the robust fixed-structure guaranteed cost controller synthesis 
framework of [10, 13, 52] for systems with structured parametric uncertainty is ex- 
tended to address the design of non-fragile or robust resilient fixed-order (i.e., full- 
and reduced-order) dynamic compensation. For flexibility in controller synthesis, we 
adopt the approach of fixed-structure controller design which allows consideration 
of arbitrary controller structure, including order, internal structure, and decentral- 
ization. Specifically, using quadratic Lyapunov bounds, a rigorous development of 
sufficient conditions for robust stability and worst-case %2 performance via fixed- 
order dynamic compensation is presented for uncertain feedback systems wherein the 
controller can tolerate multiplicative or additive gain variations with respect to its 
nominal coefficients. These sufficient conditions are in the form of a coupled system 
of algebraic Riccati equations that characterize robust resilient reduced-order con- 
trollers. Hence, the proposed robust resilient controllers guarantee robust stability 
and robust performance in the face of both system uncertainty and controller er- 
rors. The proposed approach is applied on several numerical examples which clearly 
demonstrate the need for robust resilient control. 

6.2. Robust Stability and Performance 

In this section we state the robust stability and performance problem. This prob- 
lem involves a set U C R nxn of uncertain plant perturbations A A of the nominal dy- 
namics A and a set U c C R ncXn ° x R” cXi x R mxnc of uncertain controller perturbations 
(AA C , AB C , AC c ) of the nominal controller gain matrices (A c , B c ,C c ). The objective 
of this problem is to determine a fixed-order, strictly proper dynamic compensator 
(A c , B c ,C c ) that stabilizes the plant for all variations in U x U c and minimizes the 


worst-case H 2 norm of the closed-loop system. In this section and the following sec- 
tion no explicit structure is assumed for the elements of U x U c . In Section 6.4 and 
Section 6.7, two specific structures of the variations in U x U c will be introduced. 

Robust Stability and Performance Problem. Given the n th -order stabiliz- 
able and detectable uncertain system 


x(t) = (A + AA)x(t) + Bu(t) + Diiv(t), t G [0, oo), (6.1) 

y(t) = Cx(t) + D 2 w(t), (6.2) 


determine an n‘ h -order robust resilient dynamic compensator ( A c ,B c ,C c ) such that 
the closed-loop system consisting of (6.1), (6.2), and controller dynamics 


i c (t) = (A c + AA c )x c (t) + (B c + AB c )y(t), (6.3) 

u(t) = (C c + AC c )x c (t), (6.4) 


is asymptotically stable for all allowable plant uncertainties and controller gain vari- 
ations (A.4, A^c, A B c , AC c ) G U x U c and the performance criterion 


J{A C , B c ,C c ) = sup 

(AA,AAcABc,&Cc)eUxU c 


1 

limsup -E / [x T (s)f?ia:(s) + u T (s)R. 2 u(s)] ds, 
t — KX) t J 0 

(6.5) 


is minimized. 

For each uncertain plant and controller variation (A.A, A^4 C , A B c , AC c ) G U x U c , 
the closed-loop system (6.1)-(6.4) can be written as 


x(t) = (A + A A)x(t) + (D + AD)w(t), t G [0, oo), (6.6) 


where 


A 

AA 

BACr 1 

A 

0 


A B c C 

i 

u 

<1 

, AD = 

ab c d 2 
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and where the closed-loop disturbance (D + AD)w(t) has nonnegative definite inten- 
sity 


V A = ( D+AD){D+AD) T = 


Vi 0 

0 B c V 2 Bj + A B c V 2 Bj + B c V 2 ABj + AB c V 2 AB c T j 


6.3. Sufficient Conditions for Robust Stability and 
Performance 

In practice, steady-state performance is only of interest when the undisturbed 
closed-loop system is robustly stable over U x U c . The following result is immediate. 
For convenience, define 


~ A 


R A = {E + AE) t {E+AE) = 
where 


Ri 0 

0 CjR 2 C c + ACjR 2 C c + Cj R 2 AC c + A Cj R 2 AC c 


A E = [0 E 2 AC c ] . 


Lemma 6.1. Let ( A c , B c , C c ) be given and assume that ^4 -l- A^4 is asymptotically 
stable for all plant and controller gain variations (AA, AA C , AB C , AC c ) € UxU c . 
Then 


J(A c ,B c ,C c )= sup tr Q a R a , 

(AA,AAcABcAC c )eUxUc 

where Q A € N" , is the unique nonnegative definite solution to 

0 = (i + AA)Q a + Q a (A + A A) T + V A . 


(6.7) 


( 6 . 8 ) 


The key step in guaranteeing robust stability and performance is to bound the 
uncertain terms R A in the cost function (6.7) and AAQ a + <5aA4 t and V A in the 
Lyapunov equation (6.8) by smooth bounding functions £)*(•), i = 1,2,3. For the 
statement of the next result, define the notation 


AR = R a -R = 


0 0 

0 A CTR 2 C c + CjR 2 AC c + ACJR 2 AC c 
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AV = V A - V = 


0 0 
0 A B c V 2 Bj + B c V 2 ABj + AB c V 2 ABJ 


Theorem 6.1. Let A c 6 R ncXnc , B c € R ncX ', and C c € R mx " c be given and let 

fi,: FfX R*c*n c x R*cxZ x R mxn c gn q* . R n c xi gn an( J : R mXric -► S fl , be 

such that, for (AA, A/l c , A B c , A C c ) eU xU c , and Q £ N”, 

AAQ + QAA^ < fli {Q, A c , B c , C c ), 

AV < Q 2 (B c ), 

A R < n 3 (C c ), 


and suppose there exists Q € N" satisfying 

0 = AQ + QA T + n AQ, A c , B c , C c ) + V + D 2 (B C ). 


(6.9) 

( 6 . 10 ) 

( 6 . 11 ) 


( 6 . 12 ) 


Then 

(A + AA, D + AD) is stabilizable for all (Ai4, AA C , A B c , AC c ) eWx U c , (6.13) 

if and only if A + A^4 is asymptotically stable for all (A,4, AA C , A B c , A C c ) 6 U x U c . 
In this case, 

Qa<Q , (AA,AA c ,AB c ,AC c )eUxU c , (6.14) 

where Qa is given by (6.8), and 

J(A C , J5 C , C c ) < J{A C , B c , C c ) = tr Q [r + fi 3 (C c )] . (6.15) 

Proof. We stress that in (6.9), Q denotes an arbitrary element of N”, whereas 
in (6.12), Q denotes a specific solution of the modified Lyapunov equation (6.12). 
This minor abuse of notation considerably simplifies the presentation. Now, for 
(AA,AAc,AB c , AC c ) € ZY xU c , (6.12) is equivalent to 

0 = (A + A A)Q + Q{A + A A) T 

+Qi(Q, A c , Be, Cc) - (A AQ + QAi T ) + V A + D 2 (B C ) - AV. (6.16) 
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Hence, by assumption, (6.16) has a solution Q € N* for all (/\A, AA C , AB C , AC c ) € 
U x U c and, by (6.9) and (6.10), fii(Q, A c , B c ,C c ) — (AAQ+QAA t ) and 0 2 {B C ) — AV 
are nonnegative definite. Now, if the stabilizability condition (6.13) holds for all plant 
uncertainties and controller gain variations (AA, AA C , A B c , A C c ) Gl/x U c , it follows 
from Theorem 3.6 of [115] that 

(i + A A, [V A + f A c , B c , C c ) - (A AQ + QAA t ) + Q 2 (B C ) - AF]*) 

is stabilizable for all (AA, AA C , AB c ,AC c ) eU xU c . It now follows from (6.16) and 
Lemma 12.2 of [115] that A + A A is asymptotically stable for all plant uncertainties 
and controller gain variations (AA, AA c , A B c , AC c ) GWx U c . Conversely, if A + AA 
is asymptotically stable for all (AA, AA C , AB C , AC c ) € li xU c , then (6.13) holds. 
Next, subtracting (6.8) from (6.16) yields 

0 = {A + AA)(Q - Qa ) + (Q — Qa)(A + AA) T 

+Qi (Q, A c , B c , C c ) - (A AQ + QAA T ) + f) 2 (B c ) - AV, 

or, equivalently, since A + AA is asymptotically stable for all (AA, AA C , A B c , AC c ) € 
U x U c , 

Q -Q a = (Q, A. Bo Cc) - ( A AJ+QA/T) +n 2 (B c ) - AV']e^ +A ^) T ‘ d! 
> 0 . 

which implies (6.14). The performance bound (6.15) is now an immediate consequence 
of (6.14) and (6.11) by noting that 

J(A c ,B c ,C c ) = sup trQA^A < sup trQB A = sup trQ(R+ AR) < tvQ(R+Q 3 (C c )). 

AAeUxUc AAeUxUc AAeUxUc 


Remark 6.1. In applying Theorem 6.1 it may be convenient to replace Condition 


(6.13) with the stronger condition 


A + AA, 


„ ^ i l/2\ 

+ fix (Q, A c , B c , C c ) - (A AQ + QAA T ) + Q 2 (B C ) - A Vj ) 

is stabilizable for all (AA, AA c , A B c , AC c ) 6Wx U c , (6.17) 


which is easier to verify in practice. Clearly, (6.17) is satisfied if ( Q , A c , B c , C c ) 

— (AAQ + QAA T ) + Cl 2 {B c ) — AV is positive definite for all plant uncertainties and 
controller gain variations (AA, AA C , A5 C , AC c ) G U x U c - This will be the case, for 
example, if either V is positive definite or strict inequalities hold in (6.9) or (6.10). 

6.4. Multiplicative Controller Uncertainty Struc- 
ture and Guaranteed Cost Bound 


Having established the theoretical basis for our approach, we now assign an explicit 
structure to the sets U and U c and the bounding functions fii(Q, A c , B c , C c ), Q 2 (B C ), 
and fi 3 (C c ). Specifically, the uncertainty set U capturing parametric plant uncertainty 
and the uncertainty set Uc capturing multiplicative (relative) controller gain variations 
are defined by 

U = {AA : A A = B 0 FC 0 , F r F < 7 ~ 2 / r }, (6.18) 

U c = {(AA C , AR C , AC c ) ; AA C = SA C , A B c = 6B C , A C c = 5C c , |<$| < 7 C _1 },(6-19) 

where Bo G R nxs , Cq G R rxn are fixed matrices denoting the structure of the plant 
uncertainty, F G R* xr is an uncertain matrix, 5 is an uncertain real parameter, and 7 , 
7 C are given positive numbers. With this uncertainty characterization, the closed-loop 
system (6.6) has structured uncertainty of the form 

A A = B 1 FC 1 + 5B 2 C 2 , (6.20) 
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where 


V 









0 

A, 1 

Bi = 

' Bo ' 
0 

, C x = [ Co 0 ] , B 2 = 

0 OB' 
. In c In c 0 

,c^ 

7 *. 

o 

1 


. (6.21) 


For the structure of U and U c as specified by (6.18) and (6.19), the bounding functions 
f^l(Q) ^C! B c , C c ), fi 2 (B c ), and fi 3 (C c ) can now be given a concrete form. For the 
statement of the next result, define R Cl = and V Cl = D Cl Dj y , where 


B Cl = [ 0 E 2 C c ] , D Cl = 


0 

B c D 2 


Proposition 6.1. Let a, ct c > 0. Furthermore, let U and U c be defined by (6.18) 
and (6.19), respectively, and define fii(-), fl 2 (-), and fi 3 (-) by 


f^i(Q) A c , B c , C c ) 
fi 2 (B c ) 
ft 3 (C c ) 


aB x Bj + o c B 2 B 2 t + Q{oT l ')~ 2 C x C\ + a; 1 % 2 CjC 2 )Q, (6.22) 

( 7c - 2 + 2 7 c - 1 )K 1 , (6-23) 

( 7c - 2 + 2 7 c " 1 )4 1 . (6-24) 


Then (6.9)— (6. 1 1) are satisfied. 
Proof. Note that 


0 < 


a2B x - a-iQCjF r 


1 _ _ 


q2 B x -a~2 QCjF T 


r 1 


+ 


i r i 


Oi 2 B 2 - etc 2 6QC 2 


-i T 


a 2 B 2 - etc 2 SQC 2 


< ctB x Bj + ct c B 2 Bj + Q(a~ 1 y~ 2 CjCi + ot^% 2 CjC 2 )Q - (A AQ + QAA T ), 

which proves (6.9) with U and$ c given by (6.18) and (6.19), respectively. Next, note 
that 


AV = 


0 0 
0 {8 2 + 2S)B c V 2 BJ 


< 


0 0 
o (y c ~ 2 + 2%')B c V,Bj 


= fl 2 (B c ), 
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which proves (6.10) with U c given by (6.19). Finally, a similar construction proves 
( 6 . 11 ). □ 
Next, using Theorem 6.1 and Proposition 6.1 we have the following immediate 
result. 


Theorem 6.2. Let ct, a c > 0, and suppose there exists a matrix QgN" satisfy- 
ing 


0 = AQ + QA T + Q{a~ x 'y~ 2 CjC\ + a; l j; 2 CjC 2 )Q + V 

+aB\Bj + a c B 2 Bj + (7 C 2 + 2j c 1 )Kn* (6.25) 


Then (A + A A, D + A D) is stabilizable for all (A A, AA C , A B c , AC c ) € U x U c if and 
only if A + A^4 is asymptotically stable for all (AA, AA C , AB C , AC c ) eWx U c . In 
this case, 

Qa < Q, (A A, AA C , AB c ,AC c ) eUxU c , (6.26) 

where Qa satisfies (6.8), and 


J(Ac, B c , C c ) < tr Q R + ( 7 C 2 + 27 c 


(6.27) 


6.5. Decentralized Static Output Feedback Formu- 
lation 


In this section we use the fixed-structure control framework discussed in Chapter 2 
to transform the Robust Stability and Performance Problem to a decentralized static 
output feedback setting. Specifically, note that for every dynamic controller (6.3), 
(6.4) with gain variations (AA C , A B Ci A C c ) 6 U c given by (6.19), the uncertain closed- 
loop system (6.1)-(6.4) can be written as 


' x(t) ' 


x c (t) 



A + AA (1 + S)BC c 
(1 + S)B c C (1 + S)A C 


x{t) 

x c (t) 


+ 


Dx 

(1 + 6)B c D 2 


w(t). (6.28) 
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Furthermore, by treating A c , B c , and C c as decentralized static output feedback 
gains with multiplicative (relative) uncertainty as shown in Figure 6.1, we can pull 
the uncertainty into the plant model obtaining 



Figure 6.1: Decentralized static output feedback: Multiplicative controller uncer- 
tainty 


k(t) = {A + B l FC 1 )x(t) + (1 + 5)^2 $*«»(*) + B*w(0> t e [0,oo), (6.29) 


1 


yi(t) = C yi x{t ) + V yWi w{t), i = 1,2,3, 

3 

z(t ) = C z x(t) + (1 + (5) ^2 Vzu.Uiit), 


(6.30) 

(6.31) 


i=l 


where Ui(t), i = 1, 2, 3, are given by 


ui(t) = A c y!{t), u 2 (t ) = B c y 2 {t ), u 3 (t) = C c y 3 {t). 


(6.32) 


Finally, by rewriting the decentralized control signals in the compact form 


u(t) - ICy{t), 


(6.33) 






where 


and K, is given by 



Wl(0 


y i(0 


_ u 3 (t) _ 

<11 

CO 

1 



0 0 
0 B c 0 
0 0 C c 


the uncertain closed-loop system is given by 

i(t) = + A/i) x(t) + (f) 4- A Dj w(t), 

z(t ) = (^D + AE'j x(t), 


where 


(6.34) 


(6.35) 


(6.36) 

(6.37) 


A A = B X FC\ + 5B u ICCy, A D = 6B u JCV yw , AE = 5V zu KC y . 


Note that B 2 — B u and C 2 — JCC y . 

We can now recast the Robust Stability and Performance Problem as the following 
Auxiliary Optimization Problem. 

Auxiliary Optimization Problem. For given a , cv c > 0, determine the block- 
diagonal controller matrix 1C G R( 2n c+m)x(2n c +() t h a t minimizes 


J{K) = tr Q[R + ( 7c - 2 + 2 7 - 1 )R Cl ] , 


(6.38) 


where Q G N” satisfies (6.25). 

It follows from Theorem 6.2 that the satisfaction of (6.25) for <2 G N” along with 
the generic stabilizability condition (A + AA, D + AD) leads to closed-loop robust 
stability along with robust , h 2 performance. 
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6.6. Sufficient Conditions for Fixed-Order Resil- 
ient Compensation with Multiplicative Un- 
certainty 

In this section we state sufficient conditions for characterizing dynamic output 
feedback controllers guaranteeing robust stability and robust H 2 performance with 
respect to system plant uncertainty and multiplicative controller gain variations. 

Theorem 6.3. Let a, a c > 0. Suppose there exist matrices Q, P € N” satisfying 



0 = AQ+QA t +Q (a~ l r 2 CjC x +a: l % 2 CjC 2 ) Q+V 


V 

AaBiBj +a c B 2 Bj + (7,. 2 +2% ') V Cl , 

(6.39) 

V 

0 = [A+Q( K a- l r 2 CjC 1 +a; l ^ 2 CjC 2 ')] T P 

+P [A+Q (a- l r 2 CjCi+a;h c 2 CJC 2 )] +R+{ 7c“ 2 +27c _1 ) 

(6.40) 


and let (A c ,B c ,C c ) satisfy 



0 = BlPQCl + a~ l 'y~ 2 A c C y QPQCj l , 

(6.41) 

V 

0 = BlPDVl 2 + BlPQCl +a;^B c C v QPQCl 



+ ( 7e - 2 + 2 7c -‘) BlPB u KV yw Vl^, 

(6.42) 

w 

0 = BlPQCl + VlMCl + ^'ifC^QPQCl 



+ ( 7c - 2 + 2 7c -') Vl,V,.KC,QC^. 

(6.43) 


Then (^4 + AA, D + A D) is stabilizable for all (A.4, AA c , A B c , AC c ) € U x U c if and 
only if A + AA is asymptotically stable for all (A.A, AA c , AB c , AC c ) € U xU c . In 
this case, the worst-case 7{ 2 performance of the closed-loop system (6.7) satisfies the 
bound 

J(A c ,B c ,C c ) < trQ[fl+ ( 7 - 2 +2 7 -‘) 4,] ■ (6.44) 
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Proof. First we obtain necessary conditions for the Auxiliary Optimization Prob- 
lem and then show, by construction, that these conditions serve as sufficient conditions 
for closed-loop stability and robust %% performance. Thus, to optimize (6.38) subject 
to (6.25), form the Lagrangian 

C(K, P, A) = tr {A [QR + ( 7c ~ 2 + 27c" 1 ) QRc] 

+P [AQ + QA t + Q (a-'r'cjc, + Q 

+V + aB,Bj + a c B 2 Bj + ( 7c ~ 2 + 2^) V Cl ] } , 

where the Lagrange multipliers A > 0 and P € are not both zero. By viewing 
K and Q as independent variables, we obtain 

Jg = [i+Q(a-S- 2 C?'C 1 +a c - 1 7 c - 2 C 2 T C' 2 )] T p 

+a[r+('7 2 +27 c ‘) Pc,]- (6.45) 

If A+Q (a-V’^Ci + a c 1 'y'c 7*^2 Cij is Hurwitz, then A = 0 implies P = 0. Hence, 
it can be assumed without loss of generality that A = 1. Furthermore, note that P is 
nonnegative definite. Thus the stationary conditions with A = 1 are given by 

= [i + Q(o- l 7- 2 C?'<J,+ac- 1 7c- 2 C?C 2 )] T P 

+P ].4 -r Q (a~'~f 2 C^C\ + 7; 1 yfCj C‘/j ] + P + (7 C 2 + 27 c 1 ) Pc, = 0, 

~ = Bl PQCl + a :'^A c C,QPQCl = 0 , 

H = B T „PDVl 2 + BlPQCl + a -'% 2 B c C,QPQCl 
+ (7c“ 2 + 27c' 1 ) BlPBuKV^V^ = 0, 

H = BlPQCl+VlMCl+c*;'% 2 C c C y QPQCl 

+ (if + 27c- 1 ) Vl,t> m KX,QCl = 0, 

which are equivalent to (6.40)-(6.43). Equation (6.39) is a restatement of (6.25). It 
now follows from Theorem 6.2 that the stabilizability condition (A + A A, D -I- AD) 


98 



for all (A.A, A A c , A B c , A C c ) G U x U c is equivalent to the asymptotic stability of 
A + A A for all (AA, A A c , A B c , A C c ) G U x U c . Finally, the Hi performance bound 
(6.44) is a restatement of (6.38). O 

Equations (6.39)-(6.43) provide constructive sufficient conditions that yield dy- 
namic controllers for robust resilient fixed-order (i.e., full- and reduced-order) output 
feedback compensation. In the design equations (6.40)-(6.42), one can view a and 
a c as free parameters and optimize the performance criterion (6.38) with respect to 


a and a c . In particular, setting = 0 and = 0 yields 

1 tr PQCj C X Q~\ * _ 1 [tr PQCJC 2 qV 

tr PBiBj ’ a ° ~ 7c tr PB 2 Bj 


(6.46) 


It is important to note that a and ot c given by (6.46) are implicit since Q and P are 
functions of a and a c . However, the optimal robust reduced-order controller gains 
and the scaling parameters a and a c can be determined simultaneously within a 
numerical optimization algorithm using ^ = 0 and J-^ = 0. For details of this fact, 
see Section 6.10. 


6.7. Additive Controller Uncertainty and Guaran- 
teed Cost Bound 


In this section we assign a different structure to the uncertainty set U. c and conse- 
quently the bounding functions = 1,2,3. Specifically, the uncertainty set U c 

is assumed to be of the form 

U c = {(A Ac, A B c , A C c ) : AA C = Sl Ac , A B c = Sl Bc , A C c = 5X Cc , \S\ < % 1 }, (6.47) 

where Ta c ,Zb c > an d Xc c are ones matrices of dimension R n c xn ' j R n « x/ , and R mXTlc , re- 
spectively, 5 is an uncertain real parameter, and 7 C is a given positive number. Note 
that, unlike the multiplicative uncertainty characterization addressed in Section 6.4, 
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the uncertainty characterization given by (6.47) can capture controller gain varia- 
tions with zero entries in the nominal gain matrices (A c , B c ,C c ). With this additive 
(absolute) uncertainty characterization, the closed-loop system (6.6) has structured 
uncertainty of the form 

A A = B 0 FCo, (6.48) 


where 


B 0 = 


0 I, 


H c 


0 B 
Z Bc 0 


F = 


' F 

0 

0 

0 ' 


'Co 

0 ' 

0 

SI n c 

0 

0 

Cn — 

0 

Za, 

0 

0 

Sit 

0 

> ^0 — 

c 

0 

0 

0 

0 

SIrn 


0 

1 


(6.49) 


For the structure of U x U c as specified by (6.18) and (6.47), the bounding functions 
fl 2 (B c ), and ^(Cc) can now be given a concrete form. For the statement of 


A ^ 


the next result, define R C2 = Ej E C2 and V C2 = Dj 2 D C2 , where 


K = [ o e 2 i Cc ] , d C2 


0 

Fb c E>2 


Furthermore, to enforce the block-structure of the uncertainty matrix F, define the 
set of compatible scaling matrices V by 


V = {Do > 0 : FD 0 ~ D 0 F, F T F < N}, 


where 

' 7 2 E 0 

0 7c^(«c +i+m) 

The condition FD 0 = D 0 F in V is analogous to the commuting assumption 
between the D-scales and A blocks in /z-analysis and synthesis which accounts for the 
structure in the uncertainty F. It is easy to see that there always exists such a matrix 
D 0 even if F is neither diagonal nor symmetric. For example, if F = fl r , where / is a 
scalar uncertainty, then D 0 can be an arbitrary positive definite matrix. Alternatively, 
if F € K rxr is nondiagonal, then one can always choose D 0 = block-diag[d/ r , Dq], 
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where d is a scalar and D 0 G R("c+f+m)x(n c +<+m) is an arbitrary matrix. Of course, F 
and Do may have more intricate structure, for example, they may be block-diagonal 
with commuting blocks situated on the diagonal. 

Proposition 6.2. Let oc\ > 0, a? > 0, and let Dq G T>. Furthermore, let U and 
U c be defined by (6.18) and (6.47), respectively, and define Qi(-), an d ^M') by 

fli(Q) = QCjDoNDoCoQ + B q Do 2 BJ, (6.50) 

n 2 (B c ) = ct 1 l % 2 V Cl + («! + 7c” 2 )K 2 , (6-51) 

th(C c ) = «2 J 7c r 2 ^c, + («2 + 7c- 2 )^. (6.52) 

Then (6.9)— (6. 11) are satisfied. 

Proof. Note that with D 0 G V 

0 < [b 0 Do 1 - QCjD 0 F T ] [Po^o -1 - QCjD 0 F T Y 

< B 0 Do 2 Bj + QCjDoNDoCoQ - (A AQ + QAA T ), 

which proves (6.9) with U and U c given by (6.18) and (6.47), respectively. Next, note 
that 

AV = 6D Cl £>l + SD C2 Dl + S 2 D C2 Dj 2 . 

Now, since 5 2 D C2 Dj 2 < 7 ~ 2 D C2 Dj 2 and 

■ 1 _! _ I r I _ 1 _ i T 

0 < a 2 D Cl ~ «i 2 A D Cl a 2 D C2 - 2 AD Cl 

< «i D«DZ + aT l % 2 D a Dl - S(D c ,t>l + 4.6J), 

it follows that 

AV” < (oi + 7 C 2 )V C2 + a, X 7 C 2 V Cl , 

which proves (6.10) with U c given by (6.47). Finally, a similar construction proves 

( 6 . 11 ). □ 
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Next, using Theorem 6.1 and Proposition 6.2 we have the following immediate 


result. 


Theorem 6.4. Let , a 2 > 0, D 0 6 V, and suppose there exists a matrix Q € N* 
satisfying 

0 = (6.53) 


Then ( A + A A, D + AD) is stabilizable for all (A A, A A c , A B c , A C c ) G U x U c if and 
only if A + AA is asymptotically stable for all (AA, AA C , AJ3 c , AC c ) 6 Wx U c . In 
this case, 

Qa<Q , {AA,AA c ,AB c ,AC c )eUxU c , (6.54) 

where Qa satisfies (6.8), and 


J(A C , B c ,C c ) < tr Q 


R + c%2 1 


7 C 2 


Rci + («2 + 7c 2 )^C 2 


(6.55) 


6.8. Decentralized Static Output Feedback Formu- 
lation 


As in Section 6.5, note that for every dynamic controller (6.3), (6.4) with gain 
variations (AA C , AB C , AC c ) € U c given by (6.47), the closed-loop system (6. 1)— (6.4) 
can be written as 


' x{t)' 


■ A + B 0 FC 0 

BC c + SBl Cc ' 

' x{t ) ' 

1 

L>i 

x c (t)_ 


B c C + 5I Bc C 

A c 4- 5Xj 4 C 

X c (t) 

r 

B C D2 + 5Ib c D2 


w(t). (6.56) 

Furthermore, by treating A c , B c , and C c as decentralized static output feedback gains 
with additive uncertainty as shown in Figure 6.2, we can pull the uncertainty into 
the plant model obtaining 


x(t) = (A + A A)x(t) + y ^B Ui Uj(t) + (B w + A t3 w )w(t), <g[0,oo), (6.57) 


i= 1 
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Figure 6.2: Decentralized static output feedback: Additive controller uncertainty 


Vi(t) = C yi x(t)+V yWi w(t) r i = 1,2,3, 

3 

z(t) = ( C z + A C z )x(t) + 'D ZUi Uj{i) i 

i= i 

where Uj(t), i — 1, 2,3, are given by (6.32) and 


AA 


BqFC 0 SB1 Cc 
51 Bc C 61 Ac J’ 


AB W 


A 


0 

SI Bc D 2 J ’ 


(6.58) 

(6.59) 


AC* = [ 0 5E 2 1 Cc ] • 


The uncertain closed-loop system is now given by 

£(t) = fA + AA) x(t) +(d + ad) w(t), t e [0, oo), (6.60) 

z{t) = (E + AD) x(t), (6.61) 

where AA = A A, AD = A B w , and A E = AC*. 

Now, as in the multiplicative controller uncertainty case, we introduce an Auxiliary 
Optimization Problem by considering 


J(IC) = tr Q R + a 2 2 R Cl + (of 2 + 7 C 2 )^c 2 ] > 


(6.62) 
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with Q e N” satisfying (6.53), and proceed by determining controller gains that 
minimize J(fC). 


6.9. Sufficient Conditions for Fixed-Order Resil- 
ient Compensation with Additive Uncertainty 

In this section we state sufficient conditions for characterizing dynamic output 
feedback controllers guaranteeing robust stability and robust performance with 
respect to system plant uncertainty and additive controller gain variations. 

Theorem 6.5. Let a 1? a 2 > 0, and let D 0 e V. Suppose there exist matrices Q , 
PeN" satisfying 

0 = (6.63) 


0 = (A + QCjD 0 ND 0 C o y P + p(A + QCjD 0 ND 0 Co) 

+R + c*2 1 % 2 R Cl + (02 + 7c 2 )-Rc 2 i , (6.64) 

and let ( A c , B c , C c ) satisfy 

0 = BlPQCl, ( 6 . 65 ) 

0 = BlPDVl a + BlPQCl + aT'i?BlPB»ICV m Vl 2 , (6.66) 

0 = Bl,PQCl + VlJQCl + a;'% 2 Vl 3 V,„t:C,QCl. ( 6 . 67 ) 


Then (A -f AA y D + AD) is stabilizable for all (A.4, AA C , A B c , A C c ) € U x U c if and 
only if A + AA is asymptotically stable for all (AA, AA C , AB C) AC c ) € U x ZY C . In 
this case, the worst-case performance of the closed-loop system (6.7) satisfies the 
bound v ' : 

J{Ac, B c , C c ) < tr Q [£ + + (a 2 + 7 c “ 2 )^c 2 ] • (6-68) 

Proof. The proof is similar to the proof of Theorem 6.1. □ 
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As in Section 6.6, one can view a^i = 1,2, and D 0 as free parameters and optimize 
the performance criterion J{K) given by (6.62) with respect to Ofi, * = 1,2, and Dq. 
In particular, setting fg = 0, i = 1,2, and ^ = 0 yields, respectively, 


_1_ trPVc, 2 
Tc tr PV Ci 



tr QR^ 
tr QRc 2 


1 

2 


(6.69) 


0 = NDoCoQPQCj - Dq'BZPBoDq 2 . (6.70) 


6.10. Quasi-Newton Optimization Algorithm 

As mentioned in Section 2.3, the optimization algorithm was initialized with an 
LQG controller for full-order controller designs, while for reduced-order control, the 
algorithm was initialized with a balanced truncated LQG controller. A large value 
was chosen for 7 and then a feasible value of 7 C was calculated. The quasi-Newton 
optimization algorithm was used to find the controller gains A c , B c , and C c . After each 
iteration, 7 C was decreased and the current values of the controller gains (A c , B c , C c ) 
were then used as the starting point for the next iteration. When 7 C could not be 
decreased any further, 7 was decreased and a feasible j c was calculated for the new 
value of 7 and the process was repeated. 


6.11. Second- Order Unstable System 


To demonstrate the design of robust resilient controllers, consider the second-order 
unstable system originally presented in [29] to illustrate the lack of guaranteed gain 
margins for LQG controllers. Specifically, the state space system is given by 


x{t) 

y(t) 


1 

0 


[1 


} i<t) + 

0 ] x(t). 



t e [0, 00) 
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The matrices D x , D 2 , E x , and E 2 are chosen to be 


Dx 


v/60 0 

v/60 0 


D 2 = [ 0 1 ] , 


Ex = 


x/60 \/60 

0 0 


Eo = 


0 

1 


Here we consider uncertainty in the (2, 1) component of the dynamics matrix. Using 
the uncertainty structure given by (6.18), the actual dynamics are given by A+B 0 fC 0 , 
where Bo = [1 0] T and C 0 = [0 1]. 

The quasi-Newton optimization algorithm discussed in Section 6.10 was used to 
compute full-order controllers (n c = 2) that minimize the cost bound J for several 
values of 7 and 7 C for both the multiplicative and additive uncertainty characteriza- 
tions. The actual 1-L 2 cost was computed for a range of values of the controller error 
parameter S and the plant uncertainty /. The cost dependence for the multiplicative 
(relative) uncertainty characterization (6.19) is shown in Figure 6.4. As j c decreases, 
the n 2 cost of the nominal closed-loop system increases while the % 2 cost of the 
perturbed closed-loop system remains near the nominal value for a larger range of 
perturbations. The LQG controller stabilizes the closed-loop system for only small 
perturbations in the controller error parameter, while the resilient controllers stabilize 
the closed-loop system and provide performance close to the optimal level for much 
larger perturbations in the controller error parameter. Hence, robust performance 
over a large range of the uncertain parameter is achieved for some increase in the / H 2 
cost above the optimal. 

The effects of both plant uncertainty and controller uncertainty can be seen in 
the parameter plot shown in Figure 6.5. In this case, the stability of the closed-loop 
system was checked over a grid of the uncertain parameters for both the LQG and a 
robust resilient controller. The dashed line shows the region of asymptotic stability of 
the LQG controller, the solid line corresponds to the robust resilient controller, and 
‘x’ marks the point corresponding to the nominal conditions. It is evident from the 


106 



figure that the robust resilient controller renders the closed-loop system more robust 
to perturbations in both the plant and controller. 

Next, Figure 6.6 shows the Nyquist plots of the loop gain of the system for the 
LQG controller and three different robust resilient controller designs. These plots 
clearly demonstrate the resiliency of the non-fragile controllers over the LQG con- 
troller in terms of their increased gain and phase margins. Figures 6. 7-6. 9 give the 
same plots for the case of additive (absolute) controller uncertainty. Once again, the 
achieved robustness of the robust resilient controllers over the LQG controller is ob- 
vious. Figure 6.7 also shows the cost dependence of an 'H 0O robust controller which 
shows that although the Hoo controller is robust against plant variations, it is highly 
fragile with respect to controller gain variations. The Nyquist plots of the loop gain of 
the Hoo robust controller and the robust resilient controller corresponding to identi- 
cal plant uncertainty levels are shown in Figure 6.10, which clearly shows the relative 
stability superiority of the robust resilient controller over the robust Hoo controller. 

6.12. Two-Mass Benchmark Problem 


Consider the two-mass system shown in Figure 6.3 with mj = m 2 = 1 and an 

*2 


u 



xi=y 

1-*- 


k 


m\ 

AAA 

m 2 


w 


Figure 6.3: Two mass oscillator 


uncertain spring stiffness k [114]. A control force acts on mass 1 and the position of 
mass 2 is measured, resulting in a noncolocated control problem. The nominal system 
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dynamics with k nom = 1 and states defined in the figure are given by 



0 

1 

0 

1 

I 

0 


' 0 ' 


0 

1 

0 ' 

x(t) = 

1 0 

t— * 

0 

1 

0 

0 

1 

0 

x(t ) + 

0 

1 

u(t) + 

0 

0 

0 

0 


1 

-1 

0 

0 

1 


1 

0 
1 


1 

0 

1 


y(t) =[0100 x(t ) + [ 0 1 w(t). 


Using the uncertainty structure given by (6.18), the actual dynamics are given by 
A + B 0 fC 0 , where B 0 = [ 0 0 -1 1 ] T and Cq = [ 1 -1 0 0 ]. The matrices 
Ei and E 2 are chosen to be 


Ei 


0 10 0 
0 0 0 0 ’ 


E 2 = 


0 

v / 0lx)l ' 


As in Section 6.11, the quasi-Newton optimization algorithm was used to com- 
pute full-order controllers (n c = 4) that minimize the cost bound J . The cost depen- 
dence for the multiplicative uncertainty characterization is shown in Figure 6.11. The 
Nyquist plots of the LQG controller and three robust resilient controllers are shown 
in Figure 6.12. For the additive uncertainty characterization, the cost dependence is 
compared to the LQG controller and an Boo robust controller in Figure 6.13. The 
asymptotic stability regions of the LQG controller and the robust resilient controller 
are shown in Figure 6.14. Finally, the Nyquist plots of the LQG controller and three 
robust resilient controllers are shown in Figure 6.15, while the Nyquist plots of the 
Boo robust controller and the robust resilient controller corresponding to identical 
plant uncertainty levels are shown in Figure 6.16. In all cases, the robust resilient 
controllers are superior in their ability to tolerate plant and controller uncertainty as 
compared to the LQG and robust Boo controllers. Furthermore, the robust resilient 
controllers possess far superior gain and phase margins. 
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6.13. Conclusion 


V 


V/ 


w 


In this chapter, we extended the robust fixed-structure guaranteed cost controller 
synthesis framework to synthesize robust resilient controllers for controller gain vari- 
ations and system parametric uncertainty. Specifically, the guaranteed cost approach 
of [10] and [13] was used to develop sufficient conditions for robust stability and 
performance via fixed-order dynamic compensation. A quasi-Newton optimization 
algorithm was used to obtain robust controllers for two illustrative examples. 



Figure 6.4: Dependence of the I-L 2 cost on the controller error parameter 
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Figure 6.6: Loop gain Nyquist plots 
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Figure 6.7: Dependence of the 7^2 cost on the controller error parameter 



Figure 6.8: Asymptotic stability regions 
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Figure 6,11: Dependence of the H 2 cost on the controller error parameter 






Figure 6.12: Loop gain Nyquist plots 
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6.13: Dependence of the cost on the controller error parameter 



Controller Error Parameter, 5 


Figure 6.14: Asymptotic stability regions 





Figure 6.15: Loop gain Nyquist plots 



Figure 6.16: Loop gain Nyquist plots 
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CHAPTER 7 


Fixed-Structure Controller Design 
for Axial Flow Compression 
Systems 

7.1. Introduction 

The desire to develop an integrated control system-design methodology for ad- 
vanced propulsion systems has led to significant activity in modeling and control of 
flow compression systems in recent years (see, for example, [7, 79, 83, 90, 98] and the 
references therein). Two of the main design constraints limiting jet engine compres- 
sion system performance are the compressor aerodynamic instabilities of rotating stall 
and surge. Rotating stall is an inherently three-dimensional 1 fluid dynamic instability 
which is characterized by regions of flow that rotate at a fraction of the compressor 
rotor speed while surge is a one-dimensional axisymmetric global compression system 
oscillation which involves axial flow oscillations, and in some cases even axial flow 
reversal, which can damage engine components and cause flameout to occur. 

A fundamental development in compression system modeling for low speed axial 

compressors is the Moore-Greitzer model given in [90]. Specifically, utilizing a one- 

1 When analyzing high hub-to-tip ratio compressors, rotating stall can be approximated as a 
two-dimensional local compression system oscillation. 




mode expansion of the disturbance velocity potential in the compression system and 
assuming a cubic characteristic for the compressor pressure-flow performance map, 
the authors in [90] develop a low-order three-state nonlinear model involving the 
mean flow in the compressor, the pressure rise, and the amplitude of the rotating 
stall. Starting from infinitesimal perturbations in the flow field, the model captures 
the development of rotating stall and surge. In particular, the model predicts the 
experimentally verified pitchfork hard subcritical bifurcation at the onset of rotating 
stall [89]. 

Using the Moore-Greitzer model, a bifurcation-based controller for rotating stall 
and surge was developed by Liaw and Abed [84], In particular, the Liaw and Abed 
static nonlinear controller is given by 

7throt(A) = 7o + kA 2 , (7.1) 


where 7u,rot(-'4) is the control throttle, A is the amplitude of the rotating stall, and 70 
and k are constants. The locally stabilizing bifurcation-based controller given by (7.1) 
essentially changes the bifurcation structure of the controlled system at the onset of 
rotating stall from a hard subcritical bifurcation to a soft supercritical bifurcation to 
soften the transition into rotating stall. However, as noted by Eveker et al. [36], even 
though the Liaw and Abed controller reduces the abrupt transition into rotating stall, 
it is ineffective for surge. Modifying the static nonlinear controller given by (7.1) to 


7throt(A, $) — 7o k\ A ^2*^5 


(7.2) 


where $ is the time rate of change of the mean flow in the compressor 2 and k\ and 

2 Even though a patented differentiation scheme for sensing the time rate of change of the mean 
flow in the compressor is given in [36], the calculation of 4> can be simply obtained by the most 
elementary equations of fluid dynamics. For example, under the assumption of one-dimensional 
flow, the unsteady axial momentum equation as applied to the bulk of the fluid in the inlet duct 
yields $ = —jjjjAP, where p is the fluid density, L is the duct length, U is the rotor wheel speed, 
and A P is the change in static pressure. 
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&2 are constants, Badmus et al. [8] considerably extended the domain of attraction 
of the Liaw and Abed controller. A fundamental shortcoming of the aforementioned 
controllers is the demanding two-dimensional sensing requirements for implementing 
these controllers. Specifically, measuring rotating stall amplitude is quite challenging, 
requiring pressure sensor arrays distributed circumferentially around the compressor 
annulus, along with discrete Fourier transform software for spatial and temporal fil- 
tering for computing the first circumferential spatial harmonic of rotating stall. As 
an alternative to the locally stabilizing nonlinear controllers developed in [8,36,84], 
the authors in [56, 78, 79] develop globally stabilizing controllers for controlling ro- 
tating stall and surge in axial flow compression systems. In particular, Lyapunov- 
based recursive backstepping globally stabilizing static full-state feedback nonlinear 
controllers requiring rotating stall amplitude measurements are constructed in [56], 
while a globally stabilizing static output feedback nonlinear controller is given in [79]. 
Specifically, the Krstic et al. [79] static output feedback controller is given by 

7.hrot($, *) = r + C ‘^ i -~ C2 ' i ’ , (7.3) 

where 4/ is the pressure rise in the compressor, 4> is the circumferencially averaged 
flow in the compressor, 3 and T, Ci, and c 2 are constants. Even though (7.3) provides 
a simplified sensing architecture over (7.1) and (7.2), the controller is static, pos- 
sessing gain at all frequencies. Furthermore, none of the above controllers have any 
disturbance rejection guarantees. 

In this chapter we develop linear , time-invariant, pressure rise feedback reduced- 
order dynamic compensators for the nonlinear Moore-Greitzer axial flow compressor 
model. Specifically, we construct a modified Riccati equation whose solution guaran- 
tees that the nonlinear closed-loop axial flow compression system is locally asymptot- 

3 Mean flow is relatively simple to measure and is usually measured using pitot probes located in 
the bell mouth of an engine. 



ically stable and the closed-loop output system energy is less than the net weighted 
input energy at any time T in the face of Li exogenous disturbances. Using the 
modified Riccati equation, constructive sufficient conditions for fixed-order (i.e., full- 
and reduced-order) pressure rise feedback dynamic compensators guaranteeing local 
asymptotic stability and disturbance rejection are developed. These sufficient con- 
ditions are in the form of a coupled system of algebraic Riccati equations providing 
explicit design equations for characterizing pressure rise feedback dynamic compen- 
sators that account for compression system nonlinearities and exogenous disturbances 
with a guaranteed domain of attraction. Unlike the nonlinear static and relative de- 
gree zero controllers possessing gain at all frequencies discussed above, the proposed 
linear dynamic compensators explicitly account for compressor performance versus 
sensor accuracy, compressor performance versus processor throughput, and compres- 
sor performance versus disturbance rejection. Furthermore, the proposed controller 
is predicated on only pressure rise measurements, providing a considerable simplifi- 
cation in the sensing architecture. 

7.2. Output Feedback Disturbance Rejection Con- 
trol for Axial Flow Compression Systems 

In this section we introduce the output feedback disturbance rejection control 
problem for axial flow compression systems. The goal of the problem is to deter- 
mine a linear, time-invariant, fixed-order dynamic output feedback compensator that 
stabilizes a given parameterization of the nonlinear Moore-Greitzer axial flow com- 
pressor model using only pressure rise measurements while guaranteeing closed-loop 
disturbance rejection and optimality of a quadratic performance criterion involving 
weighted state and control variables. 

To capture post-stall transients in axial flow compression systems, we use the 
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one-mode Galerkin approximation model for the partial differential equation charac- 
terizing the disturbance velocity potential at the compressor inlet proposed by Moore 
and Greitzer [90]. Specifically, we consider the basic compression system shown in 
Figure 7.1, consisting of an inlet duct, a compressor, an outlet duct, a plenum, and a 
control throttle. We assume that the plenum dimensions are large compared to the 
compressor-duct dimensions so that the fluid velocity and acceleration in the plenum 
are negligible. In this case, the pressure in the plenum is spatially uniform. Further- 
more, we assume that the flow is controlled by a throttle at the plenum exit. Finally, 
we assume a low speed compression system with oscillation frequencies much lower 
than the acoustic resonance frequencies so that the flow can be considered incom- 
pressible. However, we do assume that the gas in the plenum is compressible, and 
therefore acts as a gas spring. 



Figure 7.1: Compressor system geometry 


Now, invoking a momentum balance across the compression system, conservation 
of mass in the plenum, and using a Galerkin projection based on a one-mode circum- 
ferential spatial harmonic approximation for the non-axisymmetric flow disturbances 
yields [90] 

A(t) = — / \E r c ($(*) +.4(t)sin(0))sin(0)d0, t € [0, oo), (7.4) 

3tt J q 

i r i 2 * 

$(t) = -V(t) + — {${t) + A(t)sm(9))d6, (7.5) 

Jo 
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(7.6) 


*(f) = - *r(*)l, 

where $ is the circumferentially averaged axial mass flow in the compressor, '£ is 
the total-to-static pressure rise, A is the normalized stall cell amplitude of angular 
variation capturing a measure of nonuniformity in the flow, is the mass flow 
through the throttle, a and j3 are positive constant parameters, and 4> c (-) is a given 
compressor pressure-flow map. The compliance coefficient j3 is a function of the 
compressor rotor speed and plenum size. For large values of /?, a surge limit cycle 
can occur, while rotating stall can occur for any value of /?. Now, assuming that the 
compressor pressure-flow map 'FcO is analytic, the integral terms in (7.4) and (7.5) 
can be expressed in terms of an infinite Taylor expansion about the circumferentially 
averaged flow to give 

°° 1 d 2 *- l tf C (0 I 2*- 1 


2Jt-l 


3?r T kl(k - 1)! df 


<=*«] ' > 

m 


2 k 


€=*(0 


*(*) 


P : 


2 [$(t) - $T (t)]. 


t £ [0, oo), (7.7) 

(7.8) 

(7.9) 


The specific compressor pressure-flow map, 4/c(-)> which was considered in [90] is 


<M$) = *c 0 + 1 + §$ - l 2 &, (7-10) 

where ^c 0 is a constant parameter. In this case, (7.7)-(7.9) become 

A(t) = ^A(t)[l -$ 2 (t) - \A 2 {t)]+ fiiWiit), t6[0,oo), (7.11) 

m = -*(t) + * C (*W) - mt)A 2 (t) + &«*(*), (7.12) 

’F(t) = - $ T (t)} + p 3 w 3 {t), (7.13) 


to which we have added the C 2 external disturbance signals W\ (t), w 2 (t), and w 3 (t), 
t £ [0, oo), with scaling factors j3\, fi 2 , fa £ R. The proposed additive disturbance 
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model can be used to capture combustion noise and turbine speed fluctuations. For 
example, w 3 (t), t € [0, oo), might reflect back-pressure disturbances to the compressor 
from the combustor. 

Next, note that for fixed values of flow through the throttle, $t(£) = $T eq , (7.11)— 
(7.13) have an equilibrium point given by 


(^eq, $eq, ^eq) = (0. 4>T eq , ^c(^eq))- (714) 

Defining the shifted state variables X\ = A, x 2 = $ — $eqi and x 3 = ^ — ’J'eqj so 
that for a given equilibrium point on the axisymmetric branch of the compressor 
characteristic pressure-flow map the system equilibrium is translated to the origin, it 
follows that the translated nonlinear system is given by 


xi{t) = ^(1 - A 2 )an(t) - |[Jx?(t) +x 1 {t)xl(t) + 2Ax!(t)x 2 (t)] 

+3\W\ (t), t € [0, oo), (7-15) 


x 2 (t) = - x 3 (t ) + |(1 - A 2 )x 2 (<) - \x\{t) - |A xl{t) - §Ax?(t) - fz?(t):r 2 (t) 


+02W 2 {t), (7-16) 

i 3 (t) = -p[x 2 (<) ~ “(<)] + & w 3 (t), (7.17) 

where A = <J>T eq and « = $ T — A. Decomposing (7.15)-(7.17) into a linear and a 


nonlinear part, we obtain the state space model 


Zl(t) ’ 


|(1 - A 2 ) 0 

0 ' 


’ an (t) ' 


r-i o i 

x 2 (t) 

— 

0 f(l- A 2 ) 

-1 


x 2 (t) 

+ 

O 

1 

1 

*40 

•H 


o 
1 

0 


_ x 3 {t) _ 


1 

o 

o 
1 



’ 0 ' 


'01 

0 

0 ' 


’ Wi(t) ' 

+ 

0 

u(t) + 

0 

02 

0 


wait) 


1 

~w J 


0 

0 

ft 


_ w 3 {t) _ 


(HsfoM) 

(7.18) 


where 


yo{t) = 


1 0 0 
0 1 0 


Zl(<) 

x 2 (t) 

*3 (t) 




Xj(t) + Ax x (t)x\(t) + 8\xi(t)x 2 (t) 
3Ax 2 (t) + 3x 2 (t)x 2 (t) + 6Ax 2 (t) + 2 x 3 2 (t) ' 
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Now, it can be shown that the linear part of (7.18) is linearly stabilizable for A > 1, 
while for A = 1, corresponding to the maximum pressure rise equilibrium point, 
the linear part of (7.18) is linearly unstabilizable. With the system written in this 
form, we can now state the dynamic output feedback control problem for axial flow 
compression systems. Here, for generality of exposition, we present the formulation 
for an n-dimensional dynamical compression system. 

Dynamic Output Feedback Control for Axial Flow Compression Sys- 
tems. Given the n th -order stabilizable and detectable 4 nonlinear dynamical system 

x(t) = Ax(t) + B 0 <f>(y 0 (t)) + Bu(t) + D lW (t), * € [0, oo), (7.19) 

y Q (t) = C 0 x(t), (7.20) 

with output measurements 

y(t) = Cx(t) + D 2 w{t), (7.21) 

where u(t ) G R m , y(t ) G R l , y 0 (t) e R'°, t G [0,oo), (j) : R'° E m °, and where 
w(t) G R d , t G [0, oo), is an exogenous C 2 signal, each of whose components has norm 
less than one, determine an n* h -order linear, time-invariant dynamic compensator 

x c (t) = A c x c (t) + B c y(t), t G [0, oo), (7.22) 

u(t) = C c x c (t), (7.23) 

that satisfies the following design criteria: 

1. the undisturbed ( w(t ) = 0) closed-loop nonlinear system (7.19)-(7.23) is asymp- 
totically stable; 

4 Here, stabilizability and detectability are defined with respect to the linear part of the dynamical 
system (7.19)-(7.21). 
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2. the disturbed closed-loop system (7.19)-(7.23) from £ 2 disturbances w(-) to 
performance variables 

z(t) = E\x(t) + E 2 u(t), (7-24) 


satisfies the disturbance rejection constraint 

r »T rT 


f z T (s)z(s) ds < 7 ^ f w T (s)w(s) ds + y(rr(0)), T > 0, w(-) £ £ 2 , 

Jo Jo 

(7.25) 

where z(t) € R p , t € [0, oo), 7 d > 0 is a given constant, x(t) = [x T (t) xJ(f)] T , 
and V(-) is a Lyapunov function for the closed-loop system (7.19)— (7.23); and 


(7.26) 


3. the quadratic performance functional 

rOO 

J(A c ,B c ,C c )= z T (t)z(t) dt, 

Jo 

with w(t ) = 0 is minimized. 

For the three-state parameterized Moore-Greitzer model given by (7.18), the sys- 
tem matrices in (7.19) and (7.20) are given by 

A = 


(7.28) 


1(1 - V) 

0 

i 

o 


r-i 

0 


' 1 0 0 ■ 
0 1 0 

1 

o o 

1(1- A 2 ) 

1 

, Bo = 

1 

o o 

i 

4 

0 

II 


0 


’ ft 

0 

0 ' 

0 

, A = 

0 

ft 

0 

1 

~w J 


0 

0 

ft 


B = 

Furthermore, for pressure rise sensor measurements, 

C = [ 0 0 1 ] . 


(7.29) 


7.3. Sufficient Conditions for Closed-Loop Stabil- 
ity and Disturbance Rejection 


In this section we provide a Riccati equation that guarantees asymptotic stability 
of the undisturbed (t u(t) = 0) closed-loop system (7.19)-(7.23) as well as disturbance 



rejection of the disturbed closed-loop system in the face of £2 exogenous disturbances. 
First, however, note that the closed-loop system (7.19)-(7.23) has a state-space rep- 
resentation given by 


i(t) = Ax(t) -1- £? o 0(2/o (t)) + Dw(t), t e [0, 00 ), 

(7.30) 

y 0 (t) = C 0 x(t), 

(7.31) 

z(t ) = Ex(t), 

(7.32) 

Bo = f n B ° I , Co=[Co 0i oxnc ] . 

Un c x mo J 



Furthermore, we assume that the nonlinear part of (7.19), or, equivalently, (7.30), is 
such that 0( •) e where 

§ v = {0 : V R m °: 0(0) = 0, ||0(yo) II2 < 7n ? llyolll , 2/o G £>}, (7.33) 

where V C R l ° is a closed set and 7 n > 0 is given. For the statement of the main 
result of this section, define the notation n = n + n c , R = E T E, and V — DD T and 
set V = R*° . 

Theorem 7.1. Let {A c , B c , C c ) be given and suppose there exists a matrix P6P" 
and scalars e, 7 d , 7„ > 0 satisfying 

0 = A T P + PA + 7d 2 PVP + 7 ; 2 PBoBj P + CjC 0 + eP + R. (7.34) 

Then the function 

V(x) = x r Px, (7.35) 

is a Lyapunov function that guarantees that the undisturbed (w(t) = 0) closed-loop 
system (7.30)-(7.32) is globally asymptotically stable for all 0(-) € d^io- Further- 
more, the solution x(t), t € [0, cx>), of the nonlinear system (7.30)-(7.32) satisfies the 
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disturbance rejection constraint 


f z r (s)z(s) ds < 7^ f w T (s)w(s) ds + V(xo), T > 0, w(-) € £2- (7.36) 

Jo Jo 

Finally, in the case where w(t) = 0, the performance functional (7.26) satisfies the 
bound 

poo 

J(x 0 ,A c ,B c ,C c )= / z{t) T z(t)dt <V(x 0 ). (7-37) 

Jo 

Proof. First note that since P € P”, it follows that the Lyapunov function 
candidate V(x), x € R" \ {0}, given by (7.35), is positive definite. The corresponding 
Lyapunov derivative along the trajectories x(t), t € [0, 00), of the undisturbed (w(t) = 
0) closed-loop system (7.30)— (7.32) is given by 

V(i(«)) = V'(x(t))[Ax(t) + Bottvom 

= i T (()(i T P+Pi)i(() + 2^ T (j/„(i))BjPi(i), t € (0, 00 ), (7.38) 


or, equivalently, using (7.34) 

V(x(t )) = -x T (t) [eP + 7d - *PVP + R ] x (t) - 7 ~ 2 x T (t)PB 0 B^ Px(t) 

^ T WC o T C o x(O+^ T (y o W)4 T ^W+^ T (O^o0(yo(O), t€[ 0,oo). (7.39) 


Now, adding and subtracting 7n0 T (yo(i))<Kyo(£))> t e [0,oo), to and from (7.39) and 
grouping terms yields 


V(x(t)) = ~x T (t)[eP + 7 - ?PVP + R]x (t) + ^r{yo(t))Hyo{t)) ~ Vo (OltoW 

- 7n“ 1 5 0 T Pi(0-7n</»(yo(0)] T [7n" 1 ^o^W-7n<A(yo(0)], * € [0, oo). (7.40) 


Since eP is positive definite, 7 ^PVP+R is nonnegative definite, and 7n^ T( ^ — 2/o" 2/o < 
0, it follows that V(x(t)) < 0, i € [0,oo), and hence the undisturbed (w(t) = 0) 
nonlinear closed-loop system (7.30)-(7.32) is globally asymptotically stable. 



Next, to show that the disturbance rejection constraint (7.36) holds, note that for 
all w £ R d , 

0 < (7d l D T Px - 7d^) T (7j 1 D T Px - 7d w) 

= r y^ 2 x T PVPx + 7 \w T w - 2x t PDw 

= t ^ 2 x T PVPx + x t Rx + 7 2 w T w — z T z — 2x v PDw. (7-41) 


Now, let to(-) £ £2 and let x(t), t £ [0,oo), denote the solution of the nonlinear 
closed-loop system (7.30)-(7.32). Then 

V(x(t)) = x T (t)(A T P + PA)x(t) - 24 > t ( yo(t))pQ Px(t) + 2x T (t)PDw(t), t £ [ 0 ,oo), 

(7.42) 

which, using (7.34) and (7.41), implies 

V(x(t )) < 2x T (t)PDw(t) — jj 2 x T (t)PvPx(t) — x T (t)Rx(t) 

< 'ylw T (t)w(t) — z T (t)z(t), t £ [0, 00 ). (7.43) 


Now, integrating (7.43) over [ 0 , T ] yields 

V(x(T))<[ [ 7 d^ T (s)u;(s) - z 1 (s)z(s)]ds + F(f(0)), T > 0, «;(•)£ £ 2 , (7.44) 

Jo 

which, by noting that V(x(T)) > 0, T > 0, yields (7.36). 

Finally, to show that the performance functional (7.26) satisfies the bound (7.37), 
note that (7.40) implies 

V(x(t)) < x T (t)Rx(t). (7-45) 


Now, integrating (7.45) over [0, 00 ) yields 

poo roc 

/ z T (t)z(t) dt < - V(x(t))dt. 
Jo Jo 


(7.46) 
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Next, since x(t) ->• 0 as t -» oo, where x(t), t € [0,oo), satisfies (7.30) with w(t ) = 0, 
we obtain 


J(x 0 ,A c ,B c ,C c ) < V(x(0)) - lim V(x(t)) 

£—►00 

- ^(x(0)) 

= xJPxq. j-j 

Theorem 7.1 guarantees global asymptotic stability if with T> — R io . 

However, for the three-state axial compressor model given by (7.18), $(■) € is not 
satisfied for V — R l ° . Hence, to obtain a local stability result for the parameterized 
compressor model given by (7.18), we restrict "D to the set Z? c , where T> c is the smallest 
compact set given by 

Pc = {»o € R‘”: ||«%c)||l < 7„- 2 Mli >, (7.47) 

where 0(y o ), Vo € R /o , is the nonlinear part of (7.18). 

Proposition 7.1. For the axial flow compression system given by (7.18), T> c is 
not empty. 


Proof. Defining 

/(yo) = 7n 2 y?yo - 0 T (j/o)<A(yo), (7-48) 

it follows that 


f(Vo) = 7 n 2 ( x i + ^ 2 ) “ x il x i + 4x 2 (2A + x 2 )] 2 - [3x 2 (A + x 2 ) + 2x 2 (3A + x 2 )] 2 . (7.49) 

Now, since /'( 0) = 0, where /'( 0) denotes the Frechet derivative of /(•) at the origin, 
and 


m = 


27„- 2 o 

0 27 „- 2 


> 0, 


(7.50) 
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where /"( 0) denotes the Hessian at the origin, it follows that the origin is a local 
minimum of /(•). Thus, since /(•) is continuous, there exists a neighborhood of the 
origin where f{yo) > 0, yo € R io \{0}, and hence T> c is not empty. □ 

The size and shape of the set T> c for various values of the parameter 7„ with 
A = 1.1 are shown in Figure 7.2. 



Figure 7.2: Regions, V c , satisfying the sector bound (7.47) 


Next, note that J(x 0 , A c , B c , C c ) < xJPx 0 = trPx 0 xJ, which has the same form 
as the standard 7i 2 cost appearing in standard LQG theory. Hence, we replace x 0 Xq 

by 



0 

o b c v 2 bJ J ’ 


(7.51) 


where V\ G K nxn and V 2 € R lxl are arbitrary design weights such that V\ > 0 and 
V 2 > 0, and proceed by determining controller gains that minimize tr PV . Before 
proceeding, however, we shall require for technical reasons that V 2 = D 2 Dj — a 2 V2, 
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where the positive scalar a is a design variable such that a = 0 if and only if D 2 = 0. 
Next, in the spirit of [11], J(P, A c , B Ci C c ) = tr PV can be interpreted as an auxiliary 
cost which leads to the following optimization problem. 

Optimization Problem. Determine controller gains (v4 c , B c , C c ) that minimizes 
J(P, A c , B c , C c ) = tr PV with P G IP” subject to (7.34) and such that ( A c , B c , C c ) is 
minimal. 

It follows from Theorem 7.1 that by deriving necessary conditions for the Opti- 
mization Problem, we obtain sufficient conditions for characterizing dynamic output 
feedback controllers guaranteeing closed-loop system stability and disturbance rejec- 
tion to C 2 exogenous disturbances. 

7.4. Reduced-Order Dynamic Control for Axial 
Flow Compressors 

In this section we present our main theorem characterizing fixed-order disturbance 
rejection controllers for axial flow compression systems. For design flexibility, the 
compensator order, n c , may be less than the plant order, n. For convenience in 
stating this result, define the notation S = (/„ + oc 2 'y'^ 2 QP)~ l for arbitrary matrices 
Q, P £ N 1 , and A ( = A + \el n . Note that since Q, P £ N 1 are nonnegative 
definite and the eigenvalues of QP coincide with the eigenvalues of the nonnegative 
definite matrix Q l/2 PQ l/2 , it follows that QP has nonnegative eigenvalues. Thus the 
eigenvalues of I n -I- a 2 j^ 2 QP are all greater than one, so that S exists. Furthermore, 
define 

T>$ = {x £ R": C 0 x £ V c and x T Px < /?}, (7.52) 

where P £ N” satisfies (7.34) for a given compensator (yl c , B c , C c ) and /5 > 0. 

Theorem 7.2. Let n c < n and let e, 7 „, 7 d, a > 0. Furthermore, using the 



results of Lemma 5.2, suppose there exist matrices P, Q, P, Q G N” satisfying 


0 = Aj P+ PA e + P[j~ 2 B 0 Bj +'y^ 2 Vi]P + Ri+CJCq — PUP +tJP'EPt±, (7.53) 


0 = (A + [ 1 - 2 B 0 B^+^V l }[P+P))Q+Q(A^ 2 B 0 Bj+^ 2 V 1 )[P+P}) T 

+V l -SQt:QS T +rj L SQt:QS T Tl , (7.54) 

0 = (A-SQ^^BoBj+rfVjPFP+PiAt-SQX+i^BoBj+rfVjP) 

+P( 1 ; 2 B 0 Bj+% 2 [V l + a 2 SQi:QS T ])P+PEP-rlPEPT ± , (7.55) 

0 = (^ + [7 n - 2 B 0 B 0 T +7 d -V 1 -S]P)Q+Q(A £ + [ 7n - 2 J BoB 0 T + 7d - 2 V 1 -E]P) T 

+5QEQ5 T -r x 5Q2Q5 T rJ, (7.56) 

rank Q = rank P = rank QP = n c , (7.57) 

and let ( A c , B c , C c ) be given by 

i*. = T[A - SQX + ( 7 - 2 5oP 0 T + 7d 2 ^i ~ X)P]G T , (7.58) 

B c = TSQC t V 2 -\ (7.59) 

C c = -RpB T PG T . (7.60) 


Then 


P = 


P + P -PG T 
-GP GPG T 


(7.61) 


satisfies (7.34) and ( A c ,B c ,C c ) is an extremal of J(P, A c , B c ,C c ). Furthermore, the 
undisturbed (w(t) = 0) closed-loop system (7.30)-(7.32) is globally asymptotically 
stable for all <£(•) € $ r j 0 . In addition, if <p(-) e $p c , then the undisturbed (ie(t) = 0) 
closed-loop system (7.30)-(7.32) is locally asymptotically stable, and Pa defined by 


V a = {x£ R n : CqX € V c and x T Px < vp}, (7.62) 


where vp = max{/3 > 0 : T>p C R"} is a subset of the domain of attraction of the 
closed-loop system. Moreover, the solution x(t), t € [0, oo), of the disturbed closed- 
loop system (7.30)-(7.32) satisfies the disturbance rejection constraint (7.36). Finally, 
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J(P, A c , B c , C c ) = tr[(P + P)V X + PSQY,QS t }. 


(7.63) 


Proof. The proof is constructive in nature. Specifically, we first obtain necessary 
conditions for the Optimization Problem and then show by construction that these 
conditions serve as sufficient conditions for closed-loop stability and disturbance re- 
jection. For G $p c and w(t) = 0, the estimate of the domain of attraction 
X> A is immediate using maximal closed sublevel sets. For details of a similar proof, 
see [11,58]. O 

Since our design methodology yields reduced-order controllers, we are reducing 
control system complexity by assuring the implementation of a simpler controller for 
achieving disturbance rejection than the full-order controller. By “simple”, we are 
referring to a reduction in control system complexity measured by computer memory. 
In the full-order case, n c = n, set G — T = r = so that t±_ = 0. Now the last 
term in each of (7.53)-(7.56) can be deleted and G and T in (7.58)-(7.60) can be 
taken to be the identity. Furthermore, Q plays no role, so (7.56) is superfluous. If, 
alternatively, the reduced-order constraint is retained and the disturbance rejection 
constraint is sufficiently relaxed, i.e., 7d — > oo, then considerable simplification arises 
in (7.53)-(7.56). 

To solve the design equations (7.53)— (7.56), we employed the homotopic continu- 
ation method presented in [58]. Homotopy algorithms operate by first replacing the 
original problem with a simpler problem having a known solution. The desired solu- 
tion is then reached by integrating along the homotopy path that connects the starting 
problem to the original problem. The algorithm was initialized with 7 n = 7d = oo, 
and the LQG gains designed for the linear part of (7.18). For given values of the 
parameters 7„ and 7d, the algorithm was used to find (A c , B c ,C c ). After each iter- 



ation, 7 n and 7 d were decreased and the current values of (A c , B c ,C c ) were used to 
find feasible values for 7 n and 7 d which were then used as the starting point for the 
next iteration. Complete details of the algorithm are provided in [58]. 

7.5. Active Dynamic Control of an Axial Flow 
Compressor 

In this section we use the design equations (7.53)-(7.56) to design a full-order 
(n c =3) disturbance rejection, pressure rise feedback, dynamic controller for the non- 
linear Moore-Greitzer axial flow compressor model given in Section 7.2. Specifically, 
we chose the data parameter values of o = 3.6, /? = 0.356, and ^c 0 = 0.72, and set 
A = 1.1 in the parameterization given by (7.18). Note that with A = 1.1, the lin- 
ear part of (7.18) is linearly stabilizable, with (7.18) providing an equilibrium point 
close to the desired (A = 1.0) maximum pressure rise compressor operating point. 
Furthermore, we set 7 „ = 12.59 and 7 d = 2.3, and chose design weights 

V l =V l = BB t , V 2 = V 2 = 1, E x = [ 10 1 1 ] , E 2 = J . 

Using the initial conditions Aq = 1.0, $o — 1-866, and = 3.22 to capture 
system transients in the compressor, simulations of the Moore-Greitzer model were 
run in both the open-loop (7throt(0 = 7o = 0.6689) configuration and with the de- 
signed controller given by (7.58)-(7.60) in the feedback loop. Figure 7.3 shows the 
phase portrait of pressure rise versus flow in the compresssor for the case where no 
external disturbances are included in the simulation. It is seen from this plot that the 
controlled system converges to an equilibrium point close to the maximum pressure 
rise equilibrium point, whereas the constant throttle opening drives the system to a 
stalled equilibrium point. Figure 7.4 shows the phase portrait of pressure rise versus 
flow for a set of initial conditions chosen with various radii at angular intervals of 
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Figure 7.3: Phase portrait of pressure versus flow 
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7 r/6 about the maximum pressure rise operating point. Figure 7.5 shows the phase 



Figure 7.5: Phase portrait of pressure versus flow with exogenous disturbances 


portrait of pressure rise versus flow in the compressor when the exogenous £2 signal 
given by 


w(t) = 


0.9e~ f/5 sin(f + 2) 
1.2e _t/10 sin(i/2 + l) J ’ 


is included in the simulations. Once again, the closed-loop system converges to a 
point near the maximum pressure rise equilibrium point while the constant throttle 


opening drives the system to a stalled equilibrium point. 

Figures 7.6-7.9 show the time histories of the stall cell amplitude A(t), the com- 
pressor flow <F(f), the pressure rise in the compressor '&( t ), and the control throttle 
opening 7throt(0> with the exogenous £2 disturbance included in the simulations for 
the constant throttle opening, the closed-loop system with the disturbance rejection 
controller (7.58)-(7.60), the Liaw and Abed [83] controller given by (7.1) with k — 1, 
and the Badmus et al. [8] controller given by (7.2) with k\ = 1 and k 2 = 1. As seen in 
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H Cell 



Time 


Figure 7.6: Stall cell amplitude versus time 



Figure 7.7: Compressor flow versus time 




Time 


Figure 7.9: Throttle opening versus time 
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Figure 7.6, the disturbance rejection controller rejects the exogenous disturbance and 
drives the stall cell amplitude to zero, while the constant throttle opening is unable 
to reject the exogenous disturbance, driving the system to a stalled equilibrium. The 
Liaw and Abed [83] controller does stabilize the maximum pressure rise point but 
has poor disturbance rejection properties. Alternatively, even though, for the given 
disturbance, the Badmus et al. [8] controller has comparable disturbance rejection 
properties to the proposed controller, it requires a considerably more complex sens- 
ing architecture in practice. Figure 7.9 shows a comparison of the throttle opening 
for the three controllers considered as well as the constant throttle opening. It should 
be noted that the maximum throttle opening amplitude of the disturbance rejection 
controller is 1.1286, whereas the maximum throttle opening amplitude for the Liaw 
and Abed [83] controller is 1.669, and the Badmus et al. [8] controller is 5.018. The 



Figure 7.10: Integral squared performance versus time 


Badmus et al. [8] controller has a significantly larger throttle opening amplitude than 
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the other two controllers and therefore can violate actuator limitation constraints 
in practice. Finally, Figure 7.10 shows the integral squared performance response 
Jo zT (t) z {t) dt versus time for the three designs. 

7.6. Conclusion 

A linear, fixed-order (i.e. , full- and reduced-order) pressure rise feedback dynamic 
compensation framework for axial flow compression systems was developed. Unlike 
the nonlinear bifurcation-based and backstepping controllers proposed in the litera- 
ture, the proposed dynamic compensator framework explicitly accounts for compres- 
sor performance versus sensor noise, compressor performance versus controller order, 
and compressor performance versus disturbance rejection. Furthermore, the proposed 
pressure rise feedback controllers provide a considerable simplification in the sensing 
architecture required for controlling rotating stall and surge. Finally, we note that 
bifurcation-based controllers discussed in the Introduction are dependent only upon 
measured quantities as opposed to the proposed controller which requires a model for 
the performance characteristic map. However, it is important to recognize that since 
the proposed controller guarantees robust stability for all <j >{ •) 6 4>x> c , an accurate 
representation of the performance characteristic map is not required. 
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CHAPTER 8 


Fixed-Structure Controller 
Synthesis for Real and Complex 
Multiple Block-Structured 
Uncertainty 

8.1. Introduction 

The ability of the structured singular value to account for complex, real, and 
mixed uncertainty provides a powerful framework for robust stability and performance 
problems in both analysis and synthesis (see [30, 37, 97, 105, 123] and the numerous 
references therein). Since exact computation of the structured singular value is, in 
general, an intractable problem, the development of practically implementable bounds 
remains a high priority in robust control research. Recent work in this area includes 
upper and lower bounds for mixed uncertainty [37, 81, 123] as well as LMI-based 
computational techniques [18,41]. 

An alternative approach to developing bounds for the structured singular value is 
to specialize absolute stability criteria for sector-bounded nonlinearities to the case of 
linear uncertainty [53]. This approach, which has been explored in [23,46,48,49,51,52, 
57,61], demonstrates the direct applicability of the classical theory of absolute stability 
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to the modern structured singular value framework. In particular, the rich theory of 
multiplier-based absolute stability criteria due to Lure and Postnikov [1,82,94,103], 
Popov [102], Yakubovitch [119, 120], Zames and Falb [125], and numerous others 
can be seen to have a close and fundamental relationship with recently developed 
structured singular value bounds. 

The objective of this chapter is to use the absolute stability results of [53], which 
unify and extend existing structured singular value bounds for mixed uncertainty, 
to obtain fixed-structure controllers and fixed-order stability multipliers which pro- 
vide robust stability and performance. Using the results of [53], the robust controller 
synthesis technique proposed here permits the treatment of fully populated real un- 
certain blocks which may, in addition, possess internal structure. Such problems 
arise in a variety of applications, such as the study of modal dynamics, in which 
transformation to ‘standard’ diagonal form may introduce additional conservatism, 
computational complexity, as well as destroying the parameter space of the original 
uncertainty characterization [53]. The ability to address real uncertain blocks is based 
on the use of an appropriate class of multipliers whose structure is compatible with 
the real block uncertainty. Hence, tailoring the multipliers to the structure of the 
uncertainty not only leads to the ability to address more general uncertainty charac- 
terizations but can also lead to less conservative controllers than obtained from the 
standard mixed-// synthesis techniques. This more general class of multipliers has no 
counterpart in standard mixed-// theory [37,105,123]. 

8.2. Absolute Stability Criterion with Generalized 
Positive Real Stability Multipliers 

In this section we review the absolute stability criterion for multivariable systems 
with generalized positive real stability multipliers given in [53]. This criterion involves 
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a square nominal (open-loop or feedback) transfer function G(s) in a negative feedback 
interconnection with a complex, square, uncertain matrix A as shown in Figure 8.1. 
Here, we consider the set of block-diagonal matrices with possibly repeated blocks 



Figure 8.1: Standard feedback uncertainty representation 

defined by 

A bs = { A G C pxp : A = block-diag [I^ x ® A^, . . . , I^ r <8> AJI , 

i* + , 0 Ar+1. • ■ • , hr + c 0 Ar +C ] ; A r i <E R Pi xp ', i = 1, • . • , r; 

A£ € €P iXpi ,i = r + l,...,r + c}, 

where the dimension p, of each block and the number of repetitions of each block 

V 

are given such that xpiPi = p , where v = r + c is the number of distinct uneer- 

Z — 1 

tain blocks. Furthermore, define the subset A C Abs consisting of sector-bounded 
matrices 

A = {A G A bs : 2(A - M 1 )*{M 2 - Mi) -1 (A - My) < (A - My) + (A - My)*), 

where My,M 2 € A bs are Hermitian matrices such that M = M 2 — My is positive 
definite. Note that My and M 2 are elements of A. Alternate characterizations of A 
are given in [53]. 

To draw connections with the structured singular value for real and complex block- 
structured uncertainty, we specialize the set A to the case of norm-bounded, internally 
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block-structured uncertainty. Specifically, by letting Mi = — 7 l I and M 2 = 7 */, 
where 7 > 0, it follows that M — 2y~ l I so that M~ l = ^7 1. In this case, A becomes 

A 7 = {A € A bs : 7(A + j-'iyiA + l~ l I) < (A + 7 -1 /) + (A + 7 -1 ^)*}- 

Now, A G A 7 if and only if <r max (A) < 7 _1 . Therefore, A 7 is given by 

A 7 = {A G A bs ■ O/nax (A) ^ 7 }• 

Next we give the multivariable absolute stability criterion for sector-bounded un- 
certain matrices. To state this criterion, we define the sets V and M of Hermitian 
rational scaling matrix functions by 

V = {D : C -4 <C pxp : D(ju) > 0, D(ju) A = A D{ju), w G R, A G A bs }, 

N = {N : C -> C xp : N(juj) = N*(ju), N(ju)A - A w6R, A e A bs }. 

Furthermore, define the set Z of rational multiplier functions by 

Z = {Z : C -> : Z(ju) = D(ju) - juN(ju), D(-) G 2?, iV(-) G A/ - }. 

Note that if Z(-) G Z, £»(•) G 2>, and iV( ) G JV, then Z(ju) = D{ju) - juN{ju) 
if and only if D{ju) = He Z(juj) and N(ju) = ^Sh Z(ju), u ± 0. Hence, since 
D(ju>) >fl,uGRU oo, Z(-) G Z consists of generalized positive real functions [3]. 

Remark 8.1. Although the condition D(ju ) A = A D(ju) in V arises in complex 
and mixed -// analysis [37], the condition N(ju) A = A *N(ju) in Af has no counterpart 
in [37]. As shown in [53], this condition generalizes mixed-// analysis to address 
nondiagonal real matrices which are not considered in standard mixed-// theory. The 
condition N(ju)A = A*N(ju) is an extension of the condition used in [46] for Popov 
controller synthesis with constant real matrix uncertainty. 
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We now state our main robust stability result. 

Theorem 8.1 [53]. Suppose (I + G(s)Mi)~ l G(s) is asymptotically stable. If 
there exists Z(-) G Z such that 

He [Z(s)(M~ l + (I + G(s)M l )~ 1 G{s))] > 0, (8.1) 

for all s = ju, weRUoo, then the negative feedback interconnection of G(s ) and A 
is asymptotically stable for all A G A. 

8.3. Stability Multiplier Structure 

To ensure the commutability of the rational scaling matrix functions D(s) and 
N(s ) with the uncertainty set A, we must structure Z(s) = C m (sl—A m )~ l B m +D m so 
that Z(s) G Z. To assure this, we construct the multiplier Z(s) from the constituent 
multipliers D(s ) and N(s). Furthermore, instead of obtaining a realization for N(s), 
we obtain a realization for sN(s) directly since Z(s) — D(s ) — sN(s). We therefore 
choose multiplier realizations 



where A d G cM d , A n G <S^ n , B a G *Ss d , B n G Sb„, Gd G Sc d , C n G Sc n , and Dd € Sd 6 - 
The sets <S^ d , Sa„, Sb 6 , Sb u , Sc d , Sc n , and So d are chosen to enforce the diagonal 
structure of D(s) and N(s) given by 

D(s) = Cd(sl - Ad)~ l B d + D d = block-diag [D\ (s) <g> I pi , . . . , D v (s) ® I p J , 

N(s) = C n {sl - A d )~ l B n = block-diag [M (s) <8 > / P1 , . . . , N v (s) <8> I Pv ] , 

where 

0<A(s)eC* ,x *, i = 
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V. 


Ni(s ) = N*(s) e C* x ' pi ,N i (s)A i = Ni{s) A*, * = 1, . . . , 

The structure chosen for the rational functions representing D(s ) and N(s) in this 
chapter is similar to the one for the curve- fitting operation in standard //-synthesis [9]. 
In particular, we define S Ad , S An , $B n , $c d , and S Dd as follows: 

S Ad = {A d = block-diag [7^ lPl ® A dl , ■ . . , I^ vPv ® AiJ : A di e R" di x ” d< , 

A di is in controllable canonical companion form with the 1 st 
row = 0 a di2 0 a di4 ••• 0 adi,„ d . j ,sgn a dij = (— l)^/ 2 l +1 t 

i - 1, — , w, j = l,...,n d) }, 

«S>i n = {>ln = block-diag [I^ lPl ® A ni , . . . , I\p vPv ® -^n,,] : A ni 6 K *, 

A ni is in controllable canonical companion form with the 1 st 
row = 0 a ni2 0 a ni4 ••• 0 a n< | , i = 1, . . . , u}, 

S Bd = {B d = block-diag [I^ m <S> B dl , . . . , I^ vPv ® B d , ] : B di 
Bj.= [ 1 0 ••• 0]',» = 1,. 

S Bn = {B n = block-diag [I^ m ®B ni ,..., I^ vPv ® £„„] : B ni G R" n< xl , 

Bl = [ 1 0 ••• 0],i = l, 


€ R" d > xl , 


( 8 . 2 ) 


(8.3) 


(8.4) 


(8.5) 


Cdi = 

0 c d il2 0 c d i>4 0 c d . nd< 

, sgn c dij = (-1 Y /2 , 


* = 1. • 

-6 

£ 

rH 

II 

sT 


(8.6) 

[C a = [block-diag [1^ ® C ni , . . . , /^ Pv ® CnJ] : € R lxn "‘ , 



c n ; ,i 0 C n;3 ••• 0 C ninnji 0 , * = 1, . . . , v}, 

(8.7) 

{D d = block-diag [d dl I ^ lPx , . . . , d dv I^ vPv ] : d, 

3, G R,d di > 0, 


i = 1,- 



(8.8) 
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where n di and n ni are the a priori fixed even orders of the rational functions repre- 
senting the z th diagonal element of the multipliers D(s ) and N(s), respectively. Thus 
we see that n d and n n , defined as 


n d 


— ^ n d,1pii n n — ^ j^ii 


(8.9) 


i=l 


i=l 


are the total number of states describing D(s) and N(s), respectively. 

To construct Z(s ) = D(s) — sN(s), we obtain the augmented realizations 


Am 
B m 

D m 


block-diag [Ad, A n ] , 

B d 
B n 


[c d -c n ], 

D d , 


( 8 . 10 ) 

( 8 . 11 ) 

( 8 . 12 ) 

(8.13) 


where A d G S. 4 d , A n € S An , e 5 Bd , G <S B „, C d G Sc d , C n G <S Cn , and D d G 5 Dd . 
Note that A m G R nmXnm , where n m = n d + n n . Furthermore note that there is no 
contribution to D m from the sN(s) term. This is due to the fact that the rational 
function N(s ) is strictly proper and has only even powers of s. Thus N(s ) necessarily 
has a relative degree of two and hence sN(s) is strictly proper. 

Next, we note that with 

Z(s) 

as defined in (8.10)— (8.13), we obtain the necessary commutability properties, as well 
as ensuring that D(s) > 0 and N(s) = N*(s). However, to satisfy the condition 
N(s ) A = A*N(s), we require that the t th block of N(s) be zero whenever the i th 
uncertainty block A * is complex. Thus, in the case where A* G 0 ><xp< , we require 
that the realization of Ni(s ) be given by 



B m 

c m 

D m 


A Ui — oxo, B n% — [ ]oxi> Dth — [ ]ixo> 
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where [ ] 0X j is the Ox j empty matrix [96]. Thus, in this case, the * th block of N(s) 
is given by 

N i (s) = C ni (sI-A ni )~ l B l ni =0. 

Finally, we note that the stability multipliers cannot have arbitrary realizations 
and still be elements of their appropriate sets Sa 6 , S a „, S B(j , S Bn , S Cd , S Ca , and SD d - 
Thus we recast the stability multiplier matrices so that the free parameters appear 
along the diagonal of a separate matrix, /C m . The stability multipliers can then be 
constructed as 


•^m — -TmC + A mL^m ^mR) B m — flmCi 

Cm — Cyn^KlmCmRi DmRj 

where the matrices A mC , ^ m L, -4mR, Bmc , CmL, C mR , D mL , and D mR are structured 
appropriately. To illustrate the structure of these matrices, consider the scalar mul- 
tiplier 


Z(s ) = dd l + 


1,2 




a d\ % 2 


*71 1,2 


The gain and structure matrices for the stability multiplier with n n = n d = 2 are 
then given by 

0 
0 
0 
0 

J 


4mC — 






r a d \ t 2 

0 

0 

0 





0 ^^1,2 

0 

0 




fCm 

= 

0 

0 

ddi 

0 





0 

0 

0 

a n 1,2 





0 

0 

0 

0 

0 

0 

0 

o' 


'1 

0 0 

0 0 

1 

0 

0 

0 

A 

0 

0 0 

0 0 

0 

0 

0 

0 

, ^mL — 

0 

0 0 

1 0 

0 

0 

1 

0 


0 

0 0 

0 0 



‘0 

1 

0 

o' 


0 

0 

0 

0 

, AnR = 

0 

0 

0 

0 


0 

0 

0 

1 


0 

0 

0 

0 
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B m c ~ 




V 


w 








T 

0 

i 


'0 

0 

0 

o' 


'O' 


0 

1 

0 

0 


0 

? CmR — 

0 

0 

0 

0 

, AnR — 

1 

1 

_o_ 


0 

0 

0 

0 


0 


0 

0 

1 

0 


0 


C'mL =[0100 -1] T D mL = [0 0 1 0 0], 


With these definitions, we see that 


1 

o 

a d 1,2 

0 

1 

o 


1 

1— » l 

1 

1 

0 

0 

0 


0 

o 

0 

0 

a n 1,2 

) — 

1 

1 

o 

0 

1 

_ — 1 

o 


1 

o 
1 


Cm [ 0 ^<fi,2 i,i 0 ] » [ & di ] > 


and thus 


Z(s) = C„(«/ - An) _1 An + D m = d dl + 2 " M ■ 

s 2 - a dl 2 s 1 - a nx , 

8.4. Decentralized Static Output Feedback Formu- 
lation 


In this section we review the decentralized static output feedback problem for- 
mulation for fixed-structure controller synthesis [14,34]. Consider the (m + q + 1)- 
vector-input, (m + q 4- l)-vector-output decentralized system shown in Figure 8.2, 
where ej and di, i = 1 are used to account for model uncertainty, w is the 

exogenous disturbance input, z is the performance variable, and the signals yi and 
Ui, i = 1 ,...,m, are measurement and control signals, respectively. The decentral- 
ized static output feedback multi vector-input, multi vector-output system shown in 
Figure 8.2 is characterized by the dynamics 

m q 

x(t) = Ax(t) +£*,<•>(«)+£ a. d k (t)+B w w(t), £ € [0, oo), (8.14) 

j - l Jfc=i 

m q 

Vi(t) = c yi x{t) +T 'Dyu,} Uj (t) +T v ydik d k (t) + Vy Wi w(t), i = 1, 2, . . . , m, (8.15) 

j = i *=i 
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Figure 8.2: Decentralized static output feedback framework 

rn Q 

6i{t ) — C ei x(t) +^V eUtJ iij(t) 'D e d tk dk{t)+'D ew ,w(t), i = 1,2, . . . ,q, (8.16) 

j - l k = 1 

m q 

z(t) = (8-17) 

3 = 1 k = 1 

In the above formulation, model uncertainty is represented by the decentralized 
static output feedback map 

di(t) = Aje<(£), i=l,...,q, (8.18) 

where the uncertain matrices A* are not necessarily distinct. To represent decentral- 

ized static output feedback control with possibly repeated gains, we consider 

Ui(t) = JCiyi(t), i = (8.19) 

where the matrices tCi are not necessarily distinct. Reordering the variables in (8.18) 
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and (8.19) if necessary and defining 





Ui(t) 


»i(0 


'diW' 


ei (0 

u{t) = 


. y(t) = 

\ 

, d(t) = 

| 

, e(i) = 



u m {t) 


,!/raW_ 






5 - 

[B„, 

••• b„,] 

<311 




_^y m _ 



'e\ 


J 


V = 

5 ^eu — 


S. = [B* 

••• 

= [!>«., •• 

* ^2tl m ] 

T> 

^2/un 

V 

A 

^SMll 

'Dydiq 1 

: 

\ 

,v yd = 




T) ... 

^ V u m 1 

V 

^ yUmm 


'Dydmi 

’ 'Dydmq \ 

^etxn 

T> 

*^euim 

, V ed = 

T^edn 

... n> 

Cl 

•Q 

1 

} 

T) 

Vzu qm m 


J^ed q i 

'^edqq 



1 

ft 

P — 

’T) 

^ ew\ 

l^ywm _ 


V 

^ ew q 


T^zd — ‘ ' * 'Dzdq ] > T-^yw — 

(8.18) and (8.19) can be rewritten 

d(t) — Ae(t), 
u(t) = K,y(t), 

where A and /C have the form 

A = block-diag [1^ <g> Ai , . . . , I$ v ® A„] , 
K = block-diag ® K ,\ , . . . , ® /C 9 ] , 


( 8 . 20 ) 


( 8 . 22 ) 


(8.23) 


(8.24) 


(8.25) 

(8.26) 


(8.27) 

(8.28) 


where v is the number of distinct uncertainties A i 6 C p * xp * or R Pi xp ‘, tpi is the number 
of repetitions of uncertainty A*, g is the number of distinct gains fCi G R r ‘ xc ’ and 
4>i is the number of repetitions of gain K-i . Note that ICi, . . . , JC g are not necessarily 

V g 

square matrices, and ipi — Q and <f>i = m. 


i= 1 


«=i 
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With the definitions in (8.20)— (8.24), the transfer function G(s ) from [u T , d T , u; T ] 
to [y T , e T , z T ] T of the decentralized system has the realization 


G(s) 


\A_ 

i 33 

s 
1 

Cy 

Dyu t ^ yd. , T^yw 

C e 

T^eu f , T^exv 

_c z 

zu , ^ zd , zw , 


(8.29) 


which represents the linear, time-invariant dynamic system 

£(t) = Ax(t) + B u ii(t) + B d d(t) + B w w(t), t e [0, oo), (8.30) 

y(t) = C y x{t) + V yu u(t) + V yd d{t) + Vy W w(t), (8.31) 

e(t) = C e x(t) + / D eu u(t) + X? e< /d(t) + T> ew w(t ), (8.32) 

z(t ) = C z x(t) + V zu u(t) + V zd d{t) + V zw w(t), (8.33) 

which is equivalent to (8.14)-(8.17). Furthermore, by rewriting the decentralized 
control signals (8.19) in the compact form given by (8.26), the closed-loop system 
realization from [d T ,ui T ] T to [e T , z T ] T is given by 



' A 

Bo 

JLl 

G{s) ~ 

A 

Do A 1 - 


. E 

Ex 

E 0 J 


(8.34) 


where 


A = A + B u lCL^C y , Bo = B d + B u ICL^V yd , D = B w + B u JCL^V yw , 

C 0 = C e + V eu K,Ll l C y , Do = V ed + V eu )CL- K l V yd , D x = V ew + V eu JCL^V yw , 

E = C Z + V zu KLz l C y , E x = V zd + V zu ICL^V yd , E 0 = V zw + V zu ICL^V yw , 

and where L K = I - V yu K. Note that we assume det(Ljc) 7 ^ 0 for all K. given by 
(8.28) to ensure the well-posedness of the feedback interconnection. 

Finally, given the closed-loop system realization given by (8.34) with Do = 0, 
Dx =0,Ex = 0, and E 0 = 0, and the multiplier realization given by 

Z(s) 



B m 

C m 

Dm 
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V 






V 








the realization of G(s ) = Z(s)(M _1 + (I + G e d{s)Mi)~ 1 G e d(s)), where G e d(s) is the 
closed-loop transfer function from uncertainty inputs d(t) to uncertainty signals e(t), 
is given by 

G(s) 


where 


‘ A 

Bo 1 

Co 

Do J 




A — Bo Mi C 0 0 

BmCo A m 


B 0 = 


Bo 

B m M _1 


Co = [ D m C 0 C m ] , D 0 =[ D m M~' ] . 


8.5. Specialization to Centralized Strictly Proper 
Dynamic Compensation 


In order to give a more concrete illustration of the decentralized static output 
feedback framework, in this section we provide an example with a single uncertain 
block using a centralized, strictly proper dynamic controller. Specifically, consider 
the uncertain dynamical system 


x(t) = Ax(t) + Bu(t ) + B 0 d(t) + Diw(t), t £ [0, oo), (8.35) 

y(t ) = Cx(t) + Du(t ) + F\d(t) + D 2 w(t), (8.36) 

e(t) = C 0 x(t ) + F 2 u(t), (8.37) 

z(t) = E\x{t) -|- E 2 u(t), (8.38) 

with uncertain plant perturbations AA = B 0 ACo, A B = BoAF 2 , AC = BiACo, 
A D = F x AF 2 , of the nominal system matrices (A,B,C,D). 

The dynamics of the centralized, strictly proper controller are given by 


x c (t) = A c x c (t) + B c y(t), t £ [0, oo), 
u(t) = C c x c (t), 
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(8.39) 

(8.40) 



so that the closed-loop system can be written as 


x(t) = j4£(t) + Bod(t) + Dw(t), t € [0, 00 ), 

(8.41) 

e(t) = C 0 x(t), 

(8.42) 

■to 

11 

'’■to' 

N 

(8.43) 


where 

A 

B C D 2 

Co = [Co F 2 C c ], E= [A AC c ]. 


A 

'x(t) ' 

A- 

r a 

BC c " 


x c (t) 

1 si — 

B c C 

A c + B c DC c 


, A = 


A 




Writing this system in the decentralized static output feedback framework, we 
obtain 


and 


3 

x(t) = Ax(t) + Uj(t) + Bdd(t ) + B w w(t), t € [0, oo), 

i - i 


3 

Ui(t) — By { x(t) + ^ ] B yUij Uj(t ) + 2>*d(f) + By Wi w(t), i = 1,2,3, 

j'=i 

3 

e(<) = C e x(t) + ^^T> eUj Uj(t), 

3= i 
3 

z(t) = C z x(t) + T>zujUj{t), 

3 = 1 


(8.44) 

(8.45) 

(8.46) 

(8.47) 


Ui(f) — A c yi{tj, ^2 (0 — B c y 2 {t')i ^(0 


where 


11 > 

'/I o' 

0 oj ’ 

f 

0 

0 

l- 

<11 
^ s 

EQ 

1 

u 

0 ^ 

1 

<11 

tM 

£Q 

n ^ \B] » * B 0 

*-*u 3 q 1 q y 

B t D ' 
w [ 0 

C — 

jo 

Vyun — 0) 

Byu\2 = 

B yu 13 = 0 , 'Bydi ~ 0 , 

= 0 

c — 

[C 0 ], 

V =0 

Byu 22 = 0 ) 

A“23 = A Byd 2 — Fl, 

= D- 

c = 

[0 /„], 

By U3l = 0 , 

32 = 

A / u 33 = Byd 3 = 0 ? 

= ® 

C e = [Co 0 l, 

V eul = 0, 

V eu2 = 0 , 

Veu 3 = F 2 , 


O 

,^L 

<11 

** 

A „1 = 0 , 

A «2 = 0 , 

A«3 = A- 
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Next, defining 



Ui{t)‘ 

A 



U 2 (t) 


V2(t) 


U 3 {t)_ 


mi 1 ). 


B u =[B ul B u2 B u3 ] , V eu = [V eul V eu2 Veu 3 ],Vzu=[Vzul V zu2 V tu3 ], 


A 

r r i 

A 

\v 

Vyu\\ 

V 

Vyuyi 

V 1 

A 

'Bydl 

A 

1 

*^yw i 

c = 

L 'y 

c 

V — 

V 

^yu 21 

V 

V 

^yu 23 


'Dyt'i 

T) — 
? t^yw — 

V 

^yw2 


S'y 3. 


V 

L^y«3i 

V 

32 

^yu33 J 




V 

L^F3j 


and rewriting the decentralized control signals in the compact form 

u(t) = JCy(t), 


w 




where 



A c 0 0 

0 B c 0 , 

0 0 C c 


the system matrices A, B 0 , D, C 0 , and E in the closed-loop dynamics 


x(t) — Ax(t) + B 0 d(t ) + Dw(t), t £ [0, oo), 
e(t) = C 0 x(t ), 
z(t) = Ex(t), 


can now be written as 


A = 
B 0 = 

D = 
Co = 
E = 


A + B u ICL^Cy = 

B w + B u K.LfAV yw = 

B w + BuKLfc'Dyw = 

C z + T> zu K,L K l C y = [ Co F 2 C c ] , 
C z + V^ICL/CCy = [ Ei E 2 C c ] . 


Bo 

B C F X J ’ 

Di 

B c D 2 ’ 


A BC c 
B c C A c + B c DC c 


1 


(8.48) 

(8.49) 

(8.50) 


Furthermore, closing the uncertainty loop from d(t ) to e(t) yields the closed-loop 
realization 


Ga(s) ~ 


,4 + SqACo 

D ' 

E 

0 . 


v 
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8.6. Robust Stability and Performance 


In this section we present sufficient conditions for robust stability of the feedback 
system given in Figure 8.1 and a robust U 2 performance bound for ||G* w (s )||2 ■ 
Specifically, the following theorem provides a sufficient condition for ensuring that 
Z(s)[M~ 1 + [I + G ed (s)Mi]~ 1 G ed (s)] is strongly positive real [48], as well as providing 
a robust performance bound for ||(5 2U ,(s)||2 • For the statement of the next theorem 
define the notation 

E= [ E 0 ] , D = 
and recall the definition for a strongly positive real transfer function [48]. 

Theorem 8.2. Suppose there exists P G N” +Tlm satisfying 

0 = A T P + PA + (B^P- Cq) (Do + £>J)- l (BjP - C 0 ) T + E T E. (8.51) 

Then Z(s)[M~ l + [I + G ed (s)Mi}' l G ed (s)] is strongly positive real. Consequently, 
the feedback interconnection of G(s) and A is asymptotically stable for all A € A. 
Furthermore, 

\\G zw (s)\\l<tvPDD T . (8.52) 

To apply Theorem 8.2 to controller synthesis, we use the modified Riccati equation 
(8.51) to guarantee that the closed-loop system is robustly stable. This leads to the 
following optimization problem. 

Optimization Problem. Determine gain matrices /C €E and 

JC m e R( nm +«) x ( n m+v) minimize 

J()C,IC m )=tTpDD T , (8.53) 

where Pelf satisfies (8.51), and where h = h + n m . 


D 

0 
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8.7. Sufficient Conditions for Fixed- Order Robust 
Compensation with Dynamic Multipliers 


In this section we state sufficient conditions for characterizing dynamic output 
feedback controllers and dynamic stability multipliers guaranteeing robust stability 
and robust performance. For the statement of the next theorem, partition the 
matrices P and Q as 


P = 



Pu 

P 22 


A _ Qll Q 12 

Q ~ [ QL Q22 


where P U ,Q 11 € N” and P 22 M 22 6 N 71 "’, and define 


Ori^] xri 


0 Cl * xci 

0r 2 </>2xri 


0c 2 lfaXCi 

Ori-i^-ixn 


Oci_l<Ai-lXCi 

0 

X 

1 

. Q*, = 

0 Cl (j-])xc, 




0 r;+1 <j!> i+l xr t 


Qci + l<t>i+lXCi 

^ r v <p v xri 


_ ^Cu^ifXCj . 


(8.54) 


where and C{ are the dimensions of the i th controller gain, ICj € R ri XCi , i = 1, . . . , v, 


and j = 


Theorem 8.3. Suppose there exists n x h nonnegative definite matrices P and 
Q satisfying 

0 =A T P + PA + {BjP - C 0 )(£)o + D^)- l {BjP - C 0 ) T + E T E, (8.55) 

0= (A-Bo{Do+D^)-PB^P-Co))q+q(A-B 0 (Do+DJ)- 1 {BJP-Co)) \dD t , (8.56) 


and let K-i satisfy 


4>i 


0 = 2 { I + T,T yuL~K T EA) [B^PQr+PnQlAC^-CjM^+B^PB^ 
j = 1 


,T 

yw 
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+Vj u C z QiCy + ^Ju[Bl(PnQi + P 2 QI 2 ) - M l Bj{P l Q 1 + PnQ] 2 )]tf 
+B%[(PiQi + PnQu){PiPd + PnB m M~ l - CjD m ) 

+(PiQi2 + Pi2Q2)(PlB d + P 2 B m M~ x - Cl)](D m M~ l + M~ x D m )~ x T>y d 
-Vj u D m (D m M~ x + M~ x D m )~ x [(P\B d + P l2 B m M~ x - CjD m ) T Qx 

+(P? 2 B d + P 2 B m M~ x - ClfQj 2 ]Cj}L-jQl, (8.57) 

Furthermore, let /C m satisfy 


0 — diag{2i4^ L (Pj^(5i2 + 

- 2C^ L (D m M -1 + M~ l D m )~ l (Tj Q\ 2 + T 2 Q 2 )CmR 
-2 Dl L (D m M~' + M~ l D m )~ x {T^Qi + TjQj 2 )CjDl R 
-Dl h (D m M- x + M~ l D m )~ x QiT x + TjQ{TxM~ x 

-\-M X T 2 Q 2 T 2 + TjQ 2 T 2 M 1 + 2M X T? Q\ 2 T 2 4- 2T^Q l2 TiM ) 
•(D m M~ x + M~ l D m )~ l D^r}, (8-58) 


where 

Tj = P l B d + P\ 2 B m M~ l - Cj D m , T 2 = PlB d + P 2 B m M~ x - C £ 

Then Z(s)[M~ x + [I + G ed {s)M x )~ l G ed {s )] is strongly positive real. Thus the closed- 
loop system from w(t) to z(t) is asymptotically stable for all A e A. Furthermore, 
the worst-case % 2 performance of the closed-loop system satisfies the bound 

\\G zw (s)\\l<trPDD T . (8.59) 

Proof. First we obtain necessary conditions for the Optimization Problem and 
then show, by construction, that these conditions serve as sufficient conditions for 
closed-loop stability and robust 'h 2 performance. Thus, to optimize (8.53) subject to 
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V 


V. 








(8.51), form the Lagrangian 


£(/C,/C m ,Q, A) = tr [APM T + <0(i T P + Pi 

+(P 0 T P - C 0 )(Do + Dor'i^P - C 0 f + E T P)] , 


where the Lagrange multipliers A > 0 and Q £ R nxn are not both zero. By viewing 
/C, /C m , and P as independent variables, we obtain 


H = (a - Bo (L>0 + A t ) - 1 (Ap - A)) Q 

DP v 7 

+Q (i - B 0 (A> + D^)-\BJP - C 0 )) T + A DD t . (8.60) 

If ^4 — Bq(Dq+ D j)~ l (B q P — Co)) is Hurwitz, then A = 0 implies Q — 0. Hence, it 
can be assumed without loss of generality that A = 1. Furthermore, note that Q is 
nonnegative definite. Thus the stationary conditions with A = 1 are given by 


dC 

dP 


dC 

dICi 


dC 
dIC m 


= (A- Bo(D 0 + Dq)~ 1 (Bq P - Co)) Q 

+Q (A - Bo(Do + &S)- 1 {bJp - c o) y + DD t = 0 , 

= 2 E Ql» ( 7 + V v» A T ^ T ) [fif {PiQi + P12QDK - CjMlVyd) 

3 = 1 

+1 B] ZP,B w V* yw + VJAQK + Vl\Bl{PlQi + P 2 QJ 2 ) 
-M 1 Bj(P 1 Q l WnQl)K+Bl[(PxQi+Pi2Qj 2 )(PiB^P l 2B m M- l -CjD m ) 


+(PiQi2 + PnQ 2 ){PlB d + P 2 B m M~ l - Cl)\(D m M~ x + M~ x D m )~ x V^ d 
-Vj u D m (D m M~ l + M~ l D m )- l [(PiB d + PnB m M- x - CjD^Q, 
+(PnBd + P 2 B m M~ l - ClfQlWWQlj = 0, 


= diag{2A£ L (Pj^Qi2 + PzQzM^r 

-2Cl h {D m M~ l + M- l D m )- l {T7Qi2 +TjQ 2 )Cl R 
-2 Dl h (D m M-' + M- x D m )- x (T?Q l+ TjQj 2 )CjDl K 
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- Dl L (D m M + M~'D m )~ l (M-'TjQ.T, + T?Q,T,M~' 

+M-'TjQ 2 T 2 + TJQ 2 T 2 M - 1 + 2 M~'TjQ n T 2 + 2T?Qj 2 T,M~') 
■(D m M-' + M-'D m )-'Dl R )= 0, 

which are equivalent to (8.56)-(8.58). Equation (8.55) is a restatement of (8.51). It 
now follows from Theorem 8.2 that if Z(s)[M" 1 + [I + G' e d(s)M 1 ]- 1 G' e( i(s)] is posi- 
tive real then A + B 0 AC 0 is asymptotically stable for all A € A. Finally, the K 2 
performance bound (8.59) is a restatement of (8.53). □ 

Equations (8.55)-(8.58) provide constructive sufficient conditions that yield dy- 
namic controllers for robust fixed-order (i.e. , full- and reduced-order) output feedback 
compensation. By using these equations within a numerical optimization algorithm, 
the optimal robust fixed-order controllers and stability multipliers can be determined 
simultaneously, thus avoiding D, N — K iterations. 

8.8. Quasi-Newton Optimization Algorithm 

A general-purpose BFGS quasi-Newton algorithm [26] can be used to calculate the 
controller gains and the stability matrices, as described in Section 2.3. One require- 
ment of gradient-based optimization algorithms is an initial stabilizing design. For 
full-order controller design, the algorithm can be initialized with an LQG controller, 
while for reduced-order control, the algorithm can be initialized with a balanced trun- 
cation of an LQG controller. Small values should be chosen for Mi and M 2 so that 
the design equations (8.55) and (8.56) can be solved. The quasi-Newton optimiza- 
tion algorithm can then be used to find the controller gains A c , B c , and C c and 
the multiplier gains A m , B m , C m and D m . After each iteration, M\ and M 2 are in- 
creased and the current values of the controller gains (A c , B c ,C c ) and multiplier gains 
(A m , B m , C m , D m ) are then used as the starting point for the next iteration. 
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8.9. Conclusion 


This chapter used the absolute stability results of [53] to obtain fixed-structure 
controllers and fixed-order stability multipliers which provide robust stability and 
performance. By satisfying certain commutability properties with the uncertainty 
set, the stability multipliers designed by the robust controller synthesis technique 
proposed here permit the treatment of fully populated real and complex uncertain 
blocks which may, in addition, possess internal structure. Hence, tailoring the mul- 
tipliers to the structure of the uncertainty not only leads to the ability to address 
more general uncertainty characterizations but can also lead to less conservative con- 
trollers than obtained from the standard mixed-// synthesis techniques. Furthermore, 
since the numerical optimization routine optimizes over the free parameters in the 
controller and the stability multipliers simultaneously , this methodology avoids the 
D, N — K iterations of standard mixed-// synthesis techniques. 





CHAPTER 9 


Concluding Remarks and 
Recommendations for Future 
Research 

9.1. Conclusions 

In this report, we introduced a decentralized static output feedback framework 
for fixed-structure dynamic controller synthesis. As a special case of this framework, 
we showed how a centralized dynamic output feedback control problem can be trans- 
formed to a decentralized static output feedback form. By using this format, a nu- 
merical optimization scheme can be used to optimize the controller gains with respect 
to a given cost function and constraint equation. Furthermore, we demonstrated its 
effectiveness on the ACTEX flight experiment. 

Next, we used the decentralized static output feedback framework to synthesize 
stable ^-optimal controllers by including the W2 cost of the controller in the La- 
grangian and using a multiobjective optimization technique. It was numerically shown 
that for some systems, namely minimum phase, open-loop unstable or non-minimum 
phase, open-loop stable plants, a stable controller can rival the performance of an 
unstable %2-optinial LQG controller and yet not be constrained by the loop mar- 
gins of unstable controllers. For other systems, however, there could be a significant 



degradation in performance by requiring the controller to remain stable, although 
this technique provided controllers yielding the minimal %2 closed-loop cost for all 
stable linear controllers. 

By exploiting the ability of the decentralized static output feedback framework 
to a priori fix the structure of the controller, we were able to synthesize % 2 -optimal 
relative degree two controllers. This was accomplished by cascading two strictly 
proper dynamic controllers in the feedback loop and optimizing over the free controller 
parameters. It was shown that constraining the controller to have a relative degree of 
at least two only marginally increased the H 2 cost of the closed-loop system, though it 
was noted that changing the order of the separate cascaded controllers in the feedback 
loop does significantly affect the 7f 2 cost of the closed-loop system and the natural 
frequencies and break frequencies of the controller dynamics. 

Next we focused on robust control by extending the implicit small gain guaranteed 
cost bound [54] to controller synthesis. Specifically, the implicit small gain guaranteed 
cost bound was used to address the problem of robust stability and 7f 2 performance 
via fixed-order dynamic compensation, and a quasi-Newton optimization algorithm 
was used to solve the coupled nonlinear design equations. 

We then extended the robust fixed-structure guaranteed cost controller synthesis 
framework to synthesize robust resilient controllers, or controllers which are robust 
in the face of system parametric uncertainty and variations in the controller gains 
themselves. Specifically, the guaranteed cost approach of [10] and [13] was used to 
develop sufficient conditions for robust stability and H 2 performance via fixed-order 
dynamic compensation. 

The robust fixed-structure control framework was then used to develop linear, 
fixed-order (i.e., full- and reduced-order) pressure rise feedback dynamic compen- 
sators for axial flow compression systems. Unlike the nonlinear bifurcation-based and 
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backstepping controllers proposed in the literature, the proposed dynamic compen- 
sator framework explicitly accounts for compressor performance versus sensor noise, 
compressor performance versus controller order, and compressor performance versus 
disturbance rejection. Furthermore, the proposed pressure rise feedback controllers 
provide a considerable simplification in the sensing architecture required for control- 
ling rotating stall and surge. 

Finally, we used the absolute stability results of [53] to obtain fixed-structure 
controllers and fixed-order stability multipliers which provide robust stability and 
performance. By tailoring the multipliers to the structure of the uncertainty, we 
could permit the treatment of fully populated real and complex uncertain blocks 
which may, in addition, possess internal structure. This not only leads to the ability 
to address more general uncertainty characterizations but can also lead to less con- 
servative controllers than obtained from the standard mixed-// synthesis techniques. 
Furthermore, by using a numerical optimization routine, we could optimize over the 
free parameters in the controller and the stability multipliers simultaneously, thus 
avoiding the D, N — K iterations of standard mixed-// synthesis techniques. 

9.2. Recommendations for Future Research 

The decentralized static output feedback formulation and the quasi-Newton op- 
timization algorithm discussed in Chapter 2 are programmed into a MATLAB tool- 
box [35]. However, in its present form, it is not the most user-friendly program, and 
making it easier to use would most likely spread its use. Also, as it now stands, the 
quasi-Newton optimization algorithm has only one search method. Adding routines 
such as “double dog-leg” searches or locally constrained optimal “hook” steps [26] to 
the Armijo-type routine already programmed would enhance the robustness of the 



optimization routine. Furthermore, adding numerical gradient computations instead 
of relying solely on analytical gradient computations will extend the range of the 
problems the routine will be able to solve. 

In Chapter 4, we designed controllers with a minimum guaranteed (vector) rel- 
ative degree of two. However, the theory itself is very general, and we can use the 
same formulation to design controllers of (vector) relative degree r. Furthermore, the 
methodology used to obtain relative degree two controllers was solely based on the 
structure of the controllers, and was not constrained to a particular notion of opti- 
mality. In Chapter 4, we designed H 2 -optimal controllers, but we could just as easily 
design T^-optimal (vector) relative degree two or even relative degree two robust 
controllers. 

In Chapter 5, robustness was guaranteed by bounding the effects of the uncertain 
terms in the Lyapunov equation. However, this was just one particular bound, and it 
may not be the best bound for certain applications. Further work may be able to find 
new bounds which can reduce conservatism and therefore give better performance for 
particular problems. 

As noted in Chapter 6, almost all robust control theory available in the literature 
address either the issue of (structured or unstructured) plant uncertainty or the issue 
of uncertain exogenous disturbance rejection. However, the stability of closed-loop 
systems could be very sensitive to numerical errors in the controller gains themselves. 
Hence, an important issue in robust control should be the robustness of compensators 
with respect to implementation errors arising from floating point accuracy in the 
processor. While we presented bounds to deal with controller gain variations for 
linear control problems, extending the results of Chapter 6 to address the robustness of 
nonlinear compensators is of utmost importance and remains essentially an untapped 
area of future research. 
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The general fixed-structure positive real stability multipliers developed for real 
parameter uncertainty in Chapter 8 are structured to handle real block uncertainty, 
but the diagonal structure of Z(s ) is not the most general form for the case of repeated 


uncertainties. For example, if the uncertainty matrix is 

' S l 0 0 O' 

0 <*! 0 0 

0 0 S 2 0 ’ 

0 0 0 < 5 2 


then the commutative property AZ)(s) = D(s) A will hold for a multiplier of the form 


D(s) 


d n (s) d x2 (s) 0 0 

dx 2 (s) d x3 (s) 0 0 

0 0 d 2 i(s) d 22 (s) 

0 0 ^22(5) d 2 3 ( s ) 


Similarly, a generalization of the structure of N(s) can be made which will result 


in a nondiagonal multiplier which will still satisfy the condition N(s)A — A*N(s). 
Further work on this problem can generalize the multiplier structure to handle such 


a case. Furthermore, this structure for the multiplier Z(s ) only allows us to consider 
realizable multipliers. However, this precludes us from considering the famous Popov 
multiplier, Z(s) = I 4- Ns. When using the Popov multiplier, one must assume that 
(I+Ns)G s (s) is realizable, where (? s (s) = M~ l + (I+G e d(s)M x )G ec t(s). Thus, instead 
of assuming that Z(s) is realizable, we can take 


Z(s') — D r (s) — s7V r (s) •+■ sD\ T s^D 2 + ... 4- s T D T — sN\ — s^N 2 — ... — s r N r , 


where D r (s ) and N r (s) are the realizable portions of the multipliers, D x and N x are 
constant matrix scales, and r is chosen such that s r G s (s ) is realizable. 

Finally, although the theory for designing linear output feedback controllers is 
quite mature, nonlinear output feedback controller synthesis remains relatively un- 
developed. In numerous real world applications, system nonlinearities such as sat- 
uration, relay, deadzone, quantization, geometric and material nonlinearities require 



nonlinear output feedback controllers. Furthermore, for linear plants with parametric 
uncertainty and nonquadratic performance criteria, nonlinear controllers exist that 
generate superior performance over the best linear controller. A fruitful area of re- 
search is to develop a fixed-structure controller synthesis framework for nonlinear 
control. The motivation for fixed-structure nonlinear control theory is to address 
controller synthesis within a class of candidate nonlinear feedback controller struc- 
tures. Specifically, control Lyapunov functions can be used to provide a controller 
synthesis framework by assuring global or local asymptotic stability for an a priori 
fixed class of nonlinear feedback controllers. A specific controller within this class 
can now be chosen to optimize a given performance performance functional. Thus, 
this provides a constructive framework where Lyapunov theory is used to guarantee 
global or local asymptotic stability over a class of nonlinear feedback controllers while 
optimization is performed over the free controller gains so as to minimize a specific 
performance functional. 
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Appendix 


Controller Configuration 2 


For this configuration, the controller given by (2.37) can be expressed as 


0 

1 

0 

0 


'O' 


■ — 0.3cc?j 

0 

0 

B r = 

1 

0 

0 

0 

1 

1 

0 

0 

Wl 

3 

1 

— 0.3o>2 


0 


C c = [0 0 hull 0]. 




Note that there are three free parameters, namely, ui, 0 J 2 , and k 3 . Thus 


and 


I< = 


ui 1 0 

0 cj 2 
0 0 


0 

0 , 
h 


K 0 = 


b 1 0 
000 
000 
000 
000 
000 
000 
000 
000 


0 0 

0 0 

1 0 
0 0 
0 0 
0 1 
0 0 
0 0 
0 0 


0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 


0 0 
0 0 
0 0 
0 0 
0 0 , 
0 0 
0 0 
0 0 
0 0 
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Li = 


0 

0 

o' 
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Controller Configuration 3 


For this configuration, the controller given by (2.38) can be expressed as 
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Note that there are three free parameters, namely, wi, u) 2 , and k 3 . Thus 
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Controller Configuration 4 


For this configuration, the controller given by (2.39) can be expressed as 
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Note that there are four free parameters, namely, u } 2 , u > 3, and £3. Thus 
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